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

#include <stdio.h>
#include <math.h>
#include <locale.h>
float riad(float, float, float, int&);
int main()
{int vid, n, i ;
 float a, x, xn, xk, h;
 float eps, rez;
 setlocale(LC_ALL,"Russian") ;
 for(;;)
 {
 
  printf("\n                      Vid deistviya:\n");
  printf("                  1 - poluxhenie tablicy znachenii\n");
  printf("                  2 - zavershenie programmy\n");
  printf("                  Vybirite vid deystviya ->");
  scanf("%d",&vid);
  if (vid == 1)
  {
 
   printf("\n Vvedite znachenie osnovaniya ") ;
   scanf("%f", &a) ;
   printf(" Vvedite nachalnoe znachenie stepeni ");
   scanf("%f", &xn);
   printf(" Vvedite konechnoe znachenie stepeni ");
   scanf("%f", &xk);
   if (xk <= xn)
      {printf("\n Oshibka! Znachenie xk d.b. > xn");
       return 1;
      }
   printf(" Vvedite shag izmeneniya stepeni ");
   scanf("%f", &h);
   if (h <= 0)
      {printf("\n Oshibka! Znachenie d.b. >0");
       return 1;
      }
   printf(" Vvedite tochnost vychisleniya ");
   scanf("%f", &eps);
   if ((eps <= 0) || (eps > 0.1))
      {printf("\n Oshibka! Znachenie d.b. >0 i <=0.1");
       return 1 ;
      }
 
   printf("\nTablica znachenii %4.1f v stepeni x:\n", a) ;
   for(i = 0 ; i < 50 ; i++) printf("-") ;
   printf("\n|Stepen |  Summa |Kol.slag.|Kontrolnoe znachenie|\n") ;
   for(i = 0 ; i < 50 ; i++) printf("-") ;
   x = xn ;
   do{n = 0 ;
      rez = riad(x, a, eps, n) ;
      printf("\n|%6.2f  |%8.4f|%6i   |%14.4f      |\n", x, rez, n, 2*(pow((x-1), (2*n+1)))/((2*n+1)*(pow((x+1),(2*n+1)))));
      for(i = 0 ; i < 50 ; i++) printf("-") ;
      x += h ;
     }while(x <= xk) ;
   }
   else break ;
 }
 return 0 ;
}
 
float riad(float x, float a , float eps, int& n)
{float rez = 0 ;
 float pr = 1 ;
 n = 0 ;
 while (fabs(pr) > eps)
       {rez += pr ;
    n++ ;
    pr *= (pow((x-1), (2*n+1)))/((2*n+1)*(pow((x+1),(2*n+1))));
       }
    rez=2*rez;
 return rez ;
}