[C言語]最小二乗法を用いて数列を一次関数に近似してgnuplotで表示してみた[コード付き]

//数列を生成しファイルに入力する
//ave,verは雑音における平均と分散
void CRE1(int N, double x1, double x2, double ave, double ver){
    FILE *fp;
      
    if((fp=fopen("kk.txt","w"))==NULL){
        printf("ファイルに書き込めません");
    }else{
        for(int i=1; i<=N; i++){
              
            double Y=0;
            Y=x1+x2*((double)i/N)+ave+ver*BOX();//1次式にガウス性雑音を付加
          
            fprintf(fp,"%f\n",Y);
              
        }
        fclose(fp);
    }
}

Y=x1+x2*((double)i/N)+ave+ver*BOX();

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
//ボックスミュラー法によるガウス性雑音の生成
double BOX(void){
    double y1,y2,x1;
    y1 = (double)rand()/(double)RAND_MAX;
    y2 = (double)rand()/(double)RAND_MAX;
              
    x1 = sqrt(-2*log(y1))*cos(2*M_PI*y2);
    
    return x1;
}

//数列を生成しファイルに入力する
//ave,verは雑音における平均と分散
void CRE1(int N, double x1, double x2, double ave, double ver){
    FILE *fp;
      
    if((fp=fopen("kk.txt","w"))==NULL){
        printf("ファイルに書き込めません");
    }else{
        for(int i=1; i<=N; i++){
              
            double Y=0;
            Y=x1+x2*((double)i/N)+ave+ver*BOX();//1次式にガウス性雑音を付加
          
            fprintf(fp,"%f\n",Y);
              
        }
        fclose(fp);
    }
    
}


//正規方程式で用いる行列Aの作成　1次式
void CREA1(int N, double A[N][2]){
    double t = 1/(double)N;
    for(int i=0; i<N; i++){
        A[i][0]=1;
        A[i][1]=t*(i+1);
    }
}


//数列が入ったファイルからプログラム内の配列に出力する
void SCAN1(int N, double B[N]){
    FILE *fp;
    
    if((fp=fopen("kk.txt","r"))==NULL){
        printf("ファイルを開けません");
    }else{
        for(int i=0; i<N; i++){
            double b=0;
            fscanf(fp,"%lf",&b);
            B[i] = b;
        }
        fclose(fp);
    }
}

//正規方程式の左辺の計算
void CALA(int N, int c, double A[N][c], double L[c][c]){
    //転置したAと普通のAの積
    for(int n=0; n<c; n++){
        for(int m=0; m<c; m++){
            double S=0;
            for(int k=0; k<N; k++){
                S = S + A[k][n]*A[k][m];
            }
                L[n][m] = S; //Lに結果を入れる
          
        }
    }
}

//正規方程式の右辺の計算
void CALB(int N, int c, double A[N][c], double B[N], double M[c]){
    for(int i=0; i<c; i++){
        double S=0;
        for(int k=0; k<N; k++){
            S = S + A[k][i]*B[k];
        }
        M[i] = S; //Mに結果を入れる
    }
}

void INP(int c, double L[c][c], double M[c]){
    int i,j;
    FILE *fp;
    
    if((fp = fopen("ans1.txt","w"))==NULL)//ここのファイルは何でも良い
        printf("ファイルを開けません\n");
    else{
        for(i=0; i<c; i++){
            for(j=0; j<c; j++){
                fprintf(fp,"%f ",L[i][j]);
            }
            fprintf(fp,"%f\n",M[i]);
        }
        fclose(fp);
    }
}

int main(void) {
   int N;
    int x1,x2;
    printf("発生する点の数を入力してください:");
    scanf("%d",&N);
    printf("近似する一次式の係数を入力してください\nx1:");
    scanf("%d",&x1);
    printf("\nx2:");
    scanf("%d",&x2);
    
    double ave,ver;
    printf("発生する雑音の平均を入力してください:");
    scanf("%lf",&ave);
    printf("発生する雑音の分散を入力してください:");
    scanf("%lf",&ver);
    
    double A[N][2];
    double B[N];
    
    CRE1(N,x1,x2,ave,ver);
    SCAN1(N,B);
    CREA1(N,A);
    
    double L[2][2];
    double M[2];
       
    CALA(N,2,A,L);
    CALB(N,2,A,B,M);
       
    INP(2,L,M);
    
    return 0;
}