Suppose we have the following equation 1: $$\tag{1} A_G(x,y,z) = \frac{A_1}{q(z)} e^{-ik \frac{x^2 + y^2}{2q(z)}} $$ where $$ q(z) = z+iz_0 $$ and $i$ is equal to $\sqrt{-1}$.
Suppose we have another equation 2 (where $X(.)$, $Y(.)$, and $Z(z)$ are real functions): $$\tag{2} A(x,y,z) = X( \sqrt{2} \frac{x}{W(z)})Y(\sqrt{2} \frac{y}{W(z)})e^{iZ(z)}A_G(x,y,z) $$
Lastly, we have equation 3: $$ \Delta _T A(x,y,z) - 2ik \frac{ \partial A}{ \partial z} =0 \tag{3}$$
where $\Delta _T=\partial_{xx}+\partial_{yy}$ is the transverse Laplacian operator.
I need to show that substituting equation (2) into equation (3), given that equation (1) is also a solution of equation (3), will produce the following equation: $$ \frac{1}{X} ( \frac{\partial ^{2}X}{\partial u^{2}} - 2u \frac{\partial X}{\partial u}) + \frac{1}{Y} ( \frac{\partial ^{2}Y}{\partial v^{2}} - 2v \frac{\partial Y}{\partial v}) + kW^{2}(z) \frac{\partial Z}{\partial z}=0 $$ where $$ u=\frac{\sqrt2x}{W(z)} $$ and $$ v=\frac{\sqrt2y}{W(z)} $$ Can somebody give me some pointers as to how to begin here?