四阶龙格库塔法的基本思想_请问用四阶龙格库塔法解二阶微分方程的思想是什么?最近我遇到了一个难题!就是求一个二阶微分方程,形式如：y''=f(y),初始条件是y(0)=0,y'(0)=0,y是t的函数....

默认y的单位是弧度

k=1000;

t=0:0.001:1;

Y=[];

err=1

K=[];

Ymax=[];

xishu=1.01;

while err

X=[0 0];

k=xishu*k;

K=[K;k];

Y=[];

for i=1:1001

Y_1 = Runge_Kutta41(t(i),X,@folded_wing,0.001,k);

Y=[Y;t(i),Y_1];

X=Y_1;

end

ymax=max(Y(:,2));

Ymax=[Ymax;ymax];

if ymax>=pi/2&&ymax<=92/rad2deg(1)

break

elseif ymax>92/rad2deg(1)

xishu=0.98;

else

xishu=1.01;

end

end

plot(Y(:,1),Y(:,2),'-b','linewidth',2);

hold on

grid on

xlabel('t-time','fontsize',14)

ylabel('y(rad)','fontsize',14)

plot([0,1],[pi/2 pi/2],'-g','linewidth&#