微分方程求解

IBB-EUT

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

Uifhjvv

function yout = runge_kutta(fun,t,yin,varargin)
* l: H  o  m2 N& x) Vt_interval = t(2)-t(1);
$ J" Y, A5 c" fk1 = 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);
, K" t, v% g) ?+ M0 g% V( K* W8 S8 kyout = yin+(k1+2*k2+2*k3+k4)/6;
end

IBB-EUT

Uifhjvv 发表于 2020-8-13 13:10
6 ]; `( m6 x$ Qfunction yout = runge_kutta(fun,t,yin,varargin)
9 k# Y" d" W! |' _" C( |; s/ h: tt_interval = t(2)-t(1);
! b! M, {9 u5 Y% ik1 = t_interval*fun(t(1),y ...

4 c; L3 B. @6 D! }谢谢大佬，我先去试试

shuddkk

楼上正解，，

cichishia