一次元静電粒子コードについて・・(つづき２)

'後退代入法
A(N, NN2) = A(N, NN2) / A(N, NNP1)
For I = N To 2 Step -1
IMNN = I - NN:If IMNN < 1 Then IMNN = 1
For L = I - 1 To IMNN Step -1
A(L, NN2) = A(L, NN2) - A(I, NN2) * A(L, NNP1 + I - L)
Next L
Next I
A(1,NN2)=A(1,NN2)/A(1,NNP1)
'グラフ
Picture1.BackColor = RGB(255, 255, 255)
Picture1.Line (50, 100)-(50, 300):Picture1.Line (50, 200)-(250, 200)
Picture1.Line (350, 100)-(350, 300)
Picture1.Line (350, 200)-(550, 200)
'電子の位置AXに対する速度AYの変化　　(LEFT SIDE)
For I = 1 To NPT 'NPT=NP*(IM-1):全メッシュ数=160個
AX = X(I) * 200 / (IM - 1) + 50 :AY = 200 - V(I) * 100
Picture1.Circle (AX, AY), 1, vbBlue
Next I
'電界の初期値
AXO = 350:AYO = 200 + E(1) * 100
Picture1.Circle (AXO, AYO), 10, vbGreen
'電子の位置（格子点）に対する電界の値 (RIGHT SIDE)
For I = 2 To IM - 1
AX = 350 + (I - 1) * 200 / (IM - 1):AY = 200 + E(I) * 100
Picture1.Circle (AX, AY), 1, vbGreen
Picture1.Line (AXO, AYO)-(AX, AY), vbRed;AXO = AX :AYO = AY
Next I
Next KK
End Sub

以上です。
長ったらしく、また見にくくて申し訳ありません。　分かる方がいましたら、ぜひ教えてください。

関連URL:http://www.okweb.ne.jp/kotaeru.php3?q=449543
関連URL:http://www.okweb.ne.jp/kotaeru.php3?q=449550

No.1 ベストアンサー 回答者： jmh 回答日時： 2003/01/20 04:17

最初の「右側の図に…出力されます」と最後の方の「

　　AXO = 350:AYO = 200 + E(1) * 100
　　Picture1.Circle (AXO, AYO), 10, vbGreen
　　'…(RIGHT SIDE)
　　For I = 2 To IM - 1
　　　　AX = 350 + (I - 1) * 200 / (IM - 1):AY = 200 + E(I) * 100
　　　　Picture1.Circle (AX, AY), 1, vbGreen
　　　　Picture1.Line (AXO, AYO)-(AX, AY), vbRed:AXO = AX :AYO = AY
　　Next I
」
だけから考えて、IM, E()の値を変えればいいと思います。「IM =」と「E(*) =」を検索（パターンマッチング）したら、IM と E() が何に依存するか分かると思います。