Нужно найти ошибку или исправить в готовой программе

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
 
namespace Ovesyan_Lab_4
{
    public class vector
    {
        public double a, b;
    }
 
    public class Program
    {
        public static void Main(String[] args)
        {
            Console.WriteLine("1.Nelder-Mead:");
            NelderMead();
            Console.WriteLine("\n2.Descent gradient:");
            Descent();
            Console.WriteLine("\n3.Fletcher-Reaves:");
            FletcherReaves();
            Console.WriteLine("\n4.Newton:");
            Newton();
            Console.ReadKey();
        }
 
        public static readonly double e = 0.05;
 
        public static double f(vector v)
        {
            return Math.Pow((v.a - 12), 2) + Math.Pow((v.b - 2), 2);
        }
 
        public static vector[] sort(vector[] v)
        {
            vector x;
 
            for (int i = 0; i < 2; i++)
            {
                for (int j = 0; j < 2; j++)
                {
                    if (f(v[j]) < f(v[j + 1]))
                    {
                        x = v[j + 1];
                        v[j + 1] = v[j];
                        v[j] = x;
                    }
                }
            }
            return v;
        }
 
        public static double[] sven(int i, vector v1, vector v2)
        {
            double x0 = 0.0, d = 0.25;
            double x1 = 0;
            double x2 = 0;
            double x = 0;
            double k = 1;
            double[] segment = new double[2];
            vector fx = new vector();
            vector fx1 = new vector();
            vector fx2 = new vector();
 
            if (i == 0)
            {
                fx.a = v1.a - (x0 * v2.a);
                fx.b = v1.b - (x0 * v2.b);
                fx1.a = v1.a - (x0 - d) * v2.a;
                fx1.b = v1.b - (x0 - d) * v2.b;
                fx2.a = v1.a - (x0 + d) * v2.a;
                fx2.b = v1.b - (x0 + d) * v2.b;
            }
            else
            {
                fx.a = v1.a + (x0 * v2.a);
                fx.b = v1.b + (x0 * v2.b);
                fx1.a = v1.a + (x0 - d) * v2.a;
                fx1.b = v1.b + (x0 - d) * v2.b;
                fx2.a = v1.a + (x0 + d) * v2.a;
                fx2.b = v1.b + (x0 + d) * v2.b;
            }
 
            if ((f(fx1) >= f(fx)) && (f(fx2) >= f(fx)))
            {
                segment[0] = f(fx1);
                segment[1] = f(fx2);
                return segment;
            }
 
            if ((f(fx1) <= f(fx)) && (f(fx2) <= f(fx)))
            {
                Console.WriteLine("Impossible");
                return null;
            }
 
            if ((f(fx1) >= f(fx)) && (f(fx) >= f(fx2)))
            {
                segment[0] = x0;
                d = d;
            }
 
            if ((f(fx1) <= f(fx)) && (f(fx) <= f(fx2)))
            {
                segment[1] = x0;
                d = -d;
            }
 
            x1 = x0 + d;
 
            if (i == 0)
            {
                fx1.a = v1.a - (x1 * v2.a);
                fx1.b = v1.b - (x1 * v2.b);
            }
            else
            {
                fx1.a = v1.a + (x1 * v2.a);
                fx1.b = v1.b + (x1 * v2.b);
            }
            while (true)
            {
                x = x1;
                x1 += Math.Pow(2, k) * d;
                x2 = f(fx1);
 
                if (i == 0)
                {
                    fx1.a = v1.a - (x1 * v2.a);
                    fx1.b = v1.b - (x1 * v2.b);
                }
                else
                {
                    fx1.a = v1.a + (x1 * v2.a);
                    fx1.b = v1.b + (x1 * v2.b);
                }
 
                if (f(fx1) < x2)
                {
                    if (d > 0)
                    {
                        segment[0] = x;
                    }
                    else
                    {
                        segment[1] = x;
                    }
                }
                else
                {
                    break;
                }
 
                k++;
            }
 
            if (d > 0)
            {
                segment[1] = x1;
            }
 
            if (d < 0)
            {
                segment[0] = x1;
            }
 
            return segment;
        }
 
        public static double GoldenSectionSearch(int i, double a, double b, vector v1, vector v2)
        {
 
            int k = 0;
            double x1, x2;
            double eps = 0.01;
            vector fx1 = new vector();
            vector fx2 = new vector();
 
            while ((b - a) > eps)
            {
                x1 = a + ((3 - Math.Sqrt(5)) / 2) * (b - a);
                x2 = a + ((Math.Sqrt(5) - 1) / 2) * (b - a);
 
                if (i == 0)
                {
                    fx1.a = v1.a - x1 * v2.a;
                    fx1.b = v1.b - x1 * v2.b;
                    fx2.a = v1.a - x2 * v2.a;
                    fx2.b = v1.b - x2 * v2.b;
                }
                else
                {
                    fx1.a = v1.a + x1 * v2.a;
                    fx1.b = v1.b + x1 * v2.b;
                    fx2.a = v1.a + x2 * v2.a;
                    fx2.b = v1.b + x2 * v2.b;
                }
 
                if (f(fx1) < f(fx2))
                {
                    b = x2;
                }
                else
                {
                    a = x1;
                }
 
                k++;
            }
 
            return (a + b) / 2;
        }
 
        public static void NelderMead() {
        
        double alpha = 1;
        double beta = 0.5;
        double gamma = 2;
        double d;
        vector x4 = new vector();
        vector x5 = new vector();
        vector x6 = new vector();
        vector temp = new vector();
        vector[] x = new vector[] { new vector(), new vector(), new vector() };
        x[0].a = 8;
        x[0].b = 9;
        x[1].a = 10;
        x[1].b = 11;
        x[2].a = 8;
        x[2].b = 11;
        x = sort(x);
        x4.a = (x[2].a + x[1].a) / 2;
        x4.b = (x[2].b + x[1].b) / 2;
        
        d = Math.Sqrt((Math.Pow((f(x[0]) - f(x4)), 2) +
                       Math.Pow((f(x[1]) - f(x4)), 2) + 
                       Math.Pow((f(x[2]) - f(x4)),2)) / 3);
        
        for (int k = 0; d > e; k++) 
        {
            x = sort(x);
            x4.a = (x[2].a + x[1].a) / 2;
            x4.b = (x[2].b + x[1].b) / 2;
            x5.a = x4.a + alpha * (x4.a - x[0].a);
            x5.b = x4.b + alpha * (x4.b - x[0].b);
            
            d = Math.Sqrt((Math.Pow((f(x[0]) - f(x4)), 2) + 
                           Math.Pow((f(x[1]) - f(x4)), 2) + 
                           Math.Pow((f(x[2]) - f(x4)), 2)) / 3);
            
            if (f(x5) <= f(x[2]))
            {
                x6.a = x4.a + gamma * (x5.a - x4.a);
                x6.b = x4.b + gamma * (x5.b - x4.b);
                
                if (f(x6) < f(x[2]))
                {
                    x[0].a = x6.a;
                    x[0].b = x6.b;
                    continue;
                } 
                else
                {
                    x[0].a = x5.a;
                    x[0].b = x5.b;
                    continue;
                }
            }
            
            if ((f(x[1]) < f(x5)) && (f(x5) <= f(x[0])))
            {
                x[0].a = x4.a + beta * (x[0].a - x4.a);
                x[0].b = x4.b + beta * (x[0].b - x4.b);
                
                continue;
            }
            
            if ((f(x[2]) < f(x5)) && (f(x5) <= f(x[1])))
            {
                x[0].a = x5.a;
                x[0].b = x5.b;
                
                continue;
            }
            
            if (f(x5) > f(x[0])) 
            {
                x[0].a = x[2].a + beta * (x[0].a - x[2].a);
                x[0].b = x[2].b + beta * (x[0].b - x[2].b);
 
                x[1].a = x[2].a + beta * (x[1].a - x[2].a);
                x[1].b = x[2].b + beta * (x[1].b - x[2].b);
            }
 
            Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x[2].a, x[2].b, f(x[2]));
        }
    }
 
        public static void Descent() {
        
        vector x0 = new vector();
        vector x1 = new vector();
        vector df = new vector();
        vector diff = new vector();
        double t, absf, absdif, absdf;
        int M = 20, k = 0;
        x0.a = 12;
        x0.b = 17;
        
        df.a = 2 * (x0.a - 13);
        df.b = 2 * (x0.b - 4);
        
        while (k <= M) 
        {
            double[] segment = sven(0, x0, df);
            t = GoldenSectionSearch(0, segment[0], segment[1], x0, df);
            x1.a = x0.a - (t * df.a);
            x1.b = x0.b - (t * df.b);
 
            Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x1.a, x1.b, f(x1));
            
            diff.a = x1.a - x0.a;
            diff.b = x1.b - x0.b;
            
            absdif = Math.Sqrt(Math.Pow(diff.a, 2) + Math.Pow(diff.b, 2));
            df.a = 2 * (x1.a - 12);
            df.b = 2 * (x1.b - 2);
            
            absf = Math.Sqrt(Math.Pow(x1.a, 2) + Math.Pow(x1.b, 2));
            absdf = Math.Sqrt(Math.Pow(df.a, 2) + Math.Pow(df.b, 2));
            
            if ((f(x0) - f(x1) < e * (1 + Math.Abs(f(x1)))) && (absdif < Math.Sqrt(e) * (absf + 1)) && (absdf <= Math.Pow(e, 1 / 3) * (1 + Math.Abs(f(x1)))))
            {
                break;
            } 
            else 
            {
                x0.a = x1.a;
                x0.b = x1.b;
                df.a = 2 * (x0.a - 13);
                df.b = 2 * (x0.b - 4);
                k++;
            }
        }
    }
 
        public static void FletcherReaves() {
        
        vector diff = new vector();
        vector preX = new vector();
        vector x0 = new vector();
        vector x1 = new vector();
        vector prevdf = new vector();
        vector df = new vector();
        vector d = new vector();
        double t, absdf, absdif, absprevdf = 0, beta;
        int M = 20, k = 0, i = 0;
        x0.a = 2;
        x0.b = 9;
        
        df.a = 2 * (x0.a - 12);
        df.b = 2 * (x0.b - 2);
        absdf = Math.Sqrt(Math.Pow(df.a, 2) + Math.Pow(df.b, 2));
        
        while (k <= M)
        {
            if (absdf < e)
            {
                Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x1.a, x1.b, f(x1));
                break;
            }
            
            if (k == 0) 
            {
                d.a -= df.a;
                d.b -= df.b;
            } 
            else 
            {
                beta = Math.Pow((absdf / absprevdf), 2);
                d.a = (d.a * beta) - df.a;
                d.b = (d.b * beta) - df.b;
            }
            
            double[] segment = sven(1, x0, d);
            t = GoldenSectionSearch(1, segment[0], segment[1], x0, d);
            x1.a = x0.a + (t * d.a);
            x1.b = x0.b + (t * d.b);
            Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x1.a, x1.b, f(x1));
            diff.a = x1.a - x0.a;
            diff.b = x1.b - x0.b;
            absdif = Math.Sqrt(Math.Pow(diff.a, 2) + Math.Pow(diff.b, 2));
            
            if ((absdif < e) && (Math.Abs(f(x1) - f(x0)) < e)) 
            {
                i++;
            } 
            else 
            {
                i = 0;
                preX = x0;
                prevdf.a = 2 * (x0.a - 12);
                prevdf.b = 2 * (x0.b - 2);
                absprevdf = absdf;
                x0.a = x1.a;
                x0.b = x1.b;
                df.a = 2 * (x0.a - 12);
                df.b = 2 * (x0.b - 2);
                absdf = Math.Sqrt(Math.Pow(df.a, 2) + Math.Pow(df.b, 2));
                k++;
                
                continue;
            }
            
            if (i == 2)
            {
                break;
            }
            
            i++;
        }
 
    }
 
        public static void Newton() {
        
        int M = 20;
        int k = 0, i = 0;
        vector x0 = new vector();
        vector x1 = new vector();
        vector df = new vector();
        vector d = new vector();
        vector diff = new vector();
        double absdf, absdif;
        double t = 0;
        x0.a = 2;
        x0.b = 9;
        df.a = 2 * (x0.a - 13);
        df.b = 2 * (x0.b - 4);
        absdf = Math.Sqrt(Math.Pow(df.a, 2) + Math.Pow(df.b, 2));
        
        while (k <= M) 
        {
            
            double[] invmatrix = { 0.5, 0, 0, 0.5 };
            double det = invmatrix[0] * invmatrix[3] - invmatrix[1] * invmatrix[2];
            
            if (absdf < e)
            {
                Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x1.a, x1.b, f(x1));
                break;
            }
            
            if (invmatrix[0] > 0 && det > 0) 
            {
                d.a = -((invmatrix[0] * df.a) + (invmatrix[1] * df.b));
                d.b = -((invmatrix[2] * df.a) + (invmatrix[3] * df.b));
                t = 1;
            } 
            else 
            {
                d.a = -d.a;
                d.b = -d.b;
            }
            
            x1.a = x0.a + (d.a * t);
            x1.b = x0.b + (d.b * t);
            diff.a = x1.a - x0.a;
            diff.b = x1.b - x0.b;
            absdif = Math.Sqrt(Math.Pow(diff.a, 2) + Math.Pow(diff.b, 2));
 
            Console.WriteLine("Iteration {0}: Current Point ({1:F3},{2:F3}) f(xk)={3:F3}", k, x1.a, x1.b, f(x1));
            
            if ((absdif < e) && (Math.Abs(f(x1) - f(x0)) < e))
            {
                i++;
            } 
            else
            {
                i = 0;
                x0.a = x1.a;
                x0.b = x1.b;
                df.a = 2 * (x0.a - 13);
                df.b = 2 * (x0.b - 4);
                absdf = Math.Sqrt(Math.Pow(df.a, 2) + Math.Pow(df.b, 2));
                k++;
                
                continue;
            }
            
            if (i == 2)
            {
                break;
            }
            
            i++;
        }
    }
    }
}