function [X,~] = my564proj1(A,n,k,X0,pmu,...) % % INPUT: % A - n by n real symmetric matrix % n - size of the matrix % k - number of eigenpairs to compute % X0 - n by k matrix as a starting point % pmu - parameter mu > 0 (optional) % OUTPUT: % X - computed approximate solution