【四阶龙格库塔步长求解步长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]