微分方程求解

IBB-EUT

我最近在做四阶龙格库塔算法的微分方程求解，求大神分享代码，感谢！
9 `2 k* p( a; W* G5 c

Uifhjvv

function yout = runge_kutta(fun,t,yin,varargin)
t_interval = t(2)-t(1);
! B" F6 A8 ?) f+ b2 Ck1 = t_interval*fun(t(1),yin,varargin);
k2 = t_interval*fun(t(1)+t_interval/2,yin+k1/2varargin);
, @" }% f6 e, n" @. qk3 = t_interval*fun(t(1)+t_interval/2,yin+k2/2,varargin);
; [5 X9 _9 y5 F* H8 yk4 = t_interval*fun(t(1)+t_interval,yin+k3,varargin);
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)
t_interval = t(2)-t(1);
! c1 ^5 D3 K) L( l4 D! uk1 = t_interval*fun(t(1),y ...

谢谢大佬，我先去试试
6 t' E  S" c  I. b6 O7 }

shuddkk

楼上正解，，

cichishia