Алгоритм Свена и метод дихотомии для двумерного случая

 static double f(double[] x)
        {
            return Math.Pow(x[0] - 4, 2) + Math.Pow(x[1] - 5, 2);
        }
 
        static double f(double x1, double x2)
        {
            return Math.Pow(x1 - 4, 2) + Math.Pow(x2 - 5, 2);
        }
public static double Bisect(double[] x, double[] S)
        {
            double a = -1, b = 1;
            Svenns(x, Solution.E, ref a, ref b);
            double delta = E / 30;
            double na = 0, nb = 0;
            do
            {
                na = ((a + b) / 2) - delta;
                nb = ((a + b) / 2) + delta;
 
                if (f(x[0] + na * S[0], x[1] + na * S[1]) < f(x[0] + nb * S[0], x[1] + nb * S[1]))
                { b = nb; }
                else
                { a = na; }
            }
            while (Math.Abs(b - a) > E);
 
            return (a + b) / 2.0;
        }
 
public static void Svenns(double x0, double t, ref double l, ref double r)
        {
            double delta = 0, x = 0, xn = 0;
            int k = 0;
 
            if (f(x0 - t) >= f(x0) && f(x0 + t) >= f(x0))
            {
                l = f(x0 - t); r = f(x0 + t);
            }
            else if (f(x0) >= f(x0 - t) && f(x0) >= f(x0 + t))
            {
                Console.WriteLine("Случайным образом будет выбрана другая точка.");
                Random rnd = new Random();
                x0 = rnd.Next(-100, 100);
            }
            else
            {
                if (f(x0 - t) >= f(x0) && f(x0) >= f(x0 + t))
                {
                    delta = t;
                    l = x0;
                    x = x0 + t;
                    k = 1;
                }
 
                if (f(x0) >= f(x0 - t) && f(x0 + t) >= f(x0))
                {
                    delta = -t;
                    r = x0;
                    x = x0 - t;
                    k = 1;
                }
 
                xnext = x + Math.Pow(2, k) * delta;
                while (f(xnext) < f(x))
                {
                    xnext = x + Math.Pow(2, k) * delta;
                    if (f(xnext) < f(x) && delta == t)
                    {
                        l = x;
                        r = xnext;
                    }
                    if (f(xnext) < f(x) && delta == -t)
                    {
                        r = x;
                        l = xnext;
                    }
                    k++;
                }
            }
        }

        public Vector Solve()
        {
            Vector previousIterationY;
            var lambda = new Vector(n);
            do
            {
                previousIterationY = y;                
                for (int j = 0; j < n; j++)
                {
                    var previousStepY = y;
                    int i = j;
                    lambda[j] = Dichotomy(x => LambdaSolve(x, i), y[j]);
                    y += lambda[j] * d[j];
                }
            } while ((y - previousIterationY).Norm >= epsilon);
            return y;
        }

public static double fst(double[,] x, double[,] idx)
        {
            double[] S = new double[2];
            double[,] s = (double[,])idx.Clone();
            double[,] xprev = new double[2, 1];
            
            double lambda = 0;
            int k = 0;
            do
            {
                xprev = (double[,])x.Clone();
 
                lambda = Bisect(0, 50, new double[] { x[0, 0], x[1, 0] }, new double[2] { s[0, 0], s[1, 0] });
 
                x[0, 0] = xprev[0, 0] + lambda * s[0, 0];
                x[1, 0] = xprev[1, 0] + lambda * s[1, 0];
                k++;
            }
            while (Matrix.Euclid(grd(x)) >= E);
 
            return lambda;
        }