comput.m

来自「自己编写的变步长龙格库塔方法解强非线性微分方程」· M 代码 · 共 5 行

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