test_2.m

This the mathematical computational method of step-vary Gill method.

clc;
tic;
[y t] = Gill(@vdp_2,[0 5],[0 0],[.01,1e-2,10]);
Time  = toc;
figure

plot(t,y(1,:),'r-',t,y(2,:),'b:');
xlabel('Time s ');
ylabel('');
legend('y1','y2');
title(['耗时:' num2str(Time)]);