Расчет суммы ряда csc(x) (с++)

[system_error]

Дана формула косеканса, разложение косеканса в ряд.
Нужно составить программу для расчета суммы ряда и расчета по формуле и нахождения разности между этими значениями. значения х изменяются от (-1,0)&(0,1). Необходимые формулы находятся в приложенном файле.
Я составила такой цикл:

double b[]={1./6, 1./30, 1./42, 1./30, 5./66, 691./2730, 7./6, 3617./510,
43867./798, 174611./330, 854513./138};//массив чисел Бернулли
for(k=0;k<11;k++)
  {
   t=1.;
   double z=(2*(step(2, a[k])-1))/fact(a[k]+1);
   s=step(x, a[k]);
   t=t+(z*(b[k]*s));
  }
//---------------------------------------------------------------------------
int fact(int n)//нахождение факториала
{
 int q=1;
 for(int i=1;i<=n;i++)
  q=q*i;
 return q;
}

double step(double w, int n)//возведение в степень
{
 double p=1.;
 for(int j=1;j<=n;j++)
 p=p*w;
return p;
}

но при расчете получается слишком большая погрешность. Помогите, пожалуйста, найти ошибку в формуле.
Заранее благодарю.

[Sweta]

1 Где Вы учитываете значение 1/х?
2 В Вашей формуле 1/(2n)!, а не на n!.
Вам необходимо определить коэффициент изменения как А+1/Аn и тогда определять сумму ряда.

[system_error]

еще есть массив int a[]={-1,1,3,5,7,9,11,13,15,17,19};
оттуда и получается 1/x