comput.m

自己编写的变步长龙格库塔方法解强非线性微分方程

function [y1,y2]=comput(t,y,h)
y1=rk4(t,y,h);
y21=rk4(t,y,h/2);
y2=rk4(t+h/2,y21,h/2);