  function Apr14()
% If F([lambda;v]) = [(M + lambda*N + lambda^2*P)*v;v'*v-1]
% then the partial of F with respect to lambda is 
%                    [(N+2*lambda*P)*v;0]
% and the partial F with respect to v(k) is 
%                    [(M + lambda*N + lambda^2*P)*I(:,k);2v(k)].
% The Jacobian is the assembly of these vectors so
%
%  J = [N+2*lambda*v (M+lambda*N+lambda^2*P);0 2*v'];