matlab龙格库塔法求解微分方程 怎么编程
function [Y] = RK45(t,X,f,h)K1=f(t,X);K2=f(t+h/2,X+h/2*K1) ;K3=f(t+h/2,X+h/2*K2) ;K4=f(t+h,X+h*K3) ;Y=X+h/6*(K1+2*K2+2*K3+K4) ;end以上是4阶龙格库塔法的代码:自己写函数，存为f. mfunction dxdt = f (t,x)dxdt(1)=exp(x(1)*sin(t))+x(2) ;dxdt(2)=exp(x(2)*cos(t))+x(1) ; %...