test_2.m

来自「提供积分算法变步长Gill法代码」· M 代码 · 共 11 行

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)]);