Вычисление функций с использованием их разложения в степенной ряд

double sn = 0;
            double se = 0;
            double y = 0;
            int i;
            for (double x = 0.1; x < 1; x = x + 0.9)
            {
                for (double n1 = 1; n1 < 10; n1 = n1 + 1)
                {
                    sn = ----------------------;
                }
                for (double n1 = 1; n1 < 10; n1 = n1 + 1)
                {
                    se = ------------------------;
                    se = se * 10000;
                    i = Convert.ToInt32(se);
                    se = (Convert.ToDouble(i) / 10000);
                }
                y = (Math.Pow(Math.E, x) + ((Math.Pow(Math.E, -x)) / (2)));
                Console.WriteLine("X = {0}; SN = {1}; SE = {2}; Y = {3}", x, sn, se, y);
                se = 0;
                sn = 0;
            }

    class Program
    {
        static void Main(string[] args)
        {
            double sn = 0;
            double se = 0;
            double y = 0;
            double s;
            int i;
            long f;
 
            for (double x = 0.1; x < 1; x = x + 0.1)
            {
                se = 1;
                sn = 1;
                for (int n1 = 1; n1 < 10; n1++)
                {                  
                    f  = factorial(2 * n1);
                    s  = Math.Pow(x, 2 * n1) / f;
                    sn = sn + s;
                    se = se + s;
                    se = se * 10000;
                    i  = Convert.ToInt32(se);
                    se = Convert.ToDouble(i) / 10000;
                }
                y = (Math.Pow(Math.E, x) + Math.Pow(Math.E, -x)) / 2;
                Console.WriteLine("X = {0}; SN = {1}; SE = {2}; Y = {3}", x, sn, se, y);
            }
            Console.ReadLine();
        }
        static long factorial(int n)
        {
            long factorial = 1;
            for (int i = 2; i <= n; i++) factorial *= i;
            return factorial;
        }
    }