We employ the method from [1]} and [2]} to estimate the \(h_z\)  and its radial behavior. The basic idea is to fit the expression for the surface density \(\Sigma (Z)\)  derived from the Jeans equations to the measured velocity dispersions and covariance:

\(\Sigma (Z)&=&\frac{1}{2\pi G}\Big [k_1\cdot \int _0^Z\sigma _R^2\rm {d}Z+k_2\cdot \int _0^Z\sigma _\phi ^2\rm {d}Z+k_3\cdot \sigma _{RZ} \nonumber \\&&+\frac{\sigma _Z^2}{h_z}-\frac{\partial \sigma _Z^2}{\partial Z}+|Z|\sigma _{RZ}\frac{\partial }{\partial R}\Big (\frac{1}{h_z}\Big )+\int _0^Z|Z|\frac{\sigma _R^2}{R}\frac{\partial }{\partial R}\Big (\frac{1}{h_z}\Big )\rm {d}Z \nonumber \\&&+\int _0^Z|Z|\sigma _R^2\frac{\partial ^2}{\partial R^2}\Big (\frac{1}{h_z}\Big )\rm {d}Z\Big ] \nonumber \\&=&2\rho (R,\,Z=0)\cdot h_z(R)\cdot \big [1-e^{-Z/h_z(R)}\big ].\) 
