重庆邮电大学 大学物理实验 杨氏模量计算 c++实现

网上居然找不到有人做。自己写了一个。我不保证写的是对的，因为我实验报告还没交。

以下是源代码，写得很烂：

#include <iostream>
#include <cmath>
using namespace std;
#define Pi 3.14
double a[10],b[10],c[10],d[5];
double E;
double L,H,D;
double d_si0,d_si=0,d_ping;

double uncertainty(int len,double A[]){
	double sum=0;
	for(int i=0;i<len;i++){
		sum+=A[i];
	}
	sum/=len;
	double Sn=0;
	for(int i=0;i<len;i++){
		Sn+=((A[i]-sum)*(A[i]-sum));
	}
	return sqrt(Sn/(len*(len-1)));
}

int main(){
	cout<<"请依次输入L、H、D（单位mm），以空格间隔："<<endl;
	cin>>L>>H>>D;
	cout<<"请输入d0（单位mm）：";
	cin>>d_si0;
	double d_s[6];
	cout<<"请输入测量的六组金属丝直径视值d（单位mm），以空格间隔："<<endl;
	for(int i=0;i<6;i++){
		float temp;
		cin>>temp;
		d_s[i]=temp;
		d_si+=(temp-d_si0);
	}
	d_ping=d_si/6;
	cout<<"接下来，请成对输入相同力下的x值。单位mm，以空格间隔。"<<endl;
	for(int i=0;i<10;i++){
		cout<<"请输入第"<<i+1<<"组数据：";
		if (scanf("%lf%lf", &a[i], &b[i]) != 2) {
			cout << "Error: invalid input" << endl;
			return 1;
		}
		c[i]=(a[i]+b[i])/2;
	}
	for(int i=0;i<5;i++){
		d[i]=c[i+5]-c[i];
	}
	double E0[5];
	for(int i=0;i<5;i++){
		E=((8*5*9.8*(L*0.001)*(H*0.001))/(Pi*(d_ping*0.001)*(d_ping*0.001)*(D*0.001)))/(d[i]*0.001);
		E0[i]=E;
		cout<<"E"<<i+1<<"="<<E<<endl;
	}
	cout<<"d的不确定度："<<uncertainty(6,d_s)<<endl;
	cout<<"Δx的不确定度为："<<uncertainty(5,d)<<endl;
	cout<<"杨氏模量的不确定度为："<<uncertainty(5,E0)<<endl;
	getchar();
	getchar();
	return 0;
}

如图：

怎么用就不需要我说了吧。