微分方程求解

IBB-EUT

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

Uifhjvv

function yout = runge_kutta(fun,t,yin,varargin)
t_interval = t(2)-t(1);
k1 = t_interval*fun(t(1),yin,varargin);
% i3 C4 d: y& l9 ]4 Z- Zk2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);
) G' K0 ]' N1 p5 q5 Fk3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);
k4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);
. g8 l2 a, B! O& ?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)
* j6 c* p& R7 W9 H& mt_interval = t(2)-t(1);
& q/ P. K3 V$ ]) k, o2 uk1 = t_interval*fun(t(1),y ...

2 N( z+ j/ d. @* \% P9 R谢谢大佬，我先去试试
8 P9 _4 i3 M/ s% n2 S

shuddkk

楼上正解，，

cichishia