物理实验计算不确定度的小程序c语言

#include <stdio.h>
#include <math.h>
int main() {
	double a[100]= {0}; //input your numbers
	double ans=0;//answer
	int n=0,i=0;//this means how many numbers you have put
	double sum=0,average=0;
	double ins=0,k=0,ub=0,ua=0,u=0;//ins是误差允许极限，k是B类不确定度的置信因子 
	printf("请输入误差允许极限ins\n");
	scanf("%lf",&ins); 
	printf("请输入置信因子k\n");
	scanf("%lf",&k); 
	ub=ins/k;
	printf("请输入数据，输入0记为结束\n");
	do {
		scanf("%lf",&a[n]);
		n++;
	} while(a[n-1]!=0);
	for(i=0; i<=n; i++) {
		sum+=a[i];
	}
	average=sum/(n-1);
	sum=0;
	i=0;
	for(i=0;i<=n-2;i++){
		sum+=((a[i]-average)*(a[i]-average));
	}
	ans=sqrt(sum/(n-2));
	ua=ans/(sqrt(n-1));
	
	u=sqrt(ua*ua+ub*ub);
	
	printf("平均值的A类不确定度是%lf\n",ua);
	printf("平均值的B类不确定度是%lf\n",ub);
	printf("平均值的合成不确定度是%lf\n",u);
	return 0;
}