calling x,w = gaulegf(1.0, 2.0, 4)
x= [0, 1.0694318442029738, 1.3300094782075718,
       1.6699905217924282, 1.9305681557970262]
w= [0, 0.17392742256872379, 0.32607257743127321,
       0.32607257743127321, 0.17392742256872379]
area= 2.33333333333