微分方程求解

IBB-EUT

我最近在做四阶龙格库塔算法的微分方程求解，求大神分享代码，感谢！

Uifhjvv

function yout = runge_kutta(fun,t,yin,varargin)
) p* W4 }' b. R) z9 K9 Ft_interval = t(2)-t(1);
k1 = t_interval*fun(t(1),yin,varargin);
5 }, H' m) l' O$ T$ rk2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);
5 L- Y7 L0 i4 ak3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);
k4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);
* j! Z, `; ]" r# I; i  {yout = yin+(k1+2*k2+2*k3+k4)/6;
end

IBB-EUT

Uifhjvv 发表于 2020-8-13 13:10
function yout = runge_kutta(fun,t,yin,varargin)
; G" k, I$ j# J0 L; \, }, B+ Dt_interval = t(2)-t(1);
k1 = t_interval*fun(t(1),y ...

8 H6 j7 c4 @" q! p1 y% }谢谢大佬，我先去试试

shuddkk

楼上正解，，

cichishia