#include"math.h"
#include"stdio.h"
void main()
{
	double c,c0,dm,Rm,zm,de1,T,Qm,Lm;
	double P1max,P2max,P1min,P2min,Ppjmax,Ppjmin;
	double Vs,V,Vmin,Vmax,P0,T0;
	double Pi=3.1415926;Vmin=0;Vmax=0.0;
	T0=273;P0=1.01325*1e5;
	Rm=1.328802e-002;   zm=9.192528e-001;  de1=0.5624;  T=300.15;
	c0=0.03848;  dm=0.3436;    Qm=6.39*1e8/(350*24*3600);
	P1max=6.3*1e6;  P2min=2.0*1e6;  Lm=208000;
	c=Rm*zm*de1*T/(c0*c0*pow(dm,5));
	P2max=sqrt(P1max*P1max-c*Lm*Qm*Qm);
	P1min=sqrt(P2min*P2min+c*Lm*Qm*Qm);
	Ppjmin=2/3.0*(P1min+P2min*P2min/(P1min+P2min));
    Ppjmax=2/3.0*(P1max+P2max*P2max/(P1max+P2max));
	V=Pi*dm*dm/4*Lm;
	Vmin=Ppjmin*V*zm*T0/(P0*zm*T);
    Vmax=Ppjmax*V*zm*T0/(P0*zm*T);
	Vs=Vmax-Vmin;
	printf("調峰能力為Vs=%e立方米\n",Vs);
}