求解输运方程时的代码错误

更新了问题的内容，因为我自己理解并纠正了一些地方。因此有以下变量：

hx (шаг по пространству)= 0.1

ht (шаг по времени) = 0.5

Nx (Количество шагов по пространству) = 10;

Nt (Количество шагов по времени) = 12;

wx[i] - массив содержит все шаги по пространству

wt[j] - массив содержит все шаги по времени

wht[j][i]-двумерный массив по которому будет строится результирующий график;

于是发现了几个问题：

1）在代码中找到

for(int i = 0; i < Nx; i++)        
{
    wx[i+1] = wx[i] + hx;      //массив 
    wht[0][i] = fn(T, i * hx); //i * hx//Исправлено.
}

它在哪里fn：

//Функция Хэвисайда - Начальное условие(Граничное условие) а начальное 0-1
public:static double fn(int T, double x)
{
    if (x >= 0) 
        return T;
    else if(x < 0)
        return 0;  
}

我设置的初始条件不正确

就像图表上的酒一样，在耻辱中，一般解决方案“wave”在某处可见，但其余部分是错误的，因为我在代码中错误地设置了初始条件

2）循环和数组，即：

调用了数组中不存在的元素，结果出现了可怕的数字和图表中的跳跃，但我无法更改索引，因为存在特定的公式。

for(int j = 0;j<Nt;j++)
{
    for(int i = 1;i<Nx-1;i++)///Исправлено с помощью ответа пользователя Denis Protopopov(и график немного изменился)
    {
        wht[j+1][i] = 
            ((wht[j + 1][i] - wht[j][i]) / ht) + a * ((wht[j][i+1] - wht[j][i]) / hx);
        wht[j+1][i] = -a * (ht * (wht[j][i+1] + wht[j][i]) / hx) + (wht[j][i]);
    }
}

如果我写而不是i = 0,i = 2或j = 2，那么图形的绘制将很糟糕。

所有代码：

public:static double qx0(double x)//ось пространства
           {
               if (x<=0) 
                   return 0;
               else
                   return x;
           }

    public:static double qt0(double t)//ось по времени
           {
               if (t<=0) 
                   return 0;
               else
                   return t;
           }

    public:static double fn(int T,double x)//Функция Хэвисайда - Начальное условие
           {
               if (x>=0)
                   return T;
               else if(x<0)
                   return 0;  
           }

public: void drawfunc(double T)
            {
                double xmin = -5;
                double xmax = 10;
                for(double x = xmin;x<xmax;x+=0.01)
                {       
            chart1->Series["Функция Хэвисайда"]->BorderWidth=3;
            chart1->Series["Функция Хэвисайда"]->Points->AddXY(x,fn(T,x));                       
                }                                       
            }


    private: System::Void button1_Click(System::Object^  sender, System::EventArgs^  e) 
             { 
                 ///-------Переменные для работы с разностными схемами
                 double a;//скорость переноса
                 double hx,ht;//шаги по пространству и времени
                 int Nx,Nt;//На сколько частей разбивае отрезок сетки
                 int T;//Входной параметр для искомой функции(Хэвисайд)
                 double wx[20]={0};//массив точек x
                 double wt[20]={0};//массив точек t
                 double wht[20][20]={0};// массив для разностной схемы(сетки) 
                 //double res[20][20]={0};//результирующий массив

                 if(textBox1->Text!="" && textBox2->Text!="" && textBox3->Text!="" && textBox4->Text!="" && textBox5->Text!="" && textBox6->Text!="" && textBox6->Text!="")
                 {

                     ///-----Ввод данных-----
                     a=Convert::ToDouble(textBox1->Text);
                     hx=Convert::ToDouble(textBox2->Text);
                     q=Convert::ToDouble(textBox3->Text);
                     ht=Convert::ToDouble(textBox4->Text);
                     Nx=Convert::ToInt32(textBox5->Text);
                     Nt=Convert::ToInt32(textBox6->Text);
                     T=Convert::ToInt32(textBox7->Text);

                     //----Построение сетки,узлов
                     for(int i = 0;i<Nx;i++)        
                     {
                         wx[i+1]=wx[i]+hx;//массив 
                         wht[0][i]=qx0(wx[i]);                       
                     }   

                     for(int j = 0;j<Nt;j++) //
                     {                                                       
                         wt[j+1]=wt[j]+ht;
                         wht[j+1][0] = qt0(wt[j+1]);                     
                     }


                     ///-----------------Вычисление основых формул с разностной схемой

                     for(int j = 0;j<Nt;j++)
                     {
                         for(int i = 1;i<Nx-1;i++)///Вывод формул wht[j+1][i] 
                         {                           
        wht[j+1][i]=((wht[j+1][i]-wht[j][i])/ht)+a*((wht[j][i+1]-wht[j][i])/hx);
        wht[j+1][i]=-a*(ht*(wht[j][i+1] + wht[j][i])/hx) + (wht[j][i])+ht*fn(wt[j],hx); 
                         }
                     }

                     //----Для записи в простой текстовый файл                  
                     String^ fileName = "results.txt";
                     StreamWriter^ sw = gcnew StreamWriter(fileName);
                     for(int j = 0;j<Nt;j++)    
                     {
                         for(int i = 0 ;i<Nx;i++)
                         {
                             sw->Write("{0} ",wht[j][i]);
                         }
                         sw->WriteLine();
                     }
                     sw->Close(); 


                     //Тестовые функции для отображения
                     drawfunc(T);// Вызов функции для рисования

                     ///---------------Построение графика                        

                              for(int j = 0;j<Nt;j++)
                     {
                          for(int i = 0;i<Nx-1;i++)
                         {  
                        chart2->Series["Series2"]->BorderWidth=3;           
                        chart2->Series["Series2"]->Points->AddXY(i,wht[j][i]);
                         }
                     }                                           
                 }
                 else MessageBox::Show("Ошибка!Не все поля заполнены!");
             }