How do I solve
$$\int _1^{\infty} e^{-x^2} dx?$$
I solved like this.
let $I = \int _1^{\infty} e^{-x^2} dx$
then $I^2 = \int _1^{\infty} \int _1^{\infty} e^{-x^2 -y^2} dxdy$
If I change the cartetian coordinate to the polar coordinate.
$I^2 = \int _0^{\frac{\pi}{2}} \int _1^{\infty} e^{-r^2} rdrd\theta$
$= \int _0^{\frac{\pi}{2}} [\frac{e^{-r^2}}{-2}] _1^{\infty} d\theta$
$= \int _0^{\frac{\pi}{2}} \frac{1}{2e} d\theta$
$= \frac{1}{2e} \frac{\pi}{2}$
and then,
$ I = \frac{1}{2} \sqrt\frac{\pi}{e}$
I got The solution is $\frac{1}{2}\sqrt{ \frac{\pi }{ e} }$
but in the integral calculator the answer is $0.13940279...$.(https://www.integral-calculator.com/)
What did I do wrong?
Nope. If you consider $I^2$ you get the integral of $e^{-(x^2+y^2)}$ over the region $x,y\geq 1$, not the region $x^2+y^2\geq 1$.
Your integral simply equals $$ \frac{\sqrt{\pi}}{2}-\int_{0}^{1}e^{-x^2}\,dx = \frac{\sqrt{\pi}}{2}-\sum_{n\geq 0}\frac{(-1)^n}{n!(2n+1)}.$$