范德堡方程的龙格库塔算法求解(共2页).docx

精选优质文档-倾情为你奉上龙格库塔算法：function x,y=lgkt(f,x0,xf,y0,h)h=0.1;n=fix(xf-x0)/h); x(1)=x0;y(:,1)=y0;for i=1:n-1 x(i+1)=x(i)+h; k1=f(x(i),y(:,i); k2=f(x(i)+h/2,y(:,i)+h*k1/2); k3=f(x(i)+h/2,y(:,i)+h*k2/2); k4=f(x(i)+h,y(:,i)+h*k3); y(:,i+1)=y(:,i)+h*(k1+2*k2+2*k3+k4)/6;end函数化简：function f=f(t,y)f=y(2);(1-y(1).2).*y(2)-y(1);Command Window的输入：x0=input( x0=); xf=input( xf=);while x0=xf disp(输入不合逻辑，请重新输入)