Вычислить значение функции арккосинуса, заданной с помощью ряда Тейлора

//http://upload.wikimedia.org/math/4/0/f/40f9e00ece5a364cb6c4ca8df9a125a1.png
#include <stdio.h>
#include <math.h>
int f2nm1(int n){ if (n>=1){ return (2*n-1)*f2nm1(n-1);} else { return 1;}}
int f2n(int n) { if (n>=1){ return (2*n)*f2n(n-1);}else { return 1;}}
double pw(double xn, int st){if (st>=1){return xn*pow(xn,st-1);}else {return (double)1.0;}}
double arcs(double cr,double xn, double e,int h)
{
double cr1=cr;
int y=2*h+1;
cr = cr - ((double)f2nm1(h)/(double)f2n(h))*(pw(xn,y)/(double)y);
if ((cr1 - cr) < e) { return cr;}
else { return  arcs(cr,xn,e,h+1);}
}
double arcos(double xn, double e)
{
double res = M_PI/2.0;
return (arcs(res,xn,e,0));
}
int main(int argc, char*argv[])
{
double xn, xk, e, dx;
xn=0.0, xk=1.0, e=0.001, dx=0.01;
while(xn+dx<=xk)
{
printf("arccos(%f)=%f\n",xn,arcos(xn,e));
printf("math: arccos(%f)=%f\n",xn,acos(xn));
xn+=dx;
}
return(0);
}