getPsi()

function [Psi, epsilon]=getPsi(W)

Input:

• W: $\prod S \times N_s$ matrix of non-orthonormal functions minimizing E

Output:

• Psi: eigensolutions
• epsilon: eigenvalues

Provide a function of the above prototype in the file ``getPsi.m'' which computes the final solutions from the non-orthonormal $W$ according to the formulas Eqs. (6,7).

Hint: Given the matrix ``mu=Y'*H(Y)'', the code fragment

  [D, epsilon]=eig(mu); epsilon=real(diag(epsilon));

produces the matrix $D$ and the vector $\vec \epsilon$ in Eq. (9). (Again, we take a real part because rounding errors sometimes lead to tiny imaginary parts.)