●EXCEL・VBAの教えてスレ Part2●

●EXCEL・VBAの教えて ..

■コピペモード
□スレを通常表示
□オプションモード
□このｽﾚｯﾄﾞのURL
■項目テキスト

4:どうかご慈悲を
 07/01/09 17:29:47 
 ｄσｖ/ｄｘ　+４(ｋσｖtanφ+Ｃｗ)/Ｄ＝γω＾２*（Ｒo+x） 
 
変数はｘとσｖでこれらをＥｘｃｅｌＶＢＡによって倍精度ルンゲクッタ法で解くプログラムはどのようにしたらいいですか？以下の修正部分を教えてください。お願いします。物性値の依存性は無視してプログラムとしてどう修正すればよいか教えてください。 
Sub 円筒部() 
'Diff. Eq: F(X,Y)= ω^2*(Ro+X)*649.96*Y^0.0959-0.1*(1.1*(1-SIN(38.727*Y^(-0.0204)))*Y *TAN(28.638* Y ^(-0.04))+0.239* Y ^0.5389) 
'READ INITAIAL CONDITIONS 
    X = Cells(4, 1).Value: Y = Cells(4, 2).Value 
    H = Cells(4, 3).Value: N = Cells(4, 4).Value 
    L = Cells(4, 5).Value: ω = Cells(4, 6).Value 
    Ro = Cells(4, 7).Value ' 
    For I = 0 To N 
        Cells(I + 7, 2).Value = X 
        Cells(I + 7, 3).Value = Y 
        K1 = H * (ω ^ 2 * (Ro + X) * 649.96 * Y ^ 0.0959 - 0.1 * (1.1 * (1 - Sin(38.727 * Y ^ (-0.0204))) * Y * Tan(28.638 * Y ^ (-0.04)) + 0.239 * Y ^ 0.5389)) 
        K2 = H * (ω ^ 2 * (Ro + X) * 649.96 * (Y + K1 / 2) ^ 0.0959 - 0.1 * (1.1 * (1 - Sin(38.727 * (Y + K1 / 2) ^ (-0.0204))) * (Y + K1 / 2) * Tan(28.638 * (Y + K1 / 2) ^ (-0.04)) + 0.239 * (Y + K1 / 2) ^ 0.5389)) 
        K3 = H * (ω ^ 2 * (Ro + X) * 649.96 * (Y + K2 / 2) ^ 0.0959 - 0.1 * (1.1 * (1 - Sin(38.727 * (Y + K2 / 2) ^ (-0.0204))) * (Y + K2 / 2) * Tan(28.638 * (Y + K2 / 2) ^ (-0.04)) + 0.239 * (Y + K2 / 2) ^ 0.5389)) 
        K4 = H * (ω ^ 2 * (Ro + X) * 649.96 * (Y + K3) ^ 0.0959 - 0.1 * (1.1 * (1 - Sin(38.727 * (Y + K3) ^ (-0.0204))) * (Y + K3) * Tan(28.638 * (Y + K3) ^ (-0.04)) + 0.239 * (Y + K3) ^ 0.5389)) 
 
        X = X + H 
        Y = Y + (K1 + 2 * (K2 + K3) + K4) / 6 
    Next I 
 
End Sub

次ページ