This notebook and phase screens constitute a numerical experiment to calculate the error of the wavefront approximation using first N modes of the a) Zernike and b) Lukosz-Braat polynomials, c) SVD modes obtained with respect to the gradient-dot product, and the d) eigenfunctions of the Laplace operator with the Neumann boundary conditions. It's a complementary material to a paper submitted to Optics Express.