# include<iostream>
# include<cmath>
using namespace std;
int main()
{
double d1, d2,d3,d4,d5,d6,d7,d8;
double d;
double t1,t2,t3,t4,t5;
double t;
double value;
double v;
double biaozhun;
double wucha;
double Re;
double value1;
cout<<"输入小球8组直径:"<<endl;
cin>>d1>>d2>>d3>>d4>>d5>>d6>>d7>>d8;
d=(d1+d2+d3+d4+d5+d6+d7+d8)/8000;
cout<<"输出小球直径平均值:"<<d*1000<<"mm"<<endl;
double p=7.8e3,p0=0.95e3,g=9.8,L=0.2,D=0.02;
cout<<"输入第一组温度测量时间"<<endl;
while(1)
{
cin>>t1>>t2>>t3>>t4>>t5;
t=(t1+t2+t3+t4+t5)/5;
cout<<"输出所需时间平均值为:"<<t<<'s'<<endl;
v=0.2/t;
cout<<"输出速度为:"<<v<<"m/s"<<endl;
value=((p-p0)*g*d*d*t)/(18*L*(1+2.4*d/D));
cout<<"该温度下粘滞系数为"<<value<<endl;
cout<<"此温度下粘度标准值:";
cin>>biaozhun;
wucha=fabs(value-biaozhun)/biaozhun;
cout<<"相对误差为:"<<wucha*100<<'%'<<endl;
cout<<endl<<endl;
cout<<"输入下一组温度测量时间:"<<endl;
Re=((v*d*p0)/value);
cout<<"该温度下的雷诺系数:"<<Re<<endl;
value1=(value*(1/(1+3*Re/16)));
cout<<"修正:"<<value1<<endl;
}
return 0;
}