the program now fully works, including eigenfunctions computation