How Gauss solved the integral of e^(-x^2) from 0 to infinity