最小二乘法-直线 拟合 -C代码实现2

/*
    本实验根据数组x[], y[]列出的一组数据，用最小二乘法求它的拟合直线。
*/
#include <stdio.h>

#define maxn 5

int main(){
    double x[maxn] = {165, 123, 150, 123, 141};
    double y[maxn] = {187, 126, 172, 125, 148};
    double xi = 0, x2 = 0, yi = 0, xy = 0, a0, a1;

    for(int i = 0; i < maxn; i++){
        xi += x[i], x2 += x[i] * x[i], yi += y[i], xy += x[i] * y[i];
    }
    a0 = (yi * x2 - xy * xi) / (x2 * maxn - xi * xi);
    a1 = (yi * xi - xy * maxn) / (xi * xi - x2 * maxn);
    printf("P(x) = %f%+fx\n", a0, a1);

    return 0;
}