вычисление cos разложением в ряд

[bondik]

Вообщемто вот сам ряд

[bondik]

вот код:

double cosin(double x,double pogr){
double k=-1;double chlen;double n=1;
double xn=x*x;double summ=0;
    do {                     //цикл разложения в ряд
    chlen=(k*xn)/(fact(2*n));
    n++;k=-k;xn*=x*x;
    summ+=chlen;

	 } while (chlen>pogr); //пока член больше погрешности
return summ;
}

double fact(double n)
{double fac=1;
for (int i=1;i<=n;i++)
fac*=i;
return fac;

}

трассировкой выяснил что всегда вычслияет только два члена и выходит из цикла,независимо от погрешности =(

Подскажите пожайлуйста в чем ошибка

[bondik]

Сам нашел ошибку...
Надо в условие цикла член по модулю брать =(

[bondik]

Может кому пригодится разложения : cos,arcctg,arth,sech

//вычисление cos
double cosin(double x,double pogr){
double k=-1;double chlen=0;double n=1;
double xn=x*x;double summ=0;
do {                                //цикл разложения в ряд
	summ+=chlen;
    chlen=(k*xn)/fact(2*n);
    n++;k=-k;xn*=x*x;
    } while (fabs(chlen)>pogr); //пока член больше погрешности
return summ;
}


//вычисление arcctg
double arccatg(double x,double pogr){
double k=-1;double chlen=0;double n=1;
double xn=x*x*x;double summ=0;
do {                                //цикл разложения в ряд
	summ+=chlen;
    chlen=(k*xn)/((2*n)+1);
    n++;k=-k;xn*=x*x;
    } while (fabs(chlen)>pogr); //пока член больше погрешности
return 1.57-summ+x;
}


//вычисление arth
double arthh(double x,double pogr){
double chlen=0;double n=1;
double xn=x*x*x;double summ=0;
do {                                //цикл разложения в ряд
	summ+=chlen;
    chlen=xn/((2*n)+1);
    n++;xn*=x*x;
    } while (fabs(chlen)>pogr); //пока член больше погрешности
return summ+x;
}


//вычисление sech
double sechh(double x,double pogr){
double chlen=0;double n=1;
double xn=x*x;double summ=0;
do {                                //цикл разложения в ряд
	summ+=chlen;
    chlen=xn/fact(2*n);
    n++;xn*=x*x;
    } while (fabs(chlen)>pogr); //пока член больше погрешности
return 1/summ;
}



double fact(double n)
{double fac=1;
for (int i=1;i<=n;i++)
fac*=i;
return fac;

}

[vitaliz]

При вычислении косинуса
summ=0-неправильно.
правильно s=1; ведь ряд начинается с 1.