Hello, I try to solve this problem, and I think a) wasn't too hard, I have the following solution:
##H = \lambda (\frac{\vec{S^2-(\vec{S_1}^2+\vec{S_2}^2)}{2})##.
I struggle with 2. I find it very abstract. When I have H as a matrix I know how to calculate eigenvalues, but I don't know how...