ey105755的博客

原创 记录一下四阶龙格库塔用于求解一阶微分方程

f为微分方程 dy/dx=f(t,y) t为自变量，y为因变量。调用主程序 并且以 dy/dx=-2*t*y为例子。%定义具体的微分方程。%tspan为时间区间[t_strat,t_end]%y0为初始值，即t=t_strat对应的值。%t：计算在步长h下t区间内的点数组。