//2014111161
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
double Re1,Re2,V1,V2,a1,a2,b1,b2;
cout<<"请假设两管道沿程阻力系数:";
cin>>a1>>a2;
do
{b1=a1;b2=a2;V1=pow((117.68/(0.69+154.3*a1+52.67*a2)),0.5);
V2=(V1*0.36)/0.81;Re1=(V1*0.6)/0.000001;Re2=(V2*0.9)/0.000001;
a1=0.0055*(1+pow((50+1000000/Re1),1.0/3));
a2=0.0055*(1+pow((6.67+1000000/Re2),1.0/3));
}
while((b1-a1)/b1<0.001&&(b2-a2)/b2<0.001);
cout<<"两管道流速分别为:"<<V1<<endl<<V2<<endl<<"流量为:"<<0.09*3.14*V1<<endl;
return 0;
}
Ly8yMDE0MTExMTYxCiNpbmNsdWRlPGlvc3RyZWFtPgojaW5jbHVkZTxjbWF0aD4KdXNpbmcgbmFtZXNwYWNlIHN0ZDsKaW50IG1haW4oKQp7CmRvdWJsZSBSZTEsUmUyLFYxLFYyLGExLGEyLGIxLGIyOwpjb3V0PDwi6K+35YGH6K6+5Lik566h6YGT5rK/56iL6Zi75Yqb57O75pWwOiI7CmNpbj4+YTE+PmEyOwpkbwp7YjE9YTE7YjI9YTI7VjE9cG93KCgxMTcuNjgvKDAuNjkrMTU0LjMqYTErNTIuNjcqYTIpKSwwLjUpOwpWMj0oVjEqMC4zNikvMC44MTtSZTE9KFYxKjAuNikvMC4wMDAwMDE7UmUyPShWMiowLjkpLzAuMDAwMDAxOwphMT0wLjAwNTUqKDErcG93KCg1MCsxMDAwMDAwL1JlMSksMS4wLzMpKTsKYTI9MC4wMDU1KigxK3BvdygoNi42NysxMDAwMDAwL1JlMiksMS4wLzMpKTsKfQp3aGlsZSgoYjEtYTEpL2IxPDAuMDAxJiYoYjItYTIpL2IyPDAuMDAxKTsKY291dDw8IuS4pOeuoemBk+a1gemAn+WIhuWIq+S4ujoiPDxWMTw8ZW5kbDw8VjI8PGVuZGw8PCLmtYHph4/kuLo6Ijw8MC4wOSozLjE0KlYxPDxlbmRsOwpyZXR1cm4gMDsKfQ==