Массивы. Программа на поиск минимального Хи квадрата и вывода соответствующей модели

#include "stdafx.h"
#include <fstream>
#include <iostream>
#include <math.h>
#include <string>
using namespace std;
double pow(double,int);
double sqrt(double);
 
int main()
{
    int n1,n2,i,j,g;
    cout<<"Vvedite kollichestvo strok: n1="<<endl;
    cin>>n1;
    cout<<"Vvedite kollichestvo modelei: n2="<<endl;
    cin>>n2;
    double a[42002],b[42002][100],c[100];
    char d[100][20];
    char e[20];
    ifstream infile;
    //считывание данных с файла, содержащий интерполированные модели, и занесение в массив b[i][j] хи квадрата
    infile.open("C:\\Users\\Vita\\Desktop\\a.txt");
    for (i=0;i<=n1;i++)
    {
        for (j=0;j<=n2;j++)
        {
            infile >> a[j];
        }
        for (j=0;j<n2;j++)
        {
            b[i][j]=pow(a[n2]-a[j],2);
        }
    }
    infile.close();
    //сумма каждого столбца в массиве хи квадрат и занесение этих сумм в массив с[j]
    for (j=0;j<=n2;j++)
    {
        for (i=0;i<=n1;i++)
        {
            c[j]+=b[i][j];
        }
    }
    // поиск наименьшнго элемента из массива сумм хи квадратов с(j), и вывод строки минимального элемента g
    double minimal=c[0];
    for (i=0;i<=n2;i++)
    {
      if (minimal>c[i])
      {
        minimal=c[i];
        g=i;
      }
      h[i]=0.0;
    }
 
    //считывание сетки по моделям для отождествления (набить самостоятельно заранее)
    infile.open("C:\\Users\\Vita\\Desktop\\b.txt");
    for (i=0;i<=n2;i++)
    {
        infile >> e;
        strcpy(d[i], e);
    }
    //отождествление минимального хи квадрат с названием модели и вывод в файл itog
    ofstream outfile;
    outfile.open("C:\\Users\\Vita\\Desktop\\itog.dat");
    outfile<<d[g];
    outfile.close();
    return 0;
}

double a[42002],b[42002]