If we know the exact solution and approximate solution of a PDE then we can compute the error of both solutions in the discrete norm as given: $\|u-u_N\|^2_M$ = $\frac{2\pi}{M}\sum _{j=0}^{M-1}(u(z_j, t)-u_N(z_j, t))^2$. I am unable to understand this relation... kindly explain.
Reference: it is written in the book Implementing Spectral Methods for Partial Differential Equations by David A. Kopriva