Numerical integration based on simpsons rule implemented in Fortran90