[C#] Метод конечных разностей

[ardor]

Помогите разобраться с задачей. Нужно написать программу для расчета теплопроводности стержния методом конечных разностей. Схема Ричардсона.
Вот что у меня:

public void MKR()
        {
            double T0 = /*Convert.ToDouble(textBox1.Text);*/ 20;
            double TL = /*Convert.ToDouble(textBox2.Text);*/ 300;
            double TR = /*Convert.ToDouble(textBox3.Text);*/ 100;
            double EndTime = Convert.ToDouble(textBox5.Text);
            double L = Convert.ToDouble(textBox6.Text);
            double N = /*Convert.ToDouble(textBox4.Text);*/ 50;
            double lambda = 46;
            double ro = 7800;
            double c = 460;
            double a = lambda / (ro * c);
            double h = L/(N-1);
            double tau = (0.25 * h * h) / a;
            double nT = EndTime/tau;
            int nTi = Convert.ToInt32(nT);
            int Ni = Convert.ToInt32(N);
            T = new double[nTi+2, Ni+1];

            for (int i = 1; i < N - 1; i++)
            {
                T[0, i] = T0;
            }
            for (int i = 0; i < Convert.ToInt32(nT); i++)
            {
                T[i, 0] = TL;
                T[i, Ni-1] = TR;
            }
            int k = 1;
            double time = 0;
            while (time < EndTime)
            {
                time += tau;
                for (int i = 1; i < Ni-1; i++)
                {
                    //T[k, i] = T[k - 1, i] + a * tau / (h * h) * (T[k - 1, i + 1] - 2 * T[k - 1, i] + T[k - 1, i - 1]);
                    if (k == 1)
                    {
                        T[k, i] = T[k - 1, i] + a * tau / (h * h) * (T[k - 1, i + 1] - 2 * T[k - 1, i] + T[k - 1, i - 1]);
                    }
                    else
                    T[k, i] = T[k - 2, i] + 2 * a * tau / (h * h) * (T[k - 1, i - 1] - 2 * T[k - 1, i] + T[k - 1, i + 1]);
                }
                k++;
            }
dataGridView1.ColumnCount = Ni;
            dataGridView1.RowCount = nTi;
            n1 = Ni;
            n2 = nTi+1;

            for (int i = 0; i < Ni; i++)
            {
                for (int j = 0; j < nTi; j++)
                {
                    dataGridView1[i, j].Value = T[j, i].ToString();
                }
            }        
        }

Значения жутко скачут, да и выходят слишком маленькими. Решал эту же задачу с помощью явной четырехточечной схемы. Все работает. С этой схемой никак немогу разобраться =(