test_3.m

提供积分算法变步长Gill法代码

clc;
clear;
% 此题为龙格-库塔法算法的应用


A      = 10;
B      = 1;

para.A = A;
para.B = B;
tic;
[y t] = Gill(@vdp_3,[1 9],[1 5],[.01,1e-2,10],para);
Time  = toc;
figure
plot(t,y(1,:),'r-',t,y(2,:),'b:');
xlabel('Time s ');

legend('y1','y2');
title(['耗时:' num2str(Time)]);