#include <stdio.h>
#include <math.h>
#include <conio.h>



float a[20][20]={{4.7,0.3,2.0,3.74e-3,6.1,2.45,0.722,2.498},{6.5,0.2,1.4,5.4e-3,10.714,2.02,0.722,2.93},
			   {6.5,0.3,2.0,5.27e-3,7.9,2.67,0.785,3.24},{9.5,0.6,2.0,8.37e-3,11.1,3.02,0.785,3.03},
			  {9.5,0.6,4.2,7.63e-3,5.952,5.13,0.785,4.88},{4.7,0.3,2.0,3.74e-3,6.1,2.45,0.722,2.5},
			   {6.5,0.2,1.4,5.4e-3,10.714,2.02,0.833,2.93},{9.5,0.2,1.4,7.97e-3,15.0,2.13,0.885,4.1},
			  {9.5,0.2,1.7,8.21e-3,12.706,2.58,0.861,3.47},{4.7,0.3,2.0,3.47e-3,6.1,2.45,0.65,2.5},
			  {6.5,0.2,1.4,5.4e-3,10.714,2.02,0.833,2.93},{6.5,0.2,1.7,5.56e-3,9.176,2.42,0.8,2.51},
			  {6.5,0.3,0.2,5.27e-3,7.9,2.67,0.766,3.24}};
float t11,t12,t21,t22,q1,q2,eff,c,deltm,tm2,tm1,mu1,mu2,lamda1,lamda2,cp1,cp2,pr1,pr2,
	  n1,n2,l1,l2,l3,bs,lamdaf1,lamdaf2,ac1,ac2,gm1,gm2,re1,re2,nu1,nu2,alpha1,alpha2,mh1,mh2,
	  yitaf1,yitaf2,af1,af2,a1,a2,ap1,ap2,aeff1,aeff2,ap,rw,ka,k1,k2,w1,w2,cxing,ntu,yitai,
	  yita,fai,t12n,t22n,deltayita,deltap,deltas;
void main()
{
//  clrscr();
  
  float mupara(float);
  float lamdapara(float);
  float cppara(float);
  float prpara(float);
  float max(float,float);
  float min(float,float);
  float calcmass(int i,float l1,float l2,float l3);

int i;

t11=250.0;t21=90.0;q1=0.115;q2=1;eff=0.93;deltap=0.5;deltas=2;
  eff=eff+0.01;

  t12=t11-eff*(t11-t21);

  t22=t21+q1/q2*(t11-t12);
  c=q1/q2;
  deltm=((t11-t22)-(t12-t21))/(log((t11-t22)/(t12-t21)));
  if (c<0.5)
     {tm2=(t21+t22)/2;tm1=tm2+deltm;}
  else  {tm1=(t11+t12)/2;tm2=(t21+t22)/2;}

  mu1=mupara(tm1);
 mu2=mupara(tm2);
  lamda1=lamdapara(tm1);
  lamda2=lamdapara(tm2);
  cp1=cppara(tm1);
  cp2=cppara(tm2);
  pr1=prpara(tm1);
  pr2=prpara(tm2);
printf("%f,%f,%f,%f,%f,%f,%f,%f\n",mu1,mu2,lamda1,lamda2,cp1,cp2,pr1,pr2);


for (i=0;i<=8;i++)
{

n1=7;n2=8;l1=220;l2=60;bs=6;
l3=n1*(a[i][0]+2*deltap)+n2*a[i][0]+2*deltas;
lamdaf1=lamdaf2=175;
ac1=n1*(l2-1.5*bs)*1e-3*(a[i][0]-a[i][0]*a[i][1]/a[i][2])*1e-3;
ac2=n2*(l1-2*bs)*1e-3*(a[i][0]-a[i][0]*a[i][1]/a[i][2])*1e-3;

gm1=q1/ac1;gm2=q2/ac2;
re1=gm1*a[i][5]*1e-3/mu1;re2=gm2*a[i][5]*1e-3/mu2;
nu1=0.06*pow(re1,(2.0/3.0));nu2=0.06*pow(re2,(2.0/3.0));
alpha1=nu1*lamda1/a[i][5]/1e-3;alpha2=nu2*lamda2/a[i][5]/1e-3;
mh1=sqrt(2*alpha1/lamdaf1/a[i][1]/1e-3)*a[i][0]/2*1e-3;
mh2=sqrt(2*alpha2/lamdaf2/a[i][1]/1e-3)*a[i][0]/2*1e-3;
yitaf1=tanh(mh1)/mh1;yitaf2=tanh(mh2)/mh2;
af1=n1*(l2-1.5*bs)*1e-3/a[i][2]*a[i][0]*l1*2*1e-3;
af2=n2*(l1-2*bs)*1e-3/a[i][2]*a[i][0]*l2*2*1e-3;
a1=af1/a[i][6];a2=af2/a[i][6];
a1=af1/a[i][6];a2=af2/a[i][6];
ap1=a1-af1;
ap2=a2-af2;
aeff1=ap1+yitaf1*af1;aeff2=ap2+yitaf2*af2;
ap=2*n1*(l1-2*bs)*(l2-1.5*bs)*1e-6;
rw=0.5e-3/175/ap;
ka=1/(1/aeff1/alpha1+rw+1/aeff2/alpha2);
k1=ka/a1;
w1=q1*cp1*1e3;w2=q2*cp2*1e3;
cxing=min(w1,w2)/max(w1,w2);
ntu=ka/min(w1,w2);
yitai=1-exp(pow(ntu,0.22)/cxing*(exp(-cxing*pow(ntu,0.78))-1));
yita=(pow((1-cxing*yitai)/(1-yitai),2)-1)/(pow((1-cxing*yitai)/(1-yitai),2)-cxing);
fai=yita*min(w1,w2)*(t11-t21);
t12n=t11-fai/w1;t22n=t21+fai/w2;

deltayita=((yita-eff)/eff)<0?((eff-yita)/eff):((yita-eff)/eff);

if (deltayita<0.01&&calcmass(i,l1,l2,l3)<=4) a[i][8]=1;

printf("%d,%f,%f,%f,%f\n",i,yita,deltayita,calcmass(i,l1,l2,l3),a[i][8]);

}

scanf("%d",&i);


}
float mupara(float t)
{float a[100]={0,14.6,15.2,15.7,16.2,16.7,17.2,17.6,18.1,18.6,19.1,19.6,20.1,20.6,21.1,21.5,21.9,22.8,23.7,
		    24.5,25.3,26.0,27.4,29.7,31.4,33.0,36.2,39.1,41.8,44.3,46.7,49.0,51.2,53.5};
 float z;
t=0.1*t+6;
	z=a[floor(t)]+(t-floor(t))*(a[floor(t)+1]-a[floor(t)]);z=z*1e-6;
return(z);
}


float lamdapara(float t)
{float a[100]={0,2.04,2.12,2.20,2.28,2.36,2.44,2.51,2.59,2.67,2.76,2.83,2.90,2.96,3.05,3.13,3.21,3.34,3.49,3.64,3.78,3.93,4.27,4.60,4.91,5.21,5.74,6.22,6.71,7.18,7.63,8.07,8.50,9.15};
 float z;
t=0.1*t+6;
	z=a[floor(t)]+(t-floor(t))*(a[floor(t)+1]-a[floor(t)]);z=z*1e-2;
return(z);
}

float cppara(float t)
{float a[100]={0,1.103,1.013,1.013,1.009,1.009,1.005,1.005,1.005,1.005,1.005,1.005,1.005,
1.009,1.009,1.009,1.009,1.009,1.013,1.017,1.022,1.026,1.038,1.047,1.059,1.068,1.093,1.114,1.135,1.156,1.172,1.185,1.197,1.210};
 float z;
t=0.1*t+6;
	z=a[floor(t)]+(t-floor(t))*(a[floor(t)+1]-a[floor(t)]);
return(z);
}

float prpara(float t)
{float a[100]={0,0.728,0.728,0.723,0.716,0.712,0.707,0.705,0.703,0.701,0.699,0.698,0.696,0.694,0.692,
0.690,0.688,0.686,0.684,0.682,0.681,0.680,0.677,0.674,0.678,0.687,0.699,0.706,0.713,0.717,
0.719,0.722,0.724};
 float z;
t=0.1*t+6;
	z=a[floor(t)]+(t-floor(t))*(a[floor(t)+1]-a[floor(t)]);
return(z);
}
float max(float x,float y)
{float z;z=x>y?x:y;return(z);}
float min(float x,float y)
{float z;z=x<y?x:y;return(z);}

float calcmass(int i,float l1,float l2,float l3)
{float rou1,rou2,h,mfl,mcov,mb,m,mflt,l11,l21,delta,n11,mf1,mf2,mf,mps,ms1,ms2,ms,mcore;
rou1=2.7307e3;rou2=2.639e3;n11=16;
 l11=248;l21=134;delta=3;h=35;
 mf1=n1*(2*l2-3*bs)/a[i][2]*sqrt(a[i][0]*a[i][0]+a[i][2]*a[i][2])*l1*a[i][1]*rou1*1e-9;
 mf2=n2*(l1-2*bs)/a[i][2]*sqrt(a[i][0]*a[i][0]+a[i][2]*a[i][2])*2*l2*a[i][1]*rou1*1e-9;
 mf=mf1+mf2;
 mps=(2*n1*deltap+2*deltas)*l1*2*l2*rou1*1e-9;
 ms1=3*n1*l1*bs*a[i][0]*rou1*1e-9;
 ms2=2*n2*2*l2*bs*a[i][0]*rou1*1e-9;
 ms=ms1+ms2;
 mcore=mf+mps+ms;
 mfl=4*(l11+l21-2*n11)*n11*delta*rou2*1e-9;
mflt=mfl+4*(l1+l21-2*n11)*10*2*rou1*1e-9;
mcov=(2*l3*2*l2+(4*2*l2+5*l3)*h)*delta*rou2*1e-9;
mb=2*((62*62-42*42)*10+(46*46-42*42)*20)*3.141/4.0*rou2*1e-9;
m=mcore+mflt+mcov+mb;
return(m);
}