meili 发表于 2022-10-27 15:48:57

【四阶龙格库塔步长求解步长h怎么确定?选取多少结果才比较精确?】

<p>问题:【四阶龙格库塔步长求解步长h怎么确定?选取多少结果才比较精确?】
<p>答案:↓↓↓<p class="nav-title mt10" style="border-top:1px solid #ccc;padding-top: 10px;">江小芳的回答:<div class="content-b">网友采纳  假设求解初值问题:  y'=y-2x/y(0<x<1)  y(0)=1  设步长h=0.2,从x=0直到x=1用四阶龙格库塔法:  PrivateSubForm_click()  DimxAsSingle,yAsSingle  Dimk1AsSingle,k2AsSingle,k3AsSingle,k4AsSingle  Dimy1AsSingle,y2AsSingle,y3AsSingle,y4AsSingle  x=0  y=1  h=0.2  Fori=0To1Step0.2  k1=f(x,y)  y1=y+h/2*k1  x=x+h/2  k2=f(x,y1)  y2=y+h/2*k2  k3=f(x,y2)  y3=y+h*k3  x=x+h/2  k4=f(x,y3)  y=y+h/6*(k1+2*k2+2*k3+k4)  Print"x=";x;"y=";y  Next  EndSub  PrivateFunctionf(aAsSingle,bAsSingle)AsSingle  f=b-2*a/b  EndFunction
页: [1]
查看完整版本: 【四阶龙格库塔步长求解步长h怎么确定?选取多少结果才比较精确?】