Уравнение теплопроводности (явная схема)

#include<iostream>
#include<Windows.h>
 
 
using namespace std;
 
double kap(double h,double j);
 
 
int main()
{
    
    const double a = 0.0, b = 1.0, k = 1.0;
    double h,x, kappa = 1.0;
    
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    cout<<"Рассчетная область = [0, 1]"<<endl;
    int N, p=0;;
    double* T0, *T1;
    do{
       cout<<"Введите количество рассчетных узлов: ";
       cin>>N;
       h = (b - a)/N;
    }while (h > 0.02);
 
    T0 = new double [N];
    T1 = new double [N];
 
    const double tau = (k*h*h)/(2*kappa);
    double t_end, t = 0.0;
    do{
        cout<<"Введите рассчетное время: ";
        cin>>t_end;
    }while(t_end > 4);
 
    //выполнение начальных условий
    for(int i=0; i<N; i++)
    {
        x=i*h;
        if(i == 0)
        {   
            T0[i] = 0;      
        }
        else
        {
            if(i == N-1)
            {   
                T0[i] = 0;      
            }
            else
                T0[i]= 1-x; 
        }
    }
    int z = 0;
    
    cout<<z++<<" временной слой: ";
    for(int j = 0; j <= N-1; j++)
    {
        cout<<T0[j]<<' ';
    }
    
    while(t < t_end)
    {
        for(int j = 1; j <= N-2; j++)
            T1[j] = T0[j] + tau / h / h * (kap(h,j+0.5) * (T0[j+1]-T0[j]) - kap(h,j-0.5) * (T0[j]-T0[j-1]));
        T1[N-1] = T1[N-2];
        T1[0]= T1[1];
        cout<<endl; 
        cout<<z++<<" временной слой: ";
        for(int j = 0; j <= N-1; j++)
        {
            cout<<T1[j]<<' ';
        }
        for(int j = 0; j <= N-1; j++)
        {
            T0[j] = T1[j];
        }
        t += tau;
    }
    delete []T0;
    delete []T1;
    cout<<endl;
    return 0;
}
 
 
 
double kap(double h,double j)
{
    double found_kappa=2.0,x2;
    x2=h*j;
    if(x2<0.5)
        found_kappa=1.0;
    return found_kappa;
}