Анализ источника ошибки. Двумерная диффузия

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Diagnostics;
using System.Drawing;
using System.Drawing.Drawing2D;
using System.Drawing.Imaging;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
 
namespace Diffuse2d
{
    public partial class Form1 : Form
    {
        // параметры системы
        Params Params = new Params();
 
        // вычисленные значения системы для каждого тика времени
        double[][,] Values = new double[][,] { };
 
        public Form1()
        {
            InitializeComponent();
            Calculate();
            trackBar1.Maximum = Params.TimeSteps - 1;
            pictureBox1.Image = DrawBitmap(Values[0]);
        }
 
        // Вычисляет итерацию (y_i+1 по y_i)
        double[,] Tick(double[,] y)
        {
            var fx = Params.dt * Params.kx / Params.dx / Params.dx;
            var fy = Params.dt * Params.ky / Params.dy / Params.dy;
            var c = new double[,]
            {
                { 0, fy , 0 },
                { fx, -2*(fx + fy), fx },
                { 0, fy , 0 },
            };
 
            var r = new double[Params.Height, Params.Width];
            for (int row = 1; row < Params.Height - 1; row++)
            {
                for (int col = 1; col < Params.Width - 1; col++)
                {
                    // к старому значению в текущей ячейки внизу
                    r[row, col] = y[row, col];
 
                    // прибавляем соседние ячейки умноженные на коэффициенты из окна c
                    // границы не учитываем
                    for (int i = -1; i <= 1; i++)
                        for (int j = -1; j <= 1; j++)
                            r[row, col] += y[row + i, col + j] * c[1 + i, 1 + j];
                }
            }
 
            return r;
        }
 
        void Calculate()
        {
            Values = new double[Params.TimeSteps][,];
            Values[0] = Params.y0;
            for (int i = 1; i < Params.TimeSteps; i++)
                Values[i] = Tick(Values[i - 1]);
        }
 
        Bitmap DrawBitmap(double[,] values)
        {
            int dl = 40; // ширина и высота (размер) ячейки в пикселях
            var max = values.Cast<double>().Max(); // Максимальное значение ячейки (черный цвет)
            var bmp = new Bitmap((int)(Params.Width * dl), (int)(Params.Height * dl));
            var g = Graphics.FromImage(bmp);
            g.SmoothingMode = SmoothingMode.None; 
 
            for (int i = 0; i < Params.Width; i++)
                for (int j = 0; j < Params.Height; j++)
                {
                    var x = i * dl;
                    var y = j * dl;
 
                    var c = (int)(255.0 * values[i, j] / max);
                    if (c > 255) c = 255;
 
                    g.FillRectangle(
                        new SolidBrush(Color.FromArgb(c, c, c)),
                        (float)x, (float)y, (float)dl, (float)dl);
 
                    g.DrawString($"{values[i, j]:##E+00}", Font, Brushes.Red, (float)x, (float)y);
                }
 
            g.Dispose();
            return bmp;
        }
 
        private void trackBar1_ValueChanged(object sender, EventArgs e)
        {
            var old_image = pictureBox1.Image;
            pictureBox1.Image = DrawBitmap(Values[trackBar1.Value]);
            old_image.Dispose();
        }
 
        private void Form1_Load(object sender, EventArgs e)
        {
 
        }
    }
 
    class Params
    {
        public int TimeSteps { get; set; }
        public int Width => y0.GetLength(0);
        public int Height => y0.GetLength(1);
 
        public double dx { get; set; }
        public double dy { get; set; }
        public double dt { get; set; }
 
        public double kx { get; set; }
        public double ky { get; set; }
 
        public double[,] y0 { get; set; }
 
        public Params()
        {
            TimeSteps = 10000;
            dt = 5;
            dx = 100;
            dy = 100;
            kx = 0.000001;
            ky = 0.000001;
            y0 = new double[20, 20];
            y0[10, 10] = 100.0;
            y0[3, 3] = 50.0;
        }
    }
}