Не работает код решения СНАУ методом Ньютона

using System;
using System.Collections.Generic;
using System.Text;
 
namespace Sistema
{
    class Program
    {
        static void F(double[] x, ref double[] y)//уравнения
        {
            y[0] = -Math.Sin(x[1]) - 0.4;
            y[1] = 0.5*(Math.Cos(x[0] + 1)) ;
            
           
        }
        static void get_matrix(double[][] a, ref double[][] b, int n, int i)
        {
 
            int di = 0, dj = 0;
            for (int ki = 0; ki < n - 1; ki++)
            {
                if (ki == i) di = 1;
                dj = 0;
                for (int kj = 0; kj < n - 1; kj++)
                {
                    if (kj == 0) dj = 1;
                    b[ki][kj] = a[ki + di][kj + dj];
                }
            }
        }
 
        static double Det(double[][] a, int n)//детерминант матрицы а
        {
            double d = 0;//определитель
            int k = 1;
            double[][] b = new double[n][]; //массив указателей на строки
            for (int i = 0; i < n; i++) //выделение памяти под каждую строку массива
                b[i] = new double[n];
            if (n == 1) d = a[0][0];
            else if (n == 2) d = a[0][0] * a[1][1] - a[1][0] * a[0][1];
            else
                for (int i = 0; i < n; i++)
                {
                    get_matrix(a, ref b, n, i);
                    d += k * a[i][0] * Det(b, n - 1);
                    k = -k;
                }
            return d;
        }
 
        static void Jacobian(double[] x, ref double[][] a)
        {
            double det;
            //забиваем матрицу частными производными            
            a[0][0] = 0; a[0][1] = -Math.Cos(x[1]);
            a[1][0] =-0.5*(Math.Sin(x[0]+1)); a[1][1] = 0;
            
            det = -Math.Cos(x[1])*0.5*Math.Sin(x[0]+1);
            for (int i = 0; i < 2; i++)
                for (int j = 0; j < 2; j++) a[i][j] = a[i][j] / det;
       
        }
 
        static void Newthon(ref double[] x, double eps, ref int k)
        {
            bool cod;
            double[] x0 = new double[2];
            double[] y = new double[2];
            double[][] a = new double[2][];
            for (int i = 0; i < 2; i++) a[i] = new double[2];
            do
            {
                for (int i = 0; i < 2; i++) x0[i] = x[i];
                F(x, ref y);
                Jacobian(x, ref a);
                for (int i = 0; i < 2; i++)
                    for (int j = 0; j < 2; j++) x[i] = x[i] - a[i][j] * y[j];
                cod = true;
                for (int i = 0; i < 2; i++)
                    if (Math.Abs(x[i] - x0[i]) > eps) cod = false;
                k++;
            } while (cod == false);
 
          
 
        }
 
        static void Main(string[] args)
        {
            double eps = 0.001;
            double[] x = new double[2];
            x[0] = 1; x[1] = 1; 
            int k = 0;
            Newthon(ref x, eps, ref k);
            Console.WriteLine("X:");
            for (int i = 0; i < 2; i++) Console.WriteLine(x[i]);//выводим иксы...
            Console.WriteLine(k);
            double[] y = new double[2];
            F(x, ref y);
            Console.WriteLine("Y:");
          
            for (int i = 0; i < 2; i++) Console.WriteLine(y[i]);//это для проверки...
           
            Console.ReadKey();
 
        }
    }
}