The task is to derive the intensity pattern of a double slit diffraction, where light travels along the $z$-axis, the point of interest lies on the $x$-$z$ plane, the slits are wide (in $z$-direction) and infinite in extent in the $y$-direction. But how can this be derived by Fraunhofer diffraction, when one of the assumptions of Fraunhofer is, that the aperture is much smaller in size than the distance of aperture-POI? Because when using the Fraunhofer formula:
$$U(P)\propto e^{ik(s_0+d_0)} \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} t(x,y) e^{-ik(Xx+Yy)}\mathrm dx\mathrm dy$$
where $d_0, s_0$ are the distances of the source (here infinity) to the aperture and the distance of aperture - POI, and $X = x_s/s_0 + x_p/d_0$ and $Y = y_s/s_0 + y_p/d_0$. The problem comes from $y_p = 0, y_s=0 \implies Y=0$, so an integration over $y$ from $-\infty$ to $+\infty$ would diverge. How can I deal with this?