As I’ve already told you, put your code between three single quote `, otherwise it is not copy/paste-able. Make sure it’s runnable, it’s clearly not the case here, e.g., what’s sigma?
code should appear like this;
It’s wrong to put the shift in the left-hand side matrix. You have to remove - sigma u1 * u2 from the varf.
How is this post different than this one? Also, this question is not related to PDE solvers, but eigensolvers, if I’m not mistaken? What are the eigenvalues you are looking for? Closest to the shift? What are the “analytical results” you are mentioning, i.e., exact values you are supposed to get?