#include<stdio.h>
#include<math.h>
#define N 20
float f(float x)
{
return(x*x);
}
int main()
{
char ch1[80];
double a,d,m,n,s,tru,e;
tru=1.0/3.0;
for(m=2; m<N; m++){
d=1.0/m;
s=0;
for(a=0; a<m; a++){
s+=(f(a*d)+f((a+1.0)*d))*d/2.0;
}
}
}
I2luY2x1ZGU8c3RkaW8uaD4KI2luY2x1ZGU8bWF0aC5oPgojZGVmaW5lIE4gMjAKCmZsb2F0IGYoZmxvYXQgeCkKewogIHJldHVybih4KngpOwp9CgppbnQgbWFpbigpCnsKICAgIGNoYXIgY2gxWzgwXTsKICAgIGRvdWJsZSBhLGQsbSxuLHMsdHJ1LGU7CiAgICB0cnU9MS4wLzMuMDsKZm9yKG09MjsgbTxOOyBtKyspewogICAgICAgIGQ9MS4wL207CiAgICAgICAgcz0wOwogICAgICAgIGZvcihhPTA7IGE8bTsgYSsrKXsKICAgICAgICAgICAgcys9KGYoYSpkKStmKChhKzEuMCkqZCkpKmQvMi4wOwogICAgICAgIH0KICAgICAgICBhPXNxcnQoMi4wKTsKICAgICAgICBlPXNxcnQoKHMtdHJ1KSoocy10cnUpKTsKICAgICAgICBwcmludGYoIiVmXHQlZlxuIixtLGUpOwogICAgfQp9
IDQqaW50ZWdyYXRlKHNxcnQoMip4LXgqeCktMS8yLCB4LDEtc3FydCgzLzQpLDEvMik7IAo=
4*integrate(sqrt(2*x-x*x)-1/2, x,1-sqrt(3/4),1/2);