#include
using namespace std;
double power(double x,int n){
double val=1.0;
while(n--)
val*=x;
return val;
}
double arctan (double x){
int sum;
for (int i=1,z=2 ; i<=15 ; i+=2,z++){
if(z%2!=0)
{sum=((-(int)power(x,i))/i);
}
else
{sum=((int)power(x,i))/i;
}
sum+=sum;
}
return sum;
}
int main(){
double PI;
double a=16*(double)arctan(1/5.0);
double b=(double)4*arctan(1/239.0);
PI=a-b;
cout<<"PI= "<<PI<<endl;
return 0;
}