微分方程求解

IBB-EUT

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

cichishia

shuddkk

楼上正解，，

IBB-EUT

Uifhjvv 发表于 2020-8-13 13:10
: Q8 l9 F9 G- i8 p) B( R1 afunction yout = runge_kutta(fun,t,yin,varargin)
3 ~1 v& s& l' w! W) i% xt_interval = t(2)-t(1);
k1 = t_interval*fun(t(1),y ...

谢谢大佬，我先去试试

Uifhjvv

function yout = runge_kutta(fun,t,yin,varargin)
. t* @) m, b# I5 Yt_interval = t(2)-t(1);
) D! d) A4 M9 P8 s. v' [- D+ Bk1 = t_interval*fun(t(1),yin,varargin);
k2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);
k3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);
k4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);
# D3 @9 X. n# |7 K3 vyout = yin+(k1+2*k2+2*k3+k4)/6;
end