#include<iostream>
using namespace std;
#include<cstdlib>
#include<cmath>
#include<math.h>
#include<iostream>
using namespace std;
#include<cstdlib>
#include<cmath>
#include<math.h>
/*
double numer(double n,double x)
{
return pow(-1,n)*pow(x,xponent);
}
double denom(double n, double(*f)(double))
{return f(n)*(3*n+1);}
*/
double factorial(double n)
{
int prod =1;
for(int b=n; n>0; --n)
{prod*=n;}
cout << prod;
}
int xponent(int n)
{
return 2*n+1;
}
int main()
{
long NMAX=1000;
double sum=0;
for(int n=0;n<1000; n++)
{
sum
+= (pow(-1,n
)*pow(0.3,xponent
(n
))/(factorial
(n
)*(3*n
+1);}
cout
<< (2/sqrt(M_PI
))*sum
;}
I2luY2x1ZGU8aW9zdHJlYW0+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7CiNpbmNsdWRlPGNzdGRsaWI+CiNpbmNsdWRlPGNtYXRoPgojaW5jbHVkZTxtYXRoLmg+CnVzaW5nIHN0ZDo6c2luOwp1c2luZyBzdGQ6OmNvczsKI2luY2x1ZGU8aW9zdHJlYW0+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7CiNpbmNsdWRlPGNzdGRsaWI+CiNpbmNsdWRlPGNtYXRoPgojaW5jbHVkZTxtYXRoLmg+CnVzaW5nIHN0ZDo6c2luOwp1c2luZyBzdGQ6OmNvczsKCi8qCmRvdWJsZSBudW1lcihkb3VibGUgbixkb3VibGUgeCkKewogCiAgcmV0dXJuIHBvdygtMSxuKSpwb3coeCx4cG9uZW50KTsKfQoKZG91YmxlIGRlbm9tKGRvdWJsZSBuLCBkb3VibGUoKmYpKGRvdWJsZSkpCntyZXR1cm4gZihuKSooMypuKzEpO30gICAKCgoqLwoKZG91YmxlIGZhY3RvcmlhbChkb3VibGUgbikKCQkJIHsKICBpbnQgcHJvZCA9MTsKICBmb3IoaW50IGI9bjsgbj4wOyAtLW4pCiAgICB7cHJvZCo9bjt9CiAgY291dCA8PCBwcm9kOwoJCQkgfQoKaW50IHhwb25lbnQoaW50IG4pCnsKICByZXR1cm4gMipuKzE7Cn0KaW50IG1haW4oKQp7CiBsb25nICBOTUFYPTEwMDA7CiAgZG91YmxlIHN1bT0wOwogIGZvcihpbnQgbj0wO248MTAwMDsgbisrKQogICAgewogCiAgICAgIHN1bSArPSAocG93KC0xLG4pKnBvdygwLjMseHBvbmVudChuKSkvKGZhY3RvcmlhbChuKSooMypuKzEpOwp9CiAgY291dCA8PCAoMi9zcXJ0KE1fUEkpKSpzdW07Cn0K
U3ltcyB4IHkKZj0xCmF4PTAKYng9MgpheT0wCmJ5PShzcXJ0KDQteF4yKSleMy8yCmFyZWE9NCZ0aW1lczsoKGludChpbnQoZix5LGF5LGJ5KSx4LGF4LGJ4KSkp
Syms x y
f=1
ax=0
bx=2
ay=0
by=(sqrt(4-x^2))^3/2
area=4×((int(int(f,y,ay,by),x,ax,bx)))