9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
+
+
-
|
# %% Revision 1.1 2000/12/28 18:15:21 peterg
# %% To RCS
# %%
# %% Revision 1.1 1998/03/04 11:45:49 peterg
# %% Initial revision
# %%
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
big = 1e5; # Large number
# Initial states -- needed to choose an approppriate mass
P_0 = 1e5;
V_0 = 1;
T_0 = 300;
# Parameters
c_v = 718.0; # Parameter c_v for CU
gamma_g = 1.4; # Parameter gamma_g for CU
m_g = P_0*V_0/(T_0*(gamma_g-1)*c_v);# Parameter m for CU
## Removed by MTT on Wed Aug 6 17:32:59 BST 2003: pressure = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003
r = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003
|