用四阶龙格库塔法求解微分方程的Matlab方法1页.docx

新建一个函数脚本文件t0=0; tN=15; y0=5;0; h=0.01;t = t0 : h : tN; N = length (t); j = 1;for i = 1 : N t1 = t0 + h; K1 = R_K(t0, y0); K2 = R_K(t0 + h/2, y0 + h*K1/2); K3 = R_K(t0 + h/2, y0 + h*K2/2); K4 = R_K(t0 + h, y0 + h*K3); y1 = y0 + (h/6)*(K1 + 2*K2 + 2*K3 + K4); Y1(j)=y1(1); Y2(j)=y1(2); t0=t1; y0=y1; j=j+1;end plot (t, Y1, t, Y2); grid; 新建一个函数M-K.m文件function ydot = f (t, y)ydot=zeros(size(y);