fixed small bug
calculated integral exp(-x^2) * x^k