On the numerical solution of the equation $$\frac{{\partial ^2 z}}{{\partial x^2 }}\frac{{\partial ^2 z}}{{\partial y^2 }} - \left( {\frac{{\partial ^2 z}}{{\partial x\partial y}}} \right)^2 = f$$ and its discretizations, I