Метод Ньютона, зацикливаеться программа

[mishammm]

Здравствуйте. У меня стоит задача нахождение минимум функции методом Ньютона с дроблением шага:
Вот моя программа:

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace lr6
{
    class Program
    {

        public static double f(double x1, double x2)
        {
            return Math.Sqrt(Math.Pow(x1, 2) + Math.Pow(x2, 2) + 1) + 0.5 * x1 - 0.5 * x2;
        }

        public static double dx1(double x1, double x2)
        {
            return x1 / (Math.Sqrt(Math.Pow(x1, 2) + Math.Pow(x2, 2))+0.5);
            
        }

        public static double dx2(double x1, double x2)
        {
            return x2 / (Math.Sqrt(Math.Pow(x1, 2) + Math.Pow(x2, 2)) - 0.5);
            
        }

        public static double dx1x1(double x1, double x2)
        {
            return (-1 * (x1 * x1) / Math.Pow((x1 * x1 + x2 * x2 + 1), 3 / 2)) + 1 / (Math.Pow((x1 * x1 + x2 * x2 + 1), 1 / 2));
           
        }

        public static double dx1x2(double x1, double x2)
        {
            return -((x1) * x2) / (Math.Pow((x1 * x1 + x2 * x2 + 1), 3 / 2)); 
            
        }

        public static double dx2x2(double x1, double x2)
        {
            return (-1 * (x2 * x2) / Math.Pow((x1 * x1 + x2 * x2 + 1), 3 / 2)) + 1 / (Math.Pow((x1 * x1 + x2 * x2 + 1), 1 / 2));
        }

        static void Main(string[] args)
        {

            double x1 = 1, x2 = 1;
            double p1 = 0, p2 = 0;
            double alpha = 0.5, eps = 0.25, beta = 0.5, delta = 0.001;
            int k = 0;

            double a11, a12, a21, a22;
            double gr1, gr2;

            while (Math.Abs(dx1(x1, x2)) > delta && Math.Abs(dx2(x1, x2)) > delta)
            {
                a11 = dx1x1(x1, x2); a12 = dx1x2(x1, x2); a21 = dx1x2(x1, x2); a22 = dx2x2(x1, x2);
                gr1 = -dx1(x1, x2); gr2 = -dx2(x1, x2);
                p1 = (a22 * gr1 - a12 * gr2) / (a11 * a22 - a12 * a21);
                p2 = (a11 * gr2 - a21 * gr1) / (a11 * a22 - a12 * a21);

                if (f(x1 + alpha * p1, x2 + alpha * p2) - f(x1, x2) <= eps * alpha * (p1 * dx1(x1, x2) + p2 * dx2(x1, x2)))
                {
                    k++;
                }
                else
                {
                    alpha *= beta;
                    continue;
                }
                x1 += alpha * p1;
                x2 += alpha * p2;
            }
            Console.WriteLine("x1=" + x1);
            Console.WriteLine("x2=" + x2);
            Console.WriteLine("f=" + f(x1, x2));
            Console.WriteLine("k=" + k);
            Console.ReadKey();
        }
    }
}

Проблема в том что программа ничего в консоль не выводит , помогите плиз найти ошибку. Спасибо