Решение системы нелинейных уравнений методом Ньютона

   class Program
    {
        public delegate double Function(double x);
        public static double Newton(Function f, Function fprime, double x0, double epsilon)
        {
            double f0 = f(x0);
            double x = x0;
            while (Math.Abs(f(x)) > epsilon)
            {
                x -= f0 / fprime(x);  // fprime(x) - производная f(x)
                f0 = f(x);
            }
            return x;
        }
        static double F1(double x)
        {
//            return Math.Exp(x)-x * x;
            return x * x * x * x - 2 * x * x * x + x - 132;
        }
        static double F1prime(double x)  // производная F1
        {
//            return Math.Exp(x) -2 * x;
            return 4* x * x * x  - 6 * x * x + 1;
        }
        static void Main(string[] args)
        {
            try
            {
                double x = Newton(F1, F1prime, -1.0, 0.0001);
                Console.WriteLine("Solution from the Newton method: " + x.ToString());
                Console.WriteLine("Solution confirmation: f(x) = " + F1(x).ToString());
                x = Newton(F1, F1prime, 0, 0.0001);
                Console.WriteLine("Next solution :                  " + x.ToString());
                Console.WriteLine("Solution confirmation: f(x) = " + F1(x).ToString());
 
                Console.WriteLine();
 
                Console.ReadLine();
            }
            catch (Exception ex)
            {
                Console.WriteLine("Fatal error: " + ex.Message);
                Console.ReadLine();
            }
        }
    }