I^2=\int_0^{\infty}\int_0^{\infty}\sin x^2\sin y^2 dx dy=\frac{1}{2}\int_0^{\infty}\int_0^{\infty}[\cos(x^2-y^2)-\cos(x^2+y^2)]dxdy=I_1-I_2