Change to Gauss-Kronrod algorithm

Change to G-K algorithm. So we use quadgk().