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

2003/01/19 17:09

このQ&Aのポイント

'後退代入法A(N, NN2) = A(N, NN2) / A(N, NNP1)For I = N To 2 Step -1IMNN = I - NN:If IMNN < 1 Then IMNN = 1For L = I - 1 To IMNN Step -1A(L, NN2) = A(L, NN2) - A(I, NN2) * A(L, NNP1 + I - L)Next LNext IA(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) * 100Picture1.Circle (AX, AY), 1, vbBlueNext I'電界の初期値AXO = 350:AYO = 200 + E(1) * 100Picture1.Circle (AXO, AYO), 10, vbGreen'電子の位置（格子点）に対する電界の値 (RIGHT SIDE)For I = 2 To IM - 1AX = 350 + (I - 1) * 200 / (IM - 1):AY = 200 + E(I) * 100Picture1.Circle (AX, AY), 1, vbGreenPicture1.Line (AXO, AYO)-(AX, AY), vbRed;AXO = AX :AYO = AYNext INext KKEnd Sub以上です。長ったらしく、また見にくくて申し訳ありません。分かる方がいましたら、ぜひ教えてください。
一次元静電粒子コードについてのコード解説です。具体的には、後退代入法による演算やグラフの描画方法、電子の位置と速度の関係、電界の初期値の設定、電子の位置に対する電界の値の計算などが説明されています。また、コードの長さや見やすさについても謝罪があります。分かる方は教えていただきたい旨を伝えています。
一次元静電粒子コードについてのコード解説です。具体的なアルゴリズムや計算方法、グラフの描画方法などが説明されています。コードは一部省略されていますが、全体の流れや主要な処理が分かるようになっています。コードの見やすさについては謝罪があります。分かる方に教えていただきたい旨を伝えています。

corgi403
お礼率34% (26/75)

Visual Basic
回答数1
ありがとう数0

みんなの回答 （1）
専門家の回答

質問者が選んだベストアンサー

ベストアンサー

jmh
ベストアンサー率23% (71/304)

2003/01/20 04:17 回答No.1

最初の「右側の図に…出力されます」と最後の方の「　　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() が何に依存するか分かると思います。