public static void MethodTeilor(double x0, double y0, double h, int n)
 {
     DFun f = new DFun(y0);
     int l = 0;
     for (int i = -2; i <= 2; ++i)
     {
         xx[l] = x0 + i * h;
         for (int j = 0; j <= 5; ++j)
         {
             yy[l] += f.Reval(j) / (factorial(j)) * Math.Pow(xx[l] - x0, j);
         }
         ++l;
     }
     Console.WriteLine(" X        Y        R");
     for (int i = 0; i <= 4; ++i)
     {
         if (xx[i] < 0) Console.Write("{0:f6}", xx[i]); else Console.Write(" {0:f6}", xx[i]);
         if (yy[i] < 0) Console.Write("{0:f6} ", yy[i]); else Console.Write(" {0:f6} ", yy[i]);
         Console.WriteLine("{0:f8}", f.Reval(6) / factorial(6) * Math.Pow(Math.Abs(xx[i] - x0), 6));
     }
 }
 public static void Ailer3(double x0, double y0, double h, int n)
 {
     double[] x = new double[2 * n];
     double[] y = new double[2 * n];
     x[0] = x0; y[0] = y0;
     for (int i = 1; i <= n; ++i)
     {
         x[i] = x[i - 1] + h;
         Fun f1 = new Fun(y[i - 1]);
         Fun f = new Fun(y[i - 1] + h * f1.ReVal());
         y[i] = y[i - 1] + h/2 * (f1.ReVal() + f.ReVal());
     }
     DFun g = new DFun(y0);
     Console.WriteLine(" X        Y        R");
     for (int i = 0; i <= n; ++i)
     {
         double ri = Math.Abs(r[i] - y[i]);
         if (x[i] < 0) Console.Write("{0:f6}", x[i]); else Console.Write(" {0:f6}", x[i]);
         if (y[i] < 0) Console.Write("{0:f6} ", y[i]); else Console.Write(" {0:f6} ", y[i]);
         Console.WriteLine("{0:f8}", ri);
     }
 }