# ifdef _DEBUG
# define new DEBUG_NEW
# endif

# define ideal coefficient a1/8.37304456 * * parameter macro definition//* */
# define a2-3.70454304
# define a3 2.5000
# define a4 1.99427042
# define the a5 0.62105248
# define a6 0.41195293
# define a7 1.04028922
# define a8 0.08327678
# define b4 3.15163
# define b5 6.11190
# define b6 6.77708
# define b7 11.32384
# define b8 27.08792


/* * : define global array, returns an array of all numerical calculation, the main function of the import TXT text
Statement: w=nondimensional density, y=dimensionless temperature. This function w, y are from 1 to 50, variation is 1, so the partial derivative calculation result is the same (w, y need to edit difference function) * */


Float lx [50];/* array of globally defined */
Float cy [50];
Float HMHZ [50].
Float ca [50];
Float cb [50];
Float cc [50];
Float CD [50];
Float ce [50];
Float cf [50];
Float CG [50];
Float P [50];
Float s [50];
The float u [50];
Float cr [50];
Float h [50];


Void HMHZLX ();/* function header defines */
Void hmhzcy ();
Void hm ();
Float dchafen ();
Float schafen ();
Void piandao ();
Void shuchu ();

Void main ()/* */
main function{
HMHZLX ();/* executable program */
Hmhzcy ();
Hm ();



Piandao ();
Shuchu ();

Fc fa, the FILE * * fb, * and * fd, * fe * (ff, fg, * fh, * fi * fj, * fk1 - fk, fl * * FM, * fn, * fo;/* the output array into TXT text */
Fa=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ lx. TXT ", "wb");
Fb=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ cy TXT ", "wb");
Fc=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ HMHZ TXT ", "wb");
Fd=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ a.t xt ", "wb");
Fe=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandap \ \ b.t xt ", "wb");
Ff=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ c.t xt ", "wb");
Fg=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ which xt ", "wb");
Fh=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ e. xt ", "wb");
Fi=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ f.t xt ", "wb");
Fj:=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ piandao \ \ g.t xt ", "wb");
Fk1 - fk=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ P.t xt ", "wb");
Fl=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ s.t xt ", "wb");
FM=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ u.t xt ", "wb");
Fn=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ cr TXT ", "wb");
Fo=fopen (" C: \ \ Users \ \ acer \ \ Docktop \ \ jieguo \ \ h.t xt ", "wb");

int i;
For (I=0; I & lt; 50; I++)
{
Printf (" % f % f % % % % % f f f f f % % % % % f f f f f % f % f, % f \ n ", lx [I], cy [I], HMHZ [I], ca [I], cb [I], cc [I], [I], CD ce [I], cf [I], CG [I], [I] P, s [I], u [I], cr [I], h [I]);
Fprintf (fa, "% f \ n", lx [I]);
Fprintf (fb, % f \ "n", cy [I]);
Fprintf (fc, % f \ "n", HMHZ [I]);
Fprintf (fd, % f \ "n", ca [I]);
Fprintf (fe, "% f \ n", cb [I]);
Fprintf (ff, % f \ "n", cc [I]);
Fprintf (fg % f "\ n", CD [I]);
Fprintf (fh, % f \ "n", ce [I]);
Fprintf (fi, % f \ "n", cf [I]);
Fprintf (fj, % f \ "n", CG [I]);
Fprintf (fk1 - fk, % f \ "n", P [I]);
Fprintf (fl, % f \ "n", s [I]);
Fprintf (FM, % f \ "n", u [I]);
Fprintf (fn, % f \ "n", cr [I]);
Fprintf (fo, % f \ "n", h [I]);
}

The fclose (fa);
The fclose (fb);
The fclose (fc);
The fclose (fd);
The fclose (fe);
The fclose (ff);
The fclose (fg);
The fclose (fh);
The fclose (fi);
The fclose (fj);
The fclose (fk1 - fk);
The fclose (fl);
The fclose (FM);
The fclose (fn);
The fclose (fo);
}


Void HMHZLX ()
{
Float w=1, y=1;
int i;
For (I=0; I & lt; 50; I++)
{
Lx [I]=log (w) + a1 + (2 * y) + (a3 * log (y)) + (a4 * log (1 - exp (b4) - y *)) + (a5 * log (1 - exp (- y * b5))) + (a6 * log (1 - exp (b6) - y *)) + (a7 * log (1 - exp (- y * b7))) + (a8 * log (1 - exp (- y * b8)));

W++;
Y++;
}
}

Void hmlxcy ()
{
float n[42]={ 0.38856823203161e0,0.29385475942740e1,-0.55867188534934e0,-0.76753199592477e0,0.31729005580416e0,0.54803315897767e0,0.12279411220335e0,0.21658961543220e1,0.15841735109724e1,-0.23132705405503e0,0.5811691643146e-1,-0.55369137205382e0,0.48946615909422e0,-0.24275739843501e-1,0.62494790501e-1,-0.12175860225246e0,-0.37055685270086e0,-0.16775879700426e-1,-0.11960736637987e0,-0.45619362508778e-1,0.35612789270346e-1,-0.74427727132052e-2,-0.17395704902432e-2,-0.21810121289527e-1,0.24332166559236e-1,-0.37440133423463e-1,0.14338715756878e0,-0.13491969083286,-0.23151225053480e-1,0.12363125492901e-1,0.21058321970940e-2,-0.3358519026368e-3,0.55993651771592e-2,-0.30335118055646e-3,-0.21365488688320e3,0.26641569149272e5,-0.24027212204557e5,-0.28341603423999e3,0.2124728440017e3,-0.66642276540751e0,0.72608632349897e0,0.55068668612842e-1 };
Float d [39]=,1,1,1,2,2,3,1,2,4,5,5,5,6,6,6,1,1,4,4,4,7,8,2,3,3,5,5,6,7,8,10,4,8,2,2,2,3,3 {1};
Float t [39]={,0.75,2,0.75 0,0.75, 1, 2, 1.5, 1.5, 2.5, 0,1.5, 2,0,1,2,3,6,3,6,8,6,0,7,12,16,22,24,16,24,8,2,28,14,1,0,1,3,3};
Float c [27]=,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,3,3,3,4,4,4,4,4,4,5,6 {1};
Float alpha [5]=,25,25,25,20 {25};
Float beta [8]={325300300275275,0.3, 0.3, 0.3};
nullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnullnull