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

テストコードg1のコードを以下に示します。このコードは BGSP用です。第１行の 

BOND SYMBOL EXPR OUT DEFIN PARM INIT

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

g1のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SE1    1
OJ     -1    2     3
CI1    2
RE1    3
/
BC: X     EF1
BC: Y     EF2
/
SE1
PA: PE1    100.0D00
PA: PE2    100.0D00
PA: PI     3.141592654
PA: PFRQ1  250.0D00
PA: PFRQ2  10.D00
       E=PE1*sin(2*PI*PFRQ1*T)+PE2*sin(2*PI*PFRQ2*T)
EOD
RE1
PA: PR1    2.0D3
       R=Z/PR1
EOD
CI1
PA: PC1    3.2D-6
       C=Z/PC1
EOD
/
0.00              0.25               0.001             250
/

20-Simによるコード

Sine
parameters
	real PE1 = 100.0 {V};			// amplitude of the wave
	real PE2=100.0 {V};
	real PFRQ1=250.0 {Hz};
	real PFRQ2=10.0 {Hz};
variables 
 real omega1;
real omega2;
equations
	omega1=2*pi*PFRQ1;
	omega2=2*pi*PFRQ2;
	"calculate the sine wave"
	output = PE1 * sin ( omega1 * time)+PE2*sin(omega2*time);

MSE1
variables
	real flow;
equations
	p.e = effort;
	flow = p.f;

CI1
parameters
	real c = 3.2e-06 {F};
equations
	state = int(p.f);
	p.e = state / c;

RE1
parameters
	real r = 2.0e3 {ohm};
equations
	p.e = r * p.f;