 Solution

Let the source remain at $(0,0,-p)$ and let the observation point be $(X,Y,q)$. The paraxial approximation to the source-to-aperture-to-observer optical path length isCompleting the squares givesThe point $(x_0,y_0)$ is where the straight source-observer line intersects the aperture. Define the four dimensionless edge positionsOn the Cornu spiral $F(u)=C(u)+iS(u)$, draw the directed chord from $F(u_{x,-})$ to $F(u_{x,+})$, and a second directed chord from $F(u_{y,-})$ to $F(u_{y,+})$. Denote their complex differences by $D_x$ and $D_y$. The complex amplitude is proportional to their product: multiply their lengths and add their arguments, including the common propagation phase factor.

More explicitly, normalize by the unobscured field ${\psi }_0(X,Y)$ at this same observation point. The constant wave phase from completing the square is then already included in ${\psi }_0$, andEach directed chord is a one-dimensional Fresnel diffraction factor. One must multiply the two factors, rather than add the chords or use just one chord length for the square aperture. Setting $X=Y=0$ recovers $D_x=D_y=2F(w)$ and the on-axis answer.