Решение уравнения переноса

#include <iostream>
#include <fstream> 
#include <math.h>
using namespace std;
 
 
double fun(double x)
{
    if (x >= -1 && x <= 1)
        return sin(3.14*x);
    else
        return 0;
}
 
 
int main()
{
    ofstream file;
    file.open("0.3.txt");
    setlocale(LC_ALL, "Russian");
    double a = -1.0;//a область значений аргумента
    double T = 1.0;// шаг по времени
    const int N = 100, M = 100; // число разбиений
    double tau = T / M;
    double h = 0.003; // шаг
    double Q = 0.3; //Число куранта 
 
    double u[N][M];
 
    //  i = (x); j = (t)
 
    for (int i = 0; i < N; i++)
 
        u[i][0] = fun(a - i * h); //при t=0  
 
    for (int j = 0; j < M; j++)
        u[0][j] = 0;  //краевое условие 
 
    for (int j = 0; j < M; j++)
 
        u[1][j] = 0;  //краевое условие
    
 
    for (int j = 0; j < M - 1; j++)
 
        for (int i = 1; i < N; i++)
 
            u[i][j + 1] = Q * u[i - 1][j] + (1 - Q)*u[i][j];
 
    for (int j = M - 1; j >= 0; j--)
 
        for (int i = 0; i < N; i++)
            file << u[i][j] << endl;
 
    file.close();
    system("pause");
    return 0;
}