Writing python code to sovle the folowing question.

Use Monte Carlo simulation with $N$ samples to compute

${\displaystyle \underset{x=0}{\overset{x=1}{}}}{\displaystyle \underset{y=0}{\overset{y=2}{}}}{e}^x{y}^{x}dxdy$

The number $N$ is chosen such that the estimator has an error equaling to $+-0\mathrm{.}05$ that corresponds to the $95\%$ confidence interval.