计算二重积分 $\iint_D \mathrm{e}^{x^2} \mathrm{~d} x \mathrm{~d} y$ ，其中 $D$ 是第一象限中由直线 $y=x$ 和 $y=x^3$ 围成的封闭区域。