シミュレーションの演習のコード

GS12.BGSのコードを以下に示します。このコードは BGSP用です。第1行の

BOND SYMBOL EXPR OUT DEFIN PARM INIT

はコードの一部ですから省略してはいけません。また、BGSPではカラム指定が厳格で1カラムでもずれるとエラーになります。また、許される文字は英数字、コロン(:)、スラッシュ(/)、加減乗算記号(+-*)、半角空白のみであり、制御文字が入るとエラーになります。倍角の空白はエラーを生じますが、ソースをプリントしても気づきにくいので注意が必要です。

GS12のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SE1    11
ZJ     -11   33    20
OJ     -33   14    13    19 
II1    14
RE1    13
RE2    20
GY1    19    21
OJ     -21   24    23    29
II2    24
RE3    23 
SE2    29 
/
BC: W24   FL24
BC: M21   EF21
BC: I14   FL14
BC: M29   EF29
BC: I20   FL20
/
SE1
PA: PE1    2.200D+02
       E=PE1
EOD
II1
PA: PI1    1.7500D-01
        L=Z/PI1
EOD
RE1
PA: PR1    8.750D-01
        R=Z/PR1
EOD
RE2
PA: PR2    1.750D-00
        R=Z/PR2
EOD
GY1
PA: PGY1   1.0930D-02
FB: XXR2   FL20
        GY=PGY1*XXR2
EOD
II2
PA: PI2    0.8000D+00
       L=Z/PI2
EOD
RE3
PA:PR3     6.600D-02
       R=Z*PR3
EOD
SE2
PA: PS2    1.000D+02
PA: PS3    2.500D+00
PA: PS4    0.000D+00
       IF(T.GT.PS3)THEN
        E=PS2
       ELSE
        E=PS4
       ENDIF
EOD   
/
0.0000            5.000D-00         1.000D-08         4000
/