Matlab四阶龙格库塔法求解常微分方程答:2、然后用四阶龙格库塔法求其数值解 figure(2)y0=[1]; %初值y(0)=1 h=0.1;a=0;b=5;[x,y] = runge_kutta(@(x,y)ode_fun(x,y),y0,h,a,b);disp(' x y')A=[x',y']plot(x,y,'LineWidth',1.5),grid on xlabel('x'),ylabel('y(x)');3、绘制y-x的曲线...
用四级四阶龙格-库塔公式求解初值问题答:Nstep=floor(1/h);[x2,y2]=Rung_Kutta(@myfun,x0,y0,h,Nstep);plot(x2,y2)xlabel('x')ylabel('y')box off grid on 函数程序1 function y1=myfun(x,y)a self-defined function y1=sqrt(x+y);end 函数程序2 function [x,y]=Rung_Kutta(fxy,x0,y0,h,Nstep)Rung_Kutta法求解...