#!/usr/bin/env python

from gaussxw import gaussxwab
from math import exp

def f(z):
    return exp(-(z/(1-z))**2)/(1-z)**2

N = 50
a = 0.0
b = 1.0
x,w = gaussxwab(N,a,b)

s = 0.0
for k in range(N):
    s += w[k]*f(x[k])
print(s)