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

#define _USE_MATH_DEFINES
#include <iostream>
#include <cmath>
#include <math.h>
#include <fstream>
 
using namespace std;
 
// Задаем начальные условия
double L = 1.0; // Длина стержня
 
const int N = 6; // Количество точек на стержне
const int t = 11; // Количество точек по времени
 
double dx = L / (N - 1); // Шаг по координате
double dt = 6; // Шаг по времени
 
double akvd = 0.116011236; // Коэффициент теплопроводности
double lambda = 413;
 
double r = akvd * dt / (dx * dx); // постоянная
 
 
 
int i, j;
double u[N][t], x[N], T[t];
 
int main()
{
    ofstream outfile;
    outfile.open("znach.txt");
 
   // сетка
    
    for (i = 0; i < N; i++)
    {
        x[i] = i * dx;
    }
    for (j = 0; j < t; j++)
    {
        T[j] = j * dt;
    }
    
 
// нач условия
    u[0][0] = M_PI;
 
    for ( i = 1; i < N; i+=dx)
    {
        u[i][0] = (sin(M_PI * x[i])) / x[i];
    }
    
 
 
    // расчитываем внутр узлы
    
    // Выполняем явную разностную схему
    for (j = 1; j < t; j += dt)
    {
        for (i = 1; i < N ; i += dx)
        {
            u[i][j + 1] = r * (u[i - 1][j] - 2 * u[i][j] + u[i + 1][j]) + u[i][j];
        }
    }
 
    // граничные условия
 
    for ( j = 1; j < t; j += dt)
    {
        u[0][j+1] = (lambda*u[1][j+1] + 50000) / (lambda + 100*dx);
        u[N][j + 1] = 2 * (j+1)*dt;
    }
 
    // вывод в файл
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < t; j++)
        {
            outfile << u[i][j] << ' ';
        }
        outfile << endl;
    }
    outfile.close();
    return 0;
}

double u[N * t]