#include"math.h"
#include"stdio.h"
double P,d;
double fd(P,d)
{ double a;
    a=P*d/(2.0*1.0*0.5*380*1e6*1.0);                         //計算平均壁厚//
 return (a);
 }
void main() 
 { double Pq,Pz,Ppj,Px;
  double T,Qmax,Q,Q1,C,C0,z,R,de,A,B,PI,l0,l1,l2;
  double Aq1,Aq2,Azh,Pde,p0;
  double sum[4],d[4]={0.5862,0.5846,0.6346,0.6314};
  double steel[4],pump,resume,oper;
  int i,a;
  Qmax=30.0*1e8/(3600.0*24*350);
  Q=Qmax*293.15/273.15;                                   //換算到標況//
  Pz=0.9*1e6;Pq=6.0*1e6;l0=89000;l1=1615200;l2=1924000;
  C0=0.03848; z=0.9126133; R=0.0117103;de=0.68659;Px=4225714;
  A=10.0;PI=3.1415926;T=283.15;p0=7.8;Ppj=1750000  ;
  B=(A*Pz*Pz-Pq*Pq)/(Q*Q);
   for(i=0;i<4;i++)
   {C=R*z*de*T/(C0*C0*pow(d[i],5));
   Aq1=fd(Pq,d[i]);
   Aq2=fd(Px,d[i]);                                       //計算用鋼量//
   Azh=fd(Ppj,d[i]);
   printf("Aq1=%e,Aq2=%e,Azh=%e\n",Aq1,Aq2,Azh);
   steel[i]=(Aq1*l0+Aq2*l1+Azh*l2)*PI*d[i]*p0*10000.0;    //鋼材費用//
   pump=805*1e6;                                           //泵 費 用//
   Pde=(Pq-Pz);
   Q1=Qmax/(Ppj/101325);
   resume=Pde*350.0*24*Q1*0.7/1000;                        //運行費用//
   oper=5.0*1e5;
   sum[i]=steel[i]+pump+resume+oper;
   printf("steel[%d]=%e,pump=%e,resume=%e,oper=%e\n",i,steel[i],pump,resume,oper);
}
  a=0;
  for(i=0;i<4;i++)
  {if(sum[i]<sum[a]) a=i;
  printf("d[%d]=%e,sum[%d]=%e\n",i,d[i],sum[i]);
  }
}