Ошибка при работе с динамическим массивом

int main( void )
{
    
    int i,j;
    int n=2;
    long double S, Sn, e=0.000001;
    S=0;
    //double mas[n] [n + 1];
        
    for (;;)
    {
        long double**mass=new long double*[n];
        for(int i=0; i<=n;i++)  
            mass[i]=new long double[2]; 
/*      
        char* nfile;
        itoa (n, nfile, 16);
//другой способ:
std::stringstream sstrm;
sstrm << "D:\\Tab" << nfile << ".txt";
std::string str = sstrm.str();
 
//и еще один способ, не работающий:
        
        string str=string("D:\\Tab")+ string(nfile) + string(".txt");
 
*/
        ostringstream ostr;
        ostr << "D:\\data\\" << n << ".txt";
        string str = ostr.str();/**/
 
 
 
        ifstream input_file(str);
        for(i=0;i<n;i++)
            for(j=0;j<2;j++)
                input_file>>mass[i][j];
        input_file.close ();
 
        Sn=0;
        for (i=0; i<n; i++)
            for (j=0; j<n; j++)
                Sn+=mass[i][1]*mass[j][1]*exp(-2*mass[i][0]*mass[j][0]);
        
        cout << "Iteration N_" << n << "    I=" << Sn<<"\n";
        
        for (i=0;i<n;i++)
            delete[]mass[i];
        delete[]mass;
        if (fabs(S-Sn)<e)
        {
            break;
        }
        S=Sn;   
        n=n+1;
    }
    return (0);
}

for(int i=0; i<=n;i++)