分光计的不确定度的程序

#include<iostream>
#include<cmath>
#define pi 3.1415926
using namespace std;
int main()
{
	double sum1 = 0, sum2 = 0, asum1 = 0.0,asum2=0;
	double s[5], average1, average2, u1, u2;
	double U1, U2, U, B, d[5];
	int i;
	cout << "please input the dates:" << endl;
	for (i = 0; i < 10; i++)
	{
		cin >> s+;
		d+ = ((s+ - int(s+)) * 100 / 60 + int(s+)) * (pi / 180);
	}
	for (i = 0; i < 5; i++)
	{
		sum1 = d+ + sum1;
	}
	average1 = sum1 / 5;
	for (i = 5; i < 10; i++)
	{
		sum2 = d+ + sum2;
	}
	average2 = sum2 / 5;
	for (i = 0; i < 5; i++)
	{
		asum1 = pow((d+ - average1), 2) + asum1;
	}
	for (i = 5; i < 10; i++)
	{
		asum2 = pow((d+ - average2), 2) + asum2;
	}
	B = (1.0 / 60) * (pi / 180) / sqrt(3);
	u1 = sqrt(asum1 / 20);
	u2 = sqrt(asum2 / 20);
	U1 = sqrt(pow(u1, 2) + pow(0.1 / sqrt(3), 2));
	U2 = sqrt(pow(u2, 2) + pow(0.1 / sqrt(3), 2));
	U = sqrt(pow(U1, 2) + pow(U2, 2));
	cout << "a类不确定度为：A1=" << u1 << "  A2=" << u2 << endl;
	cout << "b类不确定度为：B=" << B << endl;
	cout << "c类不确定度为：C1=" << U1 <<"  C2"<<U2<< endl;
	cout << "U=" << U << endl;
	cout << "平均值"<<(average1 + average2) / 2 << endl;
	system("pause");
	return 0;
}