一次元静電粒子コードについて・・

一度物理学のところで質問したのですが、誰も答えてくれなかったので(ずうずうしいとは思いますが)今度はこちらで質問させていただきます。
これは、VBで作ってみたPICコードというものなのですが(電界と粒子の運動の関係のコードです)、このコードでは右側の図にSin(正弦)波形が出力されます。この波形を変えたいのですが、どこを変えればいいのでしょう？私は、ポテンシャルの初期値を変えればいいと思い何度もチャレンジしましたが、うまくいきません。　どなたか分かる方教えてください!! （つづきがあります。）

'一次元静電粒子コード
M = 3:IM = 2 ^ (M + 1) + 1:IM2 = IM * 2 - 1
Dim RO(18):Dim PHI(18):Dim A(18, 4):Dim E(18):Dim X(161):Dim V(161)
' X:電子の位置（Ｘ座標） V:電子の速度 RO:電荷密度
' A(I,J):ポアソン・マトリックスエレメント PHI:静電ポテンシャル
' E:静電界（X方向） NP:1セル当りの粒子数 XMAX:全メッシュ数
' KK:最大計算回数 DT:タイムステップ XMAX:X方向の最大値点
' DX: メッシュ間隔
NP = 10:NPT = NP * (IM - 1):PW = 1 / NP
XMAX = IM - 1:DX = XMAX / (IM - 1):DXP = DX / NP:DT = 0.1
' 粒子の初期設定
For IP = 1 To NPT
X(IP) = DXP * (IP - 0.5):V(IP) = 0
Next IP
' ポテンシャルの初期値
For I = 1 To IM
PHI(I) = 5! * Sin(2 * 3.14159 * (I - 1) / (IM - 1))
Next I
' 計算開始
For KK = 1 To 120
' 電場の計算
For I = 1 To IM - 1
E(I) = -(PHI(I + 1) - PHI(I)) / DX
Next I
E(IM) = E(1):E(0) = E(IM - 1):E(IM + 1) = E(2)

つづく・・・

No.1 回答者： OKWebAdmin 回答日時： 2003/01/19 21:17

管理者より：

続きの質問があるのでそちらをご参照下さい

参考URL：http://www.okweb.ne.jp/kotaeru.php3?q=449550