Разложение косинуса по формуле

[slow]

Помогите подправить программу с паскаля на с++

Задание:

Вычислить cos x по формуле разложением



с заданной точностью е.

вот что получилось, но считает не верно

#include <iostream.h>
#include <math.h>
using namespace std;
int main()
{
    setlocale(0,"Rus");
    double x,p,eps,cosin=1,n=0;
    cout<<"(eps>0): ";
    cin>>eps;
    cout<<"x: ";
    cin>>x;
    p=4*x*x/pow(M_PI,2);
        while(p>eps)
    {
        cosin=cosin*(1-p);
        n++;
        p=4*x*x/pow((2*n-1),2)*pow(M_PI,2);
    }
    cout<<cosin<<endl;
    cout<<cos(x)<<endl;
    system("PAUSE");
    return 1;
}

[Granus]

Начальное значение для n должно быть 1, а не 0

p=4*x*x/pow((2*n-1),2)*pow(M_PI,2);

- пи получилось в числителе. Нужно так

p=4*x*x/pow((2*n-1),2)/pow(M_PI,2);