Actually what I want to prove is $|\int_{D} e^{z^2+1} dz| \le 2\pi e^2$ for $D$ = the unit circle (centred at origin).
What I know is $|\int_\gamma f(z)dz| \le AB$ for $A$ the length of $\gamma$ (for whatever $\gamma$ you need eg smooth, continuously differentiable, simple/injective, Jordan, whatever) and $B \in \mathbb R$ assuming $|f(z)| \le B$ for all $z \in image(\gamma)$.
Then obviously $2\pi$ is the length. As for the $e^2$, I think of 2 ways.
Way 1: We have $|z^2+1| \le |z|^2+1 \le 1+1=2$ and then somehow $|e^{z^2+1}| \le e^{|z^2+1|} \le e^2$.
I believe in general for any $f: G \to \mathbb C$ for any subset $G$ of $\mathbb C$, we have that $|e^{f(z)}| \le e^{|f(z)|}$ because for $f=u+iv$, we just do $|e^{f(z)}|=e^{u(z)}$ and $e^{|f(z)|}=e^{\sqrt{u^2(z)+v^2(z)}}$.
Question 1: Is this correct?
Way 2: $$e^{z^2+1} = e^{x^2-y^2+1}e^{i(2xy)}$$
$$=e^{x^2-y^2+1} [\cos(2xy) + i\sin(2xy)]$$
Then
$$|e^{x^2-y^2+1} [\cos(2xy) + i\sin(2xy)]|$$
$$= |e^{x^2-y^2+1}| |\cos(2xy) + i\sin(2xy)|$$
$$= e^{x^2-y^2+1} |\cos(2xy) + i\sin(2xy)|$$
$$= e^{x^2-y^2+1} $$
$$= e^{x^2-y^2+1} \le e^2$$
because $x^2-y^2+1 \le 2$ because $x^2-y^2+1=x^2-(1-x^2)+1=2x^2 \le 2$ because $x^2 \le 1$ because $x \le 1$
Question 2: Is this correct?
What you wrote is correct. $$ |e^w| = e^{\operatorname{Re}(w)} \le e^{|w|} $$ holds for all complex numbers $w$, so $$ |e^{z^2+1}| \le e^{|z^2+1|} \le e^{|z|^2+1} = e^2 $$ for $z = 1$.
However, this is not needed to estimate the integral. According to Cauchy's integral theorem, $\int_{D} e^{z^2+1} dz$ is zero.