МНК, повреждение кучи

#include "stdafx.h"
#include "math.h"
#include "windows.h"
#include <iostream>
#include <fstream>
 
using namespace std;
 
char buffer[100];           
int n, m;
 
void main()
{
    AnsiToOem("Введи количество исходных точек: ", buffer);
    cout << buffer; cin >> n; cout << endl << endl;
    AnsiToOem("Введи желаемую степень полинома: ", buffer);
    cout << buffer; cin >> m; cout << endl << endl;
 
    long double *x=new long double [m+1], *y=new long double [m+1];
    long double **b=new long double *[m+1];
        for(int i=0; i<m+1; i++)
        b[i]=new long double[m+1];
    long double *c=new long double [m+1];
    long double *cr=new long double [m+1];
    
    ifstream inputfile("input.txt");
    if(!inputfile) 
    {
        AnsiToOem("Где входной файл??? >_<", buffer);
        cerr << buffer << endl << endl;
        AnsiToOem("Положи input.txt в одну директорию с этим exe-файлом..", buffer);
        cerr << buffer << endl << endl;
        system ("PAUSE");
        exit(1);
    }
    for(int i=0; i<n; i++)
    {
        inputfile >> x[i];
        inputfile >> y[i];
    }
 
    for(int k=0; k<m+1; k++)
    {
        for(int j=0; j<m+1; j++)
        {
            b[k][j]=0;
            for (int i=0; i<n; i++)
                b[k][j]+=pow(x[i], j+k);
        }
        c[k]=0;
        for(int i=0; i<n; i++)
            c[k]+=pow(x[i],k)*y[i];
    }
 
    for (int k=0; k<m; k++)
    {
        int imax=k;
        long double mx=b[k][k];
        for(int i=k; i<m+1; i++)        
            if(fabs(mx)<fabs(b[i][k]))
            {
                mx=b[i][k];
                imax=i;
            }
 
        mx=c[k];        
        c[k]=c[imax];
        c[imax]=mx;
 
        for(int j=k; j<m+1; j++)
        {
            mx=b[k][j];     
            b[k][j]=b[imax][j];
            b[imax][j]=mx;
        }
        
        for(int i=k+1; i<m+1; i++)          
        {
            for(int j=k; j<m+1; j++)
                b[i][j]-=b[k][j]*b[i][k]/b[k][k];
            c[i]-=c[k]*b[i][k]/b[k][k];
            b[i][k]=0;
        }
    }
 
    for(int i=0; i<m+1; i++)                    
        cr[i]=c[i];
    for(int i=m; i>=0; i--)
    {
        for(int j=i+1; j<m+1; j++)
            cr[i]-=b[i][j]*cr[j];
        cr[i]=cr[i]/b[i][i];
    }
        
    AnsiToOem("Найденные решения СЛАУ:", buffer);
    cout << endl << endl << buffer << endl << endl;
    for(int i=0; i<m+1; i++)
    {
        AnsiToOem("Корень №", buffer);
        cout << buffer << i+1 << ": " << cr[i] << endl;
    }
    cout << endl << endl;
    system ("PAUSE");
    system ("cls");
 
    delete[] x,y,b,c,cr; 
 
}