#include <stdio.h>
int factorial(int n) {
	int f = 1;
	for (int i=1; i<=n; i++)
		f *= i;
	return f;
}
int main(){
    int n;
    long double result=0;
    for(n=0; n<=5; n++){
        if(n%2==0){
            result+=(((long double) 1/(long double)(factorial(n)*(2*n+1))));
        } else {
            result-=(((long double) 1/(long double)(factorial(n)*(2*n+1))));
        }
    }        
    printf("The Gaussian integral from 0 to 1 is %Lf\n", result);
}
 

I2luY2x1ZGUgPHN0ZGlvLmg+CmludCBmYWN0b3JpYWwoaW50IG4pIHsKCWludCBmID0gMTsKCWZvciAoaW50IGk9MTsgaTw9bjsgaSsrKQoJCWYgKj0gaTsKCXJldHVybiBmOwp9CmludCBtYWluKCl7CiAgICBpbnQgbjsKICAgIGxvbmcgZG91YmxlIHJlc3VsdD0wOwogICAgZm9yKG49MDsgbjw9NTsgbisrKXsKICAgICAgICBpZihuJTI9PTApewogICAgICAgICAgICByZXN1bHQrPSgoKGxvbmcgZG91YmxlKSAxLyhsb25nIGRvdWJsZSkoZmFjdG9yaWFsKG4pKigyKm4rMSkpKSk7CiAgICAgICAgfSBlbHNlIHsKICAgICAgICAgICAgcmVzdWx0LT0oKChsb25nIGRvdWJsZSkgMS8obG9uZyBkb3VibGUpKGZhY3RvcmlhbChuKSooMipuKzEpKSkpOwogICAgICAgIH0KICAgIH0gICAgICAgIAogICAgcHJpbnRmKCJUaGUgR2F1c3NpYW4gaW50ZWdyYWwgZnJvbSAwIHRvIDEgaXMgJUxmXG4iLCByZXN1bHQpOwp9Cg==