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