СМЕШАННАЯ ЗАДАЧА ДЛЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

#include <QCoreApplication>
#include <math.h>
#include <iostream>
using namespace std;
int main()
{
    int N, M, i,j;
    double l=2, t=60, T=50, a, r=72 00, c=460, Y=45.4;
    cout<<"Vvod kol-vo uzlov X"<<endl;
    cin>>N;
    double h=l/(N);        // шаг x
    
    a= 1/(2*(Y/(c*r)));
    for (M=1;t/(M*h*h)>0.5;M++);
    cout<<"Vvod kol-vo uzlov T"<<endl<<
          "Kol-vo uzlov >="<<M<<endl;
    cin>>M;
    float p=t/(M);        // шаг т
    double y[N+1], x[N+1];
    cout<<"h="<<h<<endl<<"tau="<<p<<endl;
    //Задание начальных и граничных условий
    for(i=1;i<N;i++)
        y[i]=T*((h*i)/l);
        y[0]=0;
        y[N]=T;
    //Решение
    for(j=1;j<M;j++){
        for(i=1;i<N;i++){
           x[i]=((p/pow(h,2)*y[i-1])+
            ((1-((2*p)/pow(h,2)))*y[i])+
                   (p/pow(h,2))*y[i+1]);}
        x[0]=0;
        x[N]=T;
        for(int k=1;k<N;k++)
            y[k]=x[k];
        }
for (i=0;i<N+1;i++)
    cout<<x[i]<<endl;
    }