一次元静電粒子コードについて・・(つづき２)
'後退代入法
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
