Clear[inner, x, y];
inner[x_] = 
        Integrate[
          E^(-x^2 - y^2), {y, 0, Sqrt[16 - x^2]}];
Integrate[inner[x], {x, 0, 4}]