基于MATLAB实现固定步长的四阶龙格库塔方法求解常微分方程初值问题代码实例
代码构造函数function [yn,xn] = RKorder4(f,a,b,h,y0)%a，b为上下界；h为步长；y0为初值。%返回yn、xn向量，并绘制曲线yn=zeros(1,((b-a)/h)+1);yn(1)=y0;xn=a:h:b;for i=1:((b-a)/h) k1=h*f(xn(i),yn(i)); k2=h*f(xn(i)+h/2,yn(i)+k1/2); k3=h*f(xn(i)+h/2,yn(i)+k2/2); k4=h*f(xn(
