Calculate nen properly
Use higher quadrature, and gaussian density