文档介绍：天然气管道流量压强编程计算方法
#include <>
#include <>
void result(double q[11],double l[11],double Q1[2],double Q10[1;
}
if(x2>0&&y2>0)
{
Q1[1]=x2;
Q10[1]=y2;
}
}
}
int main()
{
double q[11]={0},l[11]={0},Q1[2]={0},Q10[2]={0},c[11]={0},f[11]={0},D,flow[10]={0},sum=0,P[11];
int i,j,t=0,m;
for(i=1;i<=10;i++)
{
scanf("%lf %lf",&q[i],&l[i]);
}
result(q,l,Q1,Q10);
t=t+1;
if(Q1[1]>0&&Q10[1]>0)
{
printf("%f %f\n",Q1[1],Q10[1]);
}
else
{
for(t=2;t<=10;t++)
{
for(i=1;i<=10;i++)
{
c[i]=q[i];
f[i]=l[i];
}
q[10]=c[1];
l[10]=f[1];
for(i=1;i<=9;i++)
{
q[i]=c[i+1];
l[i]=f[i+1];
}
result(q,l,Q1,Q10);
if(Q1[1]>0&&Q10[1]>0)
{
printf("分气点位置是%d号井\n",t);
printf("分气点顺时针方向流量是%fX10^5立方米/天\n",Q1[1],Q10[1]);
m=7-t+1;
break;
}
}
}
flow[1]=Q1[1];
flow[10]=Q10[1];
for(i=2;i<=9;i++)
{
if(i<m)
{
flow[i]=Q1[1];
for(j=2;j<=i;j++)
flow[i]=flow[i]+q[j];
}
else
{
flow[i]=Q10[1];
for(j=i+1;j<=10;j++)
flow[i]=flow[i]+q[j];
}
}