Метод Предиктор-корректор для решения задачи со стержнем

#include <iostream>
#include <math.h>
#include <iostream>
#include <fstream>
#include <iomanip>
 
using namespace std;
const int ro=2712, c=5100,N=100;//c=5100
double mass[N][N],v[N][N],sigma[N][N];
double h,tau,l,t;
int i,j,n,k;
int main(){
    cin>>l>>n;
    h=l/n;
    cout<<"h="<<h;
    tau=c*h;
    cout<<"tau="<<tau;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++)
        {
            v[0][j]=0;
            sigma[0][j]=0;
            sigma[i][0]=8;
            v[i][n]=0;
        }
    }
    for(i=0;i<n-3;i++){
        for(j=0;j<n-2;j++)
        {
        if(i%2==0&&j%2==0){
             sigma[i+2][j+1]=ro*c*(v[i + 3][j]-v[i + 1][j])+sigma[i+3][j]+sigma[i+1][j];
      sigma[i+2][j+1]=sigma[i+2][j+1]/2;
      v[i+2][j+1]=(v[i+3][j]+v[i+1][j])/2+(sigma[i+3][j]-sigma[i+1][j])/2*ro*c;
      v[i+1][j+2]=(tau*(sigma[i+2][j+1]-sigma[i][j+1]))/(ro*h)+v[i+1][j];
      sigma[i+1][j+2]=(tau*ro*c*c/h)*(v[i+2][j+1]-v[i][j+1])+sigma[i+1][j];
        }
        }
    }
    cout<<"sigma"<<endl;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++)
        {
            cout<<sigma[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<endl;
    cout<<"V"<<endl;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++)
        {
            cout<<v[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<endl;
ofstream out("rez.txt");
ofstream out2("rez2.txt");
for(i=0;i<n;i++){
        for(j=0;j<n;j++)
        {
        out<<v[i][j]<<" ";
        
        }
        out<<endl;
}
for(i=0;i<n;i++){
        for(j=0;j<n;j++)
        {
        out2<<sigma[i][j]<<" ";
        
        }
        out2<<endl;
}
 
 
    system("pause");
    return 0;
}