integrate Exp[-r^2/a^2]