Please look carefully at the demo, as norm_u_0 = u.x.norm() is defined in the code.
Eigenvalue problems are solved with SLEPc.
This has been illustrated in many posts on the forum (for instance: Eigenvalue solver problem) and in the documentation I’ve already referred you to.
Please spend some time reading through the demos and related posts before asking further questions.