1doi: "10.1103/PhysRevLett.103.263201";
2reference: "K Pernal, R Podeszwa, K Patkowski, K Szalewicz, Phys.Rev.Lett. 103, 263201 (2009) ";
3F000001:   0.0000000000000000E+00;
4F000002:   0.0000000000000000E+00;
5F000003:   0.0000000000000000E+00;
6F000004:   0.0000000000000000E+00;
7F000005:   0.0000000000000000E+00;
8F000006:   0.0000000000000000E+00;
9F000007:   0.0000000000000000E+00;
10F000008:   0.0000000000000000E+00;
11F000009:   0.0000000000000000E+00;
12F000010:   0.0000000000000000E+00;
13F000011:   0.0000000000000000E+00;
14F000012:   0.0000000000000000E+00;
15F000013:   0.0000000000000000E+00;
16/* comp: rhoa.le.0.5000000000000000E-18
17*/
18F000014(taua):= 0.2000000000000000E+01*taua;
19/* comp: F000014(taua).gt.0.1000000000000000E-17
20*/
21/* comp: rhoa.le.0.1000000000000000E-17
22*/
23/* comp: F000014(taua).le.0.1000000000000000E-17
24*/
25F000015(rhoa):= 0.2387324146378430E+00/rhoa;
26F000016(rhoa):= F000015(rhoa)^(1/3);
27F000017:   0.1000000000000000E+01;
28F000018:  -0.2000000000000000E+01;
29F000019: 0.1000000000000000E+01-F000017;
30F000020: 0.1000000000000000E+01+F000017;
31F000021: F000019^0.2000000000000000E+01;
32F000022: F000020^0.2000000000000000E+01;
33/* comp: F000019.gt.0.1000000000000000E-17
34*/
35/* comp: F000020.gt.0.1000000000000000E-17
36*/
37F000023: F000020^(1/3);
38F000024: F000020*F000023;
39F000025: F000018+F000024;
40F000026: F000025*0.1923661050931536E+01;
41F000027(rhoa):= 0.2137000000000000E+00*F000016(rhoa);
42F000028(rhoa):= 0.1000000000000000E+01+F000027(rhoa);
43F000029(rhoa):= 0.6218140000000000E-01*F000028(rhoa);
44F000030(rhoa):= -(F000029(rhoa));
45F000031(rhoa):= sqrt(F000016(rhoa));
46F000032(rhoa):= F000016(rhoa)*F000031(rhoa);
47F000033(rhoa):= 0.7595700000000000E+01*F000031(rhoa);
48F000034(rhoa):= 0.3587600000000000E+01*F000016(rhoa);
49F000035(rhoa):= F000033(rhoa)+F000034(rhoa);
50F000036(rhoa):= 0.1638200000000000E+01*F000032(rhoa);
51F000037(rhoa):= F000035(rhoa)+F000036(rhoa);
52F000038(rhoa):= 0.4929400000000000E+00*F000016(rhoa);
53F000039(rhoa):= F000038(rhoa)*F000016(rhoa);
54F000040(rhoa):= F000037(rhoa)+F000039(rhoa);
55F000041(rhoa):= 0.6218140000000000E-01*F000040(rhoa);
56F000042(rhoa):= 0.1000000000000000E+01/F000041(rhoa);
57F000043(rhoa):= 0.1000000000000000E+01+F000042(rhoa);
58F000044(rhoa):= log(F000043(rhoa));
59F000045(rhoa):= F000030(rhoa)*F000044(rhoa);
60F000046(rhoa):= 0.2054800000000000E+00*F000016(rhoa);
61F000047(rhoa):= 0.1000000000000000E+01+F000046(rhoa);
62F000048(rhoa):= 0.3109070000000000E-01*F000047(rhoa);
63F000049(rhoa):= -(F000048(rhoa));
64F000050(rhoa):= sqrt(F000016(rhoa));
65F000051(rhoa):= F000016(rhoa)*F000050(rhoa);
66F000052(rhoa):= 0.1411890000000000E+02*F000050(rhoa);
67F000053(rhoa):= 0.6197700000000000E+01*F000016(rhoa);
68F000054(rhoa):= F000052(rhoa)+F000053(rhoa);
69F000055(rhoa):= 0.3366200000000000E+01*F000051(rhoa);
70F000056(rhoa):= F000054(rhoa)+F000055(rhoa);
71F000057(rhoa):= 0.6251700000000000E+00*F000016(rhoa);
72F000058(rhoa):= F000057(rhoa)*F000016(rhoa);
73F000059(rhoa):= F000056(rhoa)+F000058(rhoa);
74F000060(rhoa):= 0.3109070000000000E-01*F000059(rhoa);
75F000061(rhoa):= 0.1000000000000000E+01/F000060(rhoa);
76F000062(rhoa):= 0.1000000000000000E+01+F000061(rhoa);
77F000063(rhoa):= log(F000062(rhoa));
78F000064(rhoa):= F000049(rhoa)*F000063(rhoa);
79F000065(rhoa):= 0.1112500000000000E+00*F000016(rhoa);
80F000066(rhoa):= 0.1000000000000000E+01+F000065(rhoa);
81F000067(rhoa):= 0.3377372788077926E-01*F000066(rhoa);
82F000068(rhoa):= -(F000067(rhoa));
83F000069(rhoa):= sqrt(F000016(rhoa));
84F000070(rhoa):= F000016(rhoa)*F000069(rhoa);
85F000071(rhoa):= 0.1035700000000000E+02*F000069(rhoa);
86F000072(rhoa):= 0.3623100000000000E+01*F000016(rhoa);
87F000073(rhoa):= F000071(rhoa)+F000072(rhoa);
88F000074(rhoa):= 0.8802600000000000E+00*F000070(rhoa);
89F000075(rhoa):= F000073(rhoa)+F000074(rhoa);
90F000076(rhoa):= 0.4967100000000000E+00*F000016(rhoa);
91F000077(rhoa):= F000076(rhoa)*F000016(rhoa);
92F000078(rhoa):= F000075(rhoa)+F000077(rhoa);
93F000079(rhoa):= 0.3377372788077926E-01*F000078(rhoa);
94F000080(rhoa):= 0.1000000000000000E+01/F000079(rhoa);
95F000081(rhoa):= 0.1000000000000000E+01+F000080(rhoa);
96F000082(rhoa):= log(F000081(rhoa));
97F000083(rhoa):= F000068(rhoa)*F000082(rhoa);
98F000084: F000017*F000017;
99F000085: F000017*F000084;
100F000086: F000017*F000085;
101F000087: F000026*F000086;
102F000088: F000026-F000087;
103F000089: 0.5848223622634647E+00*F000088;
104F000090: 0.1000000000000000E+01-F000087;
105F000091(rhoa):= F000045(rhoa)*F000090;
106F000092(rhoa):= F000064(rhoa)*F000087;
107F000093(rhoa):= F000091(rhoa)+F000092(rhoa);
108F000094(rhoa):= F000083(rhoa)*F000089;
109F000095(rhoa):= F000093(rhoa)-F000094(rhoa);
110F000096(rhoa):= rhoa*F000095(rhoa);
111F000097(sigmaaa):= 0.2500000000000000E+00*sigmaaa;
112F000098(rhoa,sigmaaa):= F000097(sigmaaa)/rhoa;
113F000099(rhoa,sigmaaa,taua):= F000014(taua)-F000098(rhoa,sigmaaa);
114F000100(rhoa):= rhoa^(8/3);
115F000101(rhoa,sigmaaa):= sigmaaa/F000100(rhoa);
116F000102(rhoa,sigmaaa):= 0.6000000000000000E-01*F000101(rhoa,sigmaaa);
117F000103(rhoa,sigmaaa):= 0.6000000000000000E-01*F000101(rhoa,sigmaaa);
118F000104(rhoa,sigmaaa):= 0.1000000000000000E+01+F000103(rhoa,sigmaaa);
119F000105(rhoa,sigmaaa):= F000102(rhoa,sigmaaa)/F000104(rhoa,sigmaaa);
120F000106(rhoa,sigmaaa):= F000105(rhoa,sigmaaa)*0.0000000000000000E+00;
121F000107(rhoa,sigmaaa):= 0.0000000000000000E+00+F000106(rhoa,sigmaaa);
122F000108(rhoa,sigmaaa):= F000105(rhoa,sigmaaa)*F000107(rhoa,sigmaaa);
123F000109(rhoa,sigmaaa):= 0.2223379300000000E+01+F000108(rhoa,sigmaaa);
124F000110(rhoa,sigmaaa):= F000105(rhoa,sigmaaa)*F000109(rhoa,sigmaaa);
125F000111(rhoa,sigmaaa):= -0.2596089700000000E+01+F000110(rhoa,sigmaaa);
126F000112(rhoa,sigmaaa):= F000105(rhoa,sigmaaa)*F000111(rhoa,sigmaaa);
127F000113(rhoa,sigmaaa):= 0.1000000000000000E+01+F000112(rhoa,sigmaaa);
128F000114(rhoa,sigmaaa,taua):= F000099(rhoa,sigmaaa,taua)/F000014(taua);
129F000115(rhoa,sigmaaa,taua):= F000114(rhoa,sigmaaa,taua)*F000113(rhoa,sigmaaa);
130F000116(rhoa,sigmaaa,taua):= F000115(rhoa,sigmaaa,taua)*F000096(rhoa);
131F000117(rhoa,sigmaaa,taua):= F000116(rhoa,sigmaaa,taua)*0.1000000000000000E+01;
132F000118(rhoa,sigmaaa,taua):= F000117(rhoa,sigmaaa,taua)*0.1000000000000000E+01;
133F000119(rhoa,sigmaaa,taua):= F000001+F000118(rhoa,sigmaaa,taua);
134/* comp: rhob.le.0.5000000000000000E-18
135*/
136F000120(taub):= 0.2000000000000000E+01*taub;
137/* comp: F000120(taub).gt.0.1000000000000000E-17
138*/
139/* comp: rhob.le.0.1000000000000000E-17
140*/
141/* comp: F000120(taub).le.0.1000000000000000E-17
142*/
143F000121(rhob):= 0.2387324146378430E+00/rhob;
144F000122(rhob):= F000121(rhob)^(1/3);
145F000123:   0.1000000000000000E+01;
146F000124:  -0.2000000000000000E+01;
147F000125: 0.1000000000000000E+01-F000123;
148F000126: 0.1000000000000000E+01+F000123;
149F000127: F000125^0.2000000000000000E+01;
150F000128: F000126^0.2000000000000000E+01;
151/* comp: F000125.gt.0.1000000000000000E-17
152*/
153/* comp: F000126.gt.0.1000000000000000E-17
154*/
155F000129: F000126^(1/3);
156F000130: F000126*F000129;
157F000131: F000124+F000130;
158F000132: F000131*0.1923661050931536E+01;
159F000133(rhob):= 0.2137000000000000E+00*F000122(rhob);
160F000134(rhob):= 0.1000000000000000E+01+F000133(rhob);
161F000135(rhob):= 0.6218140000000000E-01*F000134(rhob);
162F000136(rhob):= -(F000135(rhob));
163F000137(rhob):= sqrt(F000122(rhob));
164F000138(rhob):= F000122(rhob)*F000137(rhob);
165F000139(rhob):= 0.7595700000000000E+01*F000137(rhob);
166F000140(rhob):= 0.3587600000000000E+01*F000122(rhob);
167F000141(rhob):= F000139(rhob)+F000140(rhob);
168F000142(rhob):= 0.1638200000000000E+01*F000138(rhob);
169F000143(rhob):= F000141(rhob)+F000142(rhob);
170F000144(rhob):= 0.4929400000000000E+00*F000122(rhob);
171F000145(rhob):= F000144(rhob)*F000122(rhob);
172F000146(rhob):= F000143(rhob)+F000145(rhob);
173F000147(rhob):= 0.6218140000000000E-01*F000146(rhob);
174F000148(rhob):= 0.1000000000000000E+01/F000147(rhob);
175F000149(rhob):= 0.1000000000000000E+01+F000148(rhob);
176F000150(rhob):= log(F000149(rhob));
177F000151(rhob):= F000136(rhob)*F000150(rhob);
178F000152(rhob):= 0.2054800000000000E+00*F000122(rhob);
179F000153(rhob):= 0.1000000000000000E+01+F000152(rhob);
180F000154(rhob):= 0.3109070000000000E-01*F000153(rhob);
181F000155(rhob):= -(F000154(rhob));
182F000156(rhob):= sqrt(F000122(rhob));
183F000157(rhob):= F000122(rhob)*F000156(rhob);
184F000158(rhob):= 0.1411890000000000E+02*F000156(rhob);
185F000159(rhob):= 0.6197700000000000E+01*F000122(rhob);
186F000160(rhob):= F000158(rhob)+F000159(rhob);
187F000161(rhob):= 0.3366200000000000E+01*F000157(rhob);
188F000162(rhob):= F000160(rhob)+F000161(rhob);
189F000163(rhob):= 0.6251700000000000E+00*F000122(rhob);
190F000164(rhob):= F000163(rhob)*F000122(rhob);
191F000165(rhob):= F000162(rhob)+F000164(rhob);
192F000166(rhob):= 0.3109070000000000E-01*F000165(rhob);
193F000167(rhob):= 0.1000000000000000E+01/F000166(rhob);
194F000168(rhob):= 0.1000000000000000E+01+F000167(rhob);
195F000169(rhob):= log(F000168(rhob));
196F000170(rhob):= F000155(rhob)*F000169(rhob);
197F000171(rhob):= 0.1112500000000000E+00*F000122(rhob);
198F000172(rhob):= 0.1000000000000000E+01+F000171(rhob);
199F000173(rhob):= 0.3377372788077926E-01*F000172(rhob);
200F000174(rhob):= -(F000173(rhob));
201F000175(rhob):= sqrt(F000122(rhob));
202F000176(rhob):= F000122(rhob)*F000175(rhob);
203F000177(rhob):= 0.1035700000000000E+02*F000175(rhob);
204F000178(rhob):= 0.3623100000000000E+01*F000122(rhob);
205F000179(rhob):= F000177(rhob)+F000178(rhob);
206F000180(rhob):= 0.8802600000000000E+00*F000176(rhob);
207F000181(rhob):= F000179(rhob)+F000180(rhob);
208F000182(rhob):= 0.4967100000000000E+00*F000122(rhob);
209F000183(rhob):= F000182(rhob)*F000122(rhob);
210F000184(rhob):= F000181(rhob)+F000183(rhob);
211F000185(rhob):= 0.3377372788077926E-01*F000184(rhob);
212F000186(rhob):= 0.1000000000000000E+01/F000185(rhob);
213F000187(rhob):= 0.1000000000000000E+01+F000186(rhob);
214F000188(rhob):= log(F000187(rhob));
215F000189(rhob):= F000174(rhob)*F000188(rhob);
216F000190: F000123*F000123;
217F000191: F000123*F000190;
218F000192: F000123*F000191;
219F000193: F000132*F000192;
220F000194: F000132-F000193;
221F000195: 0.5848223622634647E+00*F000194;
222F000196: 0.1000000000000000E+01-F000193;
223F000197(rhob):= F000151(rhob)*F000196;
224F000198(rhob):= F000170(rhob)*F000193;
225F000199(rhob):= F000197(rhob)+F000198(rhob);
226F000200(rhob):= F000189(rhob)*F000195;
227F000201(rhob):= F000199(rhob)-F000200(rhob);
228F000202(rhob):= rhob*F000201(rhob);
229F000203(sigmabb):= 0.2500000000000000E+00*sigmabb;
230F000204(rhob,sigmabb):= F000203(sigmabb)/rhob;
231F000205(rhob,sigmabb,taub):= F000120(taub)-F000204(rhob,sigmabb);
232F000206(rhob):= rhob^(8/3);
233F000207(rhob,sigmabb):= sigmabb/F000206(rhob);
234F000208(rhob,sigmabb):= 0.6000000000000000E-01*F000207(rhob,sigmabb);
235F000209(rhob,sigmabb):= 0.6000000000000000E-01*F000207(rhob,sigmabb);
236F000210(rhob,sigmabb):= 0.1000000000000000E+01+F000209(rhob,sigmabb);
237F000211(rhob,sigmabb):= F000208(rhob,sigmabb)/F000210(rhob,sigmabb);
238F000212(rhob,sigmabb):= F000211(rhob,sigmabb)*0.0000000000000000E+00;
239F000213(rhob,sigmabb):= 0.0000000000000000E+00+F000212(rhob,sigmabb);
240F000214(rhob,sigmabb):= F000211(rhob,sigmabb)*F000213(rhob,sigmabb);
241F000215(rhob,sigmabb):= 0.2223379300000000E+01+F000214(rhob,sigmabb);
242F000216(rhob,sigmabb):= F000211(rhob,sigmabb)*F000215(rhob,sigmabb);
243F000217(rhob,sigmabb):= -0.2596089700000000E+01+F000216(rhob,sigmabb);
244F000218(rhob,sigmabb):= F000211(rhob,sigmabb)*F000217(rhob,sigmabb);
245F000219(rhob,sigmabb):= 0.1000000000000000E+01+F000218(rhob,sigmabb);
246F000220(rhob,sigmabb,taub):= F000205(rhob,sigmabb,taub)/F000120(taub);
247F000221(rhob,sigmabb,taub):= F000220(rhob,sigmabb,taub)*F000219(rhob,sigmabb);
248F000222(rhob,sigmabb,taub):= F000221(rhob,sigmabb,taub)*F000202(rhob);
249F000223(rhob,sigmabb,taub):= F000222(rhob,sigmabb,taub)*0.1000000000000000E+01;
250F000224(rhob,sigmabb,taub):= F000223(rhob,sigmabb,taub)*0.1000000000000000E+01;
251F000225(rhoa,rhob,sigmaaa,sigmabb,taua,taub):= F000119(rhoa,sigmaaa,taua)+F000224(rhob,sigmabb,taub);
252F000226(rhoa,rhob):= rhoa+rhob;
253/* comp: rhoa.gt.0.5000000000000000E-18
254*/
255/* comp: rhob.gt.0.5000000000000000E-18
256*/
257F000227(rhoa,rhob):= 0.2387324146378430E+00/F000226(rhoa,rhob);
258F000228(rhoa,rhob):= F000227(rhoa,rhob)^(1/3);
259F000229(rhoa,rhob):= rhoa-rhob;
260F000230(rhoa,rhob):= F000229(rhoa,rhob)/F000226(rhoa,rhob);
261F000231:  -0.2000000000000000E+01;
262F000232(rhoa,rhob):= 0.1000000000000000E+01-F000230(rhoa,rhob);
263F000233(rhoa,rhob):= 0.1000000000000000E+01+F000230(rhoa,rhob);
264F000234(rhoa,rhob):= F000232(rhoa,rhob)^0.2000000000000000E+01;
265F000235(rhoa,rhob):= F000233(rhoa,rhob)^0.2000000000000000E+01;
266/* comp: F000232(rhoa,rhob).gt.0.1000000000000000E-17
267*/
268F000236(rhoa,rhob):= F000232(rhoa,rhob)^(1/3);
269F000237(rhoa,rhob):= F000232(rhoa,rhob)*F000236(rhoa,rhob);
270F000238(rhoa,rhob):= F000231+F000237(rhoa,rhob);
271/* comp: F000233(rhoa,rhob).gt.0.1000000000000000E-17
272*/
273F000239(rhoa,rhob):= F000233(rhoa,rhob)^(1/3);
274F000240(rhoa,rhob):= F000233(rhoa,rhob)*F000239(rhoa,rhob);
275F000241(rhoa,rhob):= F000238(rhoa,rhob)+F000240(rhoa,rhob);
276F000242(rhoa,rhob):= F000241(rhoa,rhob)*0.1923661050931536E+01;
277F000243(rhoa,rhob):= 0.2137000000000000E+00*F000228(rhoa,rhob);
278F000244(rhoa,rhob):= 0.1000000000000000E+01+F000243(rhoa,rhob);
279F000245(rhoa,rhob):= 0.6218140000000000E-01*F000244(rhoa,rhob);
280F000246(rhoa,rhob):= -(F000245(rhoa,rhob));
281F000247(rhoa,rhob):= sqrt(F000228(rhoa,rhob));
282F000248(rhoa,rhob):= F000228(rhoa,rhob)*F000247(rhoa,rhob);
283F000249(rhoa,rhob):= 0.7595700000000000E+01*F000247(rhoa,rhob);
284F000250(rhoa,rhob):= 0.3587600000000000E+01*F000228(rhoa,rhob);
285F000251(rhoa,rhob):= F000249(rhoa,rhob)+F000250(rhoa,rhob);
286F000252(rhoa,rhob):= 0.1638200000000000E+01*F000248(rhoa,rhob);
287F000253(rhoa,rhob):= F000251(rhoa,rhob)+F000252(rhoa,rhob);
288F000254(rhoa,rhob):= 0.4929400000000000E+00*F000228(rhoa,rhob);
289F000255(rhoa,rhob):= F000254(rhoa,rhob)*F000228(rhoa,rhob);
290F000256(rhoa,rhob):= F000253(rhoa,rhob)+F000255(rhoa,rhob);
291F000257(rhoa,rhob):= 0.6218140000000000E-01*F000256(rhoa,rhob);
292F000258(rhoa,rhob):= 0.1000000000000000E+01/F000257(rhoa,rhob);
293F000259(rhoa,rhob):= 0.1000000000000000E+01+F000258(rhoa,rhob);
294F000260(rhoa,rhob):= log(F000259(rhoa,rhob));
295F000261(rhoa,rhob):= F000246(rhoa,rhob)*F000260(rhoa,rhob);
296F000262(rhoa,rhob):= 0.2054800000000000E+00*F000228(rhoa,rhob);
297F000263(rhoa,rhob):= 0.1000000000000000E+01+F000262(rhoa,rhob);
298F000264(rhoa,rhob):= 0.3109070000000000E-01*F000263(rhoa,rhob);
299F000265(rhoa,rhob):= -(F000264(rhoa,rhob));
300F000266(rhoa,rhob):= sqrt(F000228(rhoa,rhob));
301F000267(rhoa,rhob):= F000228(rhoa,rhob)*F000266(rhoa,rhob);
302F000268(rhoa,rhob):= 0.1411890000000000E+02*F000266(rhoa,rhob);
303F000269(rhoa,rhob):= 0.6197700000000000E+01*F000228(rhoa,rhob);
304F000270(rhoa,rhob):= F000268(rhoa,rhob)+F000269(rhoa,rhob);
305F000271(rhoa,rhob):= 0.3366200000000000E+01*F000267(rhoa,rhob);
306F000272(rhoa,rhob):= F000270(rhoa,rhob)+F000271(rhoa,rhob);
307F000273(rhoa,rhob):= 0.6251700000000000E+00*F000228(rhoa,rhob);
308F000274(rhoa,rhob):= F000273(rhoa,rhob)*F000228(rhoa,rhob);
309F000275(rhoa,rhob):= F000272(rhoa,rhob)+F000274(rhoa,rhob);
310F000276(rhoa,rhob):= 0.3109070000000000E-01*F000275(rhoa,rhob);
311F000277(rhoa,rhob):= 0.1000000000000000E+01/F000276(rhoa,rhob);
312F000278(rhoa,rhob):= 0.1000000000000000E+01+F000277(rhoa,rhob);
313F000279(rhoa,rhob):= log(F000278(rhoa,rhob));
314F000280(rhoa,rhob):= F000265(rhoa,rhob)*F000279(rhoa,rhob);
315F000281(rhoa,rhob):= 0.1112500000000000E+00*F000228(rhoa,rhob);
316F000282(rhoa,rhob):= 0.1000000000000000E+01+F000281(rhoa,rhob);
317F000283(rhoa,rhob):= 0.3377372788077926E-01*F000282(rhoa,rhob);
318F000284(rhoa,rhob):= -(F000283(rhoa,rhob));
319F000285(rhoa,rhob):= sqrt(F000228(rhoa,rhob));
320F000286(rhoa,rhob):= F000228(rhoa,rhob)*F000285(rhoa,rhob);
321F000287(rhoa,rhob):= 0.1035700000000000E+02*F000285(rhoa,rhob);
322F000288(rhoa,rhob):= 0.3623100000000000E+01*F000228(rhoa,rhob);
323F000289(rhoa,rhob):= F000287(rhoa,rhob)+F000288(rhoa,rhob);
324F000290(rhoa,rhob):= 0.8802600000000000E+00*F000286(rhoa,rhob);
325F000291(rhoa,rhob):= F000289(rhoa,rhob)+F000290(rhoa,rhob);
326F000292(rhoa,rhob):= 0.4967100000000000E+00*F000228(rhoa,rhob);
327F000293(rhoa,rhob):= F000292(rhoa,rhob)*F000228(rhoa,rhob);
328F000294(rhoa,rhob):= F000291(rhoa,rhob)+F000293(rhoa,rhob);
329F000295(rhoa,rhob):= 0.3377372788077926E-01*F000294(rhoa,rhob);
330F000296(rhoa,rhob):= 0.1000000000000000E+01/F000295(rhoa,rhob);
331F000297(rhoa,rhob):= 0.1000000000000000E+01+F000296(rhoa,rhob);
332F000298(rhoa,rhob):= log(F000297(rhoa,rhob));
333F000299(rhoa,rhob):= F000284(rhoa,rhob)*F000298(rhoa,rhob);
334F000300(rhoa,rhob):= F000230(rhoa,rhob)*F000230(rhoa,rhob);
335F000301(rhoa,rhob):= F000230(rhoa,rhob)*F000300(rhoa,rhob);
336F000302(rhoa,rhob):= F000230(rhoa,rhob)*F000301(rhoa,rhob);
337F000303(rhoa,rhob):= F000242(rhoa,rhob)*F000302(rhoa,rhob);
338F000304(rhoa,rhob):= F000242(rhoa,rhob)-F000303(rhoa,rhob);
339F000305(rhoa,rhob):= 0.5848223622634647E+00*F000304(rhoa,rhob);
340F000306(rhoa,rhob):= 0.1000000000000000E+01-F000303(rhoa,rhob);
341F000307(rhoa,rhob):= F000261(rhoa,rhob)*F000306(rhoa,rhob);
342F000308(rhoa,rhob):= F000280(rhoa,rhob)*F000303(rhoa,rhob);
343F000309(rhoa,rhob):= F000307(rhoa,rhob)+F000308(rhoa,rhob);
344F000310(rhoa,rhob):= F000299(rhoa,rhob)*F000305(rhoa,rhob);
345F000311(rhoa,rhob):= F000309(rhoa,rhob)-F000310(rhoa,rhob);
346F000312(rhoa,rhob):= F000226(rhoa,rhob)*F000311(rhoa,rhob);
347F000313(rhoa,rhob):= F000312(rhoa,rhob)-F000096(rhoa);
348F000314(rhoa,rhob):= F000313(rhoa,rhob)-F000202(rhob);
349F000315(rhoa,rhob,sigmaaa,sigmabb):= F000101(rhoa,sigmaaa)+F000207(rhob,sigmabb);
350F000316(rhoa,rhob,sigmaaa,sigmabb):= 0.3100000000000000E-02*F000315(rhoa,rhob,sigmaaa,sigmabb);
351F000317(rhoa,rhob,sigmaaa,sigmabb):= F000101(rhoa,sigmaaa)+F000207(rhob,sigmabb);
352F000318(rhoa,rhob,sigmaaa,sigmabb):= 0.3100000000000000E-02*F000317(rhoa,rhob,sigmaaa,sigmabb);
353F000319(rhoa,rhob,sigmaaa,sigmabb):= 0.1000000000000000E+01+F000318(rhoa,rhob,sigmaaa,sigmabb);
354F000320(rhoa,rhob,sigmaaa,sigmabb):= F000316(rhoa,rhob,sigmaaa,sigmabb)/F000319(rhoa,rhob,sigmaaa,sigmabb);
355F000321(rhoa,rhob,sigmaaa,sigmabb):= F000320(rhoa,rhob,sigmaaa,sigmabb)*0.0000000000000000E+00;
356F000322(rhoa,rhob,sigmaaa,sigmabb):= 0.0000000000000000E+00+F000321(rhoa,rhob,sigmaaa,sigmabb);
357F000323(rhoa,rhob,sigmaaa,sigmabb):= F000320(rhoa,rhob,sigmaaa,sigmabb)*F000322(rhoa,rhob,sigmaaa,sigmabb);
358F000324(rhoa,rhob,sigmaaa,sigmabb):= -0.1116028770000000E+02+F000323(rhoa,rhob,sigmaaa,sigmabb);
359F000325(rhoa,rhob,sigmaaa,sigmabb):= F000320(rhoa,rhob,sigmaaa,sigmabb)*F000324(rhoa,rhob,sigmaaa,sigmabb);
360F000326(rhoa,rhob,sigmaaa,sigmabb):= 0.5951530800000000E+01+F000325(rhoa,rhob,sigmaaa,sigmabb);
361F000327(rhoa,rhob,sigmaaa,sigmabb):= F000320(rhoa,rhob,sigmaaa,sigmabb)*F000326(rhoa,rhob,sigmaaa,sigmabb);
362F000328(rhoa,rhob,sigmaaa,sigmabb):= 0.1000000000000000E+01+F000327(rhoa,rhob,sigmaaa,sigmabb);
363F000329(rhoa,rhob):= 0.1000000000000000E+01*F000314(rhoa,rhob);
364F000330(rhoa,rhob,sigmaaa,sigmabb):= F000329(rhoa,rhob)*F000328(rhoa,rhob,sigmaaa,sigmabb);
365F000331(rhoa,rhob,sigmaaa,sigmabb):= F000330(rhoa,rhob,sigmaaa,sigmabb)*0.1000000000000000E+01;
366F000332(rhoa,rhob,sigmaaa,sigmabb,taua,taub):= F000225(rhoa,rhob,sigmaaa,sigmabb,taua,taub)+F000331(rhoa,rhob,sigmaaa,sigmabb);
367f(zeta,rhoa,rhob,rho,sigmaaa,sigmaab,sigmabb,sigma,taua,taub,tau):=F000332(rhoa,rhob,sigmaaa,sigmabb,taua,taub);
368F000333:   0.0000000000000000E+00;
369F000334:   0.0000000000000000E+00;
370F000335:   0.0000000000000000E+00;
371F000336:   0.0000000000000000E+00;
372F000337:   0.0000000000000000E+00;
373F000338:   0.0000000000000000E+00;
374F000339:   0.0000000000000000E+00;
375F000340:   0.0000000000000000E+00;
376F000341:   0.0000000000000000E+00;
377F000342:   0.0000000000000000E+00;
378F000343:   0.0000000000000000E+00;
379F000344:   0.0000000000000000E+00;
380F000345:   0.0000000000000000E+00;
381/* comp: rhoa.le.0.5000000000000000E-18
382*/
383F000346(taua):= 0.2000000000000000E+01*taua;
384/* comp: F000346(taua).gt.0.1000000000000000E-17
385*/
386/* comp: rhoa.le.0.1000000000000000E-17
387*/
388/* comp: F000346(taua).le.0.1000000000000000E-17
389*/
390F000347(rhoa):= 0.2387324146378430E+00/rhoa;
391F000348(rhoa):= F000347(rhoa)^(1/3);
392F000349:   0.1000000000000000E+01;
393F000350:  -0.2000000000000000E+01;
394F000351: 0.1000000000000000E+01-F000349;
395F000352: 0.1000000000000000E+01+F000349;
396F000353: F000351^0.2000000000000000E+01;
397F000354: F000352^0.2000000000000000E+01;
398/* comp: F000351.gt.0.1000000000000000E-17
399*/
400/* comp: F000352.gt.0.1000000000000000E-17
401*/
402F000355: F000352^(1/3);
403F000356: F000352*F000355;
404F000357: F000350+F000356;
405F000358: F000357*0.1923661050931536E+01;
406F000359(rhoa):= 0.2137000000000000E+00*F000348(rhoa);
407F000360(rhoa):= 0.1000000000000000E+01+F000359(rhoa);
408F000361(rhoa):= 0.6218140000000000E-01*F000360(rhoa);
409F000362(rhoa):= -(F000361(rhoa));
410F000363(rhoa):= sqrt(F000348(rhoa));
411F000364(rhoa):= F000348(rhoa)*F000363(rhoa);
412F000365(rhoa):= 0.7595700000000000E+01*F000363(rhoa);
413F000366(rhoa):= 0.3587600000000000E+01*F000348(rhoa);
414F000367(rhoa):= F000365(rhoa)+F000366(rhoa);
415F000368(rhoa):= 0.1638200000000000E+01*F000364(rhoa);
416F000369(rhoa):= F000367(rhoa)+F000368(rhoa);
417F000370(rhoa):= 0.4929400000000000E+00*F000348(rhoa);
418F000371(rhoa):= F000370(rhoa)*F000348(rhoa);
419F000372(rhoa):= F000369(rhoa)+F000371(rhoa);
420F000373(rhoa):= 0.6218140000000000E-01*F000372(rhoa);
421F000374(rhoa):= 0.1000000000000000E+01/F000373(rhoa);
422F000375(rhoa):= 0.1000000000000000E+01+F000374(rhoa);
423F000376(rhoa):= log(F000375(rhoa));
424F000377(rhoa):= F000362(rhoa)*F000376(rhoa);
425F000378(rhoa):= 0.2054800000000000E+00*F000348(rhoa);
426F000379(rhoa):= 0.1000000000000000E+01+F000378(rhoa);
427F000380(rhoa):= 0.3109070000000000E-01*F000379(rhoa);
428F000381(rhoa):= -(F000380(rhoa));
429F000382(rhoa):= sqrt(F000348(rhoa));
430F000383(rhoa):= F000348(rhoa)*F000382(rhoa);
431F000384(rhoa):= 0.1411890000000000E+02*F000382(rhoa);
432F000385(rhoa):= 0.6197700000000000E+01*F000348(rhoa);
433F000386(rhoa):= F000384(rhoa)+F000385(rhoa);
434F000387(rhoa):= 0.3366200000000000E+01*F000383(rhoa);
435F000388(rhoa):= F000386(rhoa)+F000387(rhoa);
436F000389(rhoa):= 0.6251700000000000E+00*F000348(rhoa);
437F000390(rhoa):= F000389(rhoa)*F000348(rhoa);
438F000391(rhoa):= F000388(rhoa)+F000390(rhoa);
439F000392(rhoa):= 0.3109070000000000E-01*F000391(rhoa);
440F000393(rhoa):= 0.1000000000000000E+01/F000392(rhoa);
441F000394(rhoa):= 0.1000000000000000E+01+F000393(rhoa);
442F000395(rhoa):= log(F000394(rhoa));
443F000396(rhoa):= F000381(rhoa)*F000395(rhoa);
444F000397(rhoa):= 0.1112500000000000E+00*F000348(rhoa);
445F000398(rhoa):= 0.1000000000000000E+01+F000397(rhoa);
446F000399(rhoa):= 0.3377372788077926E-01*F000398(rhoa);
447F000400(rhoa):= -(F000399(rhoa));
448F000401(rhoa):= sqrt(F000348(rhoa));
449F000402(rhoa):= F000348(rhoa)*F000401(rhoa);
450F000403(rhoa):= 0.1035700000000000E+02*F000401(rhoa);
451F000404(rhoa):= 0.3623100000000000E+01*F000348(rhoa);
452F000405(rhoa):= F000403(rhoa)+F000404(rhoa);
453F000406(rhoa):= 0.8802600000000000E+00*F000402(rhoa);
454F000407(rhoa):= F000405(rhoa)+F000406(rhoa);
455F000408(rhoa):= 0.4967100000000000E+00*F000348(rhoa);
456F000409(rhoa):= F000408(rhoa)*F000348(rhoa);
457F000410(rhoa):= F000407(rhoa)+F000409(rhoa);
458F000411(rhoa):= 0.3377372788077926E-01*F000410(rhoa);
459F000412(rhoa):= 0.1000000000000000E+01/F000411(rhoa);
460F000413(rhoa):= 0.1000000000000000E+01+F000412(rhoa);
461F000414(rhoa):= log(F000413(rhoa));
462F000415(rhoa):= F000400(rhoa)*F000414(rhoa);
463F000416: F000349*F000349;
464F000417: F000349*F000416;
465F000418: F000349*F000417;
466F000419: F000358*F000418;
467F000420: F000358-F000419;
468F000421: 0.5848223622634647E+00*F000420;
469F000422: 0.1000000000000000E+01-F000419;
470F000423(rhoa):= F000377(rhoa)*F000422;
471F000424(rhoa):= F000396(rhoa)*F000419;
472F000425(rhoa):= F000423(rhoa)+F000424(rhoa);
473F000426(rhoa):= F000415(rhoa)*F000421;
474F000427(rhoa):= F000425(rhoa)-F000426(rhoa);
475F000428(rhoa):= rhoa*F000427(rhoa);
476F000429(sigmaaa):= 0.2500000000000000E+00*sigmaaa;
477F000430(rhoa,sigmaaa):= F000429(sigmaaa)/rhoa;
478F000431(rhoa,sigmaaa,taua):= F000346(taua)-F000430(rhoa,sigmaaa);
479F000432(rhoa):= rhoa^(8/3);
480F000433(rhoa,sigmaaa):= sigmaaa/F000432(rhoa);
481F000434(rhoa,sigmaaa):= 0.6000000000000000E-01*F000433(rhoa,sigmaaa);
482F000435(rhoa,sigmaaa):= 0.6000000000000000E-01*F000433(rhoa,sigmaaa);
483F000436(rhoa,sigmaaa):= 0.1000000000000000E+01+F000435(rhoa,sigmaaa);
484F000437(rhoa,sigmaaa):= F000434(rhoa,sigmaaa)/F000436(rhoa,sigmaaa);
485F000438(rhoa,sigmaaa):= F000437(rhoa,sigmaaa)*0.0000000000000000E+00;
486F000439(rhoa,sigmaaa):= 0.0000000000000000E+00+F000438(rhoa,sigmaaa);
487F000440(rhoa,sigmaaa):= F000437(rhoa,sigmaaa)*F000439(rhoa,sigmaaa);
488F000441(rhoa,sigmaaa):= 0.2223379300000000E+01+F000440(rhoa,sigmaaa);
489F000442(rhoa,sigmaaa):= F000437(rhoa,sigmaaa)*F000441(rhoa,sigmaaa);
490F000443(rhoa,sigmaaa):= -0.2596089700000000E+01+F000442(rhoa,sigmaaa);
491F000444(rhoa,sigmaaa):= F000437(rhoa,sigmaaa)*F000443(rhoa,sigmaaa);
492F000445(rhoa,sigmaaa):= 0.1000000000000000E+01+F000444(rhoa,sigmaaa);
493F000446(rhoa,sigmaaa,taua):= F000431(rhoa,sigmaaa,taua)/F000346(taua);
494F000447(rhoa,sigmaaa,taua):= F000446(rhoa,sigmaaa,taua)*F000445(rhoa,sigmaaa);
495F000448(rhoa,sigmaaa,taua):= F000447(rhoa,sigmaaa,taua)*F000428(rhoa);
496F000449(rhoa,sigmaaa,taua):= F000448(rhoa,sigmaaa,taua)*0.1000000000000000E+01;
497F000450(rhoa,sigmaaa,taua):= F000449(rhoa,sigmaaa,taua)*0.1000000000000000E+01;
498F000451(rhoa,sigmaaa,taua):= F000333+F000450(rhoa,sigmaaa,taua);
499/* comp: rhob.le.0.5000000000000000E-18
500*/
501F000452(rhoa):= rhoa+F000341;
502/* comp: rhoa.gt.0.5000000000000000E-18
503*/
504/* comp: F000341.gt.0.5000000000000000E-18
505*/
506F000453(rhoa):= 0.2387324146378430E+00/F000452(rhoa);
507F000454(rhoa):= F000453(rhoa)^(1/3);
508F000455(rhoa):= rhoa-F000341;
509F000456(rhoa):= F000455(rhoa)/F000452(rhoa);
510F000457:  -0.2000000000000000E+01;
511F000458(rhoa):= 0.1000000000000000E+01-F000456(rhoa);
512F000459(rhoa):= 0.1000000000000000E+01+F000456(rhoa);
513F000460(rhoa):= F000458(rhoa)^0.2000000000000000E+01;
514F000461(rhoa):= F000459(rhoa)^0.2000000000000000E+01;
515/* comp: F000458(rhoa).gt.0.1000000000000000E-17
516*/
517/* comp: F000459(rhoa).gt.0.1000000000000000E-17
518*/
519F000462(rhoa):= F000459(rhoa)^(1/3);
520F000463(rhoa):= F000459(rhoa)*F000462(rhoa);
521F000464(rhoa):= F000457+F000463(rhoa);
522F000465(rhoa):= F000464(rhoa)*0.1923661050931536E+01;
523F000466(rhoa):= 0.2137000000000000E+00*F000454(rhoa);
524F000467(rhoa):= 0.1000000000000000E+01+F000466(rhoa);
525F000468(rhoa):= 0.6218140000000000E-01*F000467(rhoa);
526F000469(rhoa):= -(F000468(rhoa));
527F000470(rhoa):= sqrt(F000454(rhoa));
528F000471(rhoa):= F000454(rhoa)*F000470(rhoa);
529F000472(rhoa):= 0.7595700000000000E+01*F000470(rhoa);
530F000473(rhoa):= 0.3587600000000000E+01*F000454(rhoa);
531F000474(rhoa):= F000472(rhoa)+F000473(rhoa);
532F000475(rhoa):= 0.1638200000000000E+01*F000471(rhoa);
533F000476(rhoa):= F000474(rhoa)+F000475(rhoa);
534F000477(rhoa):= 0.4929400000000000E+00*F000454(rhoa);
535F000478(rhoa):= F000477(rhoa)*F000454(rhoa);
536F000479(rhoa):= F000476(rhoa)+F000478(rhoa);
537F000480(rhoa):= 0.6218140000000000E-01*F000479(rhoa);
538F000481(rhoa):= 0.1000000000000000E+01/F000480(rhoa);
539F000482(rhoa):= 0.1000000000000000E+01+F000481(rhoa);
540F000483(rhoa):= log(F000482(rhoa));
541F000484(rhoa):= F000469(rhoa)*F000483(rhoa);
542F000485(rhoa):= 0.2054800000000000E+00*F000454(rhoa);
543F000486(rhoa):= 0.1000000000000000E+01+F000485(rhoa);
544F000487(rhoa):= 0.3109070000000000E-01*F000486(rhoa);
545F000488(rhoa):= -(F000487(rhoa));
546F000489(rhoa):= sqrt(F000454(rhoa));
547F000490(rhoa):= F000454(rhoa)*F000489(rhoa);
548F000491(rhoa):= 0.1411890000000000E+02*F000489(rhoa);
549F000492(rhoa):= 0.6197700000000000E+01*F000454(rhoa);
550F000493(rhoa):= F000491(rhoa)+F000492(rhoa);
551F000494(rhoa):= 0.3366200000000000E+01*F000490(rhoa);
552F000495(rhoa):= F000493(rhoa)+F000494(rhoa);
553F000496(rhoa):= 0.6251700000000000E+00*F000454(rhoa);
554F000497(rhoa):= F000496(rhoa)*F000454(rhoa);
555F000498(rhoa):= F000495(rhoa)+F000497(rhoa);
556F000499(rhoa):= 0.3109070000000000E-01*F000498(rhoa);
557F000500(rhoa):= 0.1000000000000000E+01/F000499(rhoa);
558F000501(rhoa):= 0.1000000000000000E+01+F000500(rhoa);
559F000502(rhoa):= log(F000501(rhoa));
560F000503(rhoa):= F000488(rhoa)*F000502(rhoa);
561F000504(rhoa):= 0.1112500000000000E+00*F000454(rhoa);
562F000505(rhoa):= 0.1000000000000000E+01+F000504(rhoa);
563F000506(rhoa):= 0.3377372788077926E-01*F000505(rhoa);
564F000507(rhoa):= -(F000506(rhoa));
565F000508(rhoa):= sqrt(F000454(rhoa));
566F000509(rhoa):= F000454(rhoa)*F000508(rhoa);
567F000510(rhoa):= 0.1035700000000000E+02*F000508(rhoa);
568F000511(rhoa):= 0.3623100000000000E+01*F000454(rhoa);
569F000512(rhoa):= F000510(rhoa)+F000511(rhoa);
570F000513(rhoa):= 0.8802600000000000E+00*F000509(rhoa);
571F000514(rhoa):= F000512(rhoa)+F000513(rhoa);
572F000515(rhoa):= 0.4967100000000000E+00*F000454(rhoa);
573F000516(rhoa):= F000515(rhoa)*F000454(rhoa);
574F000517(rhoa):= F000514(rhoa)+F000516(rhoa);
575F000518(rhoa):= 0.3377372788077926E-01*F000517(rhoa);
576F000519(rhoa):= 0.1000000000000000E+01/F000518(rhoa);
577F000520(rhoa):= 0.1000000000000000E+01+F000519(rhoa);
578F000521(rhoa):= log(F000520(rhoa));
579F000522(rhoa):= F000507(rhoa)*F000521(rhoa);
580F000523(rhoa):= F000456(rhoa)*F000456(rhoa);
581F000524(rhoa):= F000456(rhoa)*F000523(rhoa);
582F000525(rhoa):= F000456(rhoa)*F000524(rhoa);
583F000526(rhoa):= F000465(rhoa)*F000525(rhoa);
584F000527(rhoa):= F000465(rhoa)-F000526(rhoa);
585F000528(rhoa):= 0.5848223622634647E+00*F000527(rhoa);
586F000529(rhoa):= 0.1000000000000000E+01-F000526(rhoa);
587F000530(rhoa):= F000484(rhoa)*F000529(rhoa);
588F000531(rhoa):= F000503(rhoa)*F000526(rhoa);
589F000532(rhoa):= F000530(rhoa)+F000531(rhoa);
590F000533(rhoa):= F000522(rhoa)*F000528(rhoa);
591F000534(rhoa):= F000532(rhoa)-F000533(rhoa);
592F000535(rhoa):= F000452(rhoa)*F000534(rhoa);
593F000536(rhoa):= F000535(rhoa)-F000428(rhoa);
594F000537(rhoa):= F000536(rhoa)-F000338;
595F000538(rhoa,sigmaaa):= F000433(rhoa,sigmaaa)+F000339;
596F000539(rhoa,sigmaaa):= 0.3100000000000000E-02*F000538(rhoa,sigmaaa);
597F000540(rhoa,sigmaaa):= F000433(rhoa,sigmaaa)+F000339;
598F000541(rhoa,sigmaaa):= 0.3100000000000000E-02*F000540(rhoa,sigmaaa);
599F000542(rhoa,sigmaaa):= 0.1000000000000000E+01+F000541(rhoa,sigmaaa);
600F000543(rhoa,sigmaaa):= F000539(rhoa,sigmaaa)/F000542(rhoa,sigmaaa);
601F000544(rhoa,sigmaaa):= F000543(rhoa,sigmaaa)*0.0000000000000000E+00;
602F000545(rhoa,sigmaaa):= 0.0000000000000000E+00+F000544(rhoa,sigmaaa);
603F000546(rhoa,sigmaaa):= F000543(rhoa,sigmaaa)*F000545(rhoa,sigmaaa);
604F000547(rhoa,sigmaaa):= -0.1116028770000000E+02+F000546(rhoa,sigmaaa);
605F000548(rhoa,sigmaaa):= F000543(rhoa,sigmaaa)*F000547(rhoa,sigmaaa);
606F000549(rhoa,sigmaaa):= 0.5951530800000000E+01+F000548(rhoa,sigmaaa);
607F000550(rhoa,sigmaaa):= F000543(rhoa,sigmaaa)*F000549(rhoa,sigmaaa);
608F000551(rhoa,sigmaaa):= 0.1000000000000000E+01+F000550(rhoa,sigmaaa);
609F000552(rhoa):= 0.1000000000000000E+01*F000537(rhoa);
610F000553(rhoa,sigmaaa):= F000552(rhoa)*F000551(rhoa,sigmaaa);
611F000554(rhoa,sigmaaa):= F000553(rhoa,sigmaaa)*0.1000000000000000E+01;
612F000555(rhoa,sigmaaa,taua):= F000451(rhoa,sigmaaa,taua)+F000554(rhoa,sigmaaa);
613G(rhoa,sigmaaa,taua,taub):=F000555(rhoa,sigmaaa,taua);
614