C语言一元线性回归方程

#include <stdio.h> 
#include <stdlib.h> 
void main() {
	float x[] = {0,5182,10391,15628,20791,26001,31228,36419,41683,46846,52083,57292,62263,67674,72966,83375};
	float y[] = {0,5000,10000,15000,20000,25000,30000,35000,40000,45000,50000,55000,60000,65000,70000,80000};

	int n; 
	n = sizeof(x) / sizeof(x[0]);
	float a, b, b1, mxy, sum_x, sum_y, lxy, xiSubSqr;
	a = b = mxy = sum_x = sum_y = lxy = xiSubSqr = 0.0;
	for (int i = 0; i < n; i++) {
		sum_x += x[i]; 
		sum_y += y[i];
	}
	float x_ave = sum_x / n;
	float y_ave = sum_y / n;
	for (int i = 0; i != n; i++)
	{
		lxy += (x[i] - x_ave) * (y[i] - y_ave);
		xiSubSqr += (x[i] - x_ave) * (x[i] - x_ave);
	}
	b = lxy / xiSubSqr;
	a = y_ave - b * x_ave;
	printf("y=%0.2fx+%0.2f\n", b, a);
	system("pause");
}