The correct option is D π2
Given integral is
I=∫∞0∫∞0e−x2.e−y2dxdy
=∫∞0∫∞0e−(x2+y2)dxdy
Put x=rcosθ,y=rsinθ,|J|=r (Jacobian of polar coordinate system).
Then to cover whole 1st Quadrant θ varies from 0 to π2 and r varies from 0 to ∞
So, I=∫∞0∫∞0e−(x2+y2).dxdy
=∫π/2θ=0∫∞r=0e−r2.(rdrdθ)
=π2∫πr=0r.e−r2.dr=π4 (Put r2=t)