function TestEllipsePoint()
% Checks out EllipsePoint
clc
randn('seed',0);
A = randn(2,2);
[U,S,V] = svd(A);
fprintf('   || x ||    || Ax ||     Required \n')
fprintf('---------------------------------------------\n')
for lambda = linspace(0,1,11)
   x = EllipsePointCVL(A,lambda);
   normx = norm(x);
   normAx = norm(A*x);
   value = lambda*S(1,1) + (1-lambda)*S(2,2);
   fprintf('%10.6f   %10.6f    %10.6f   \n',normx,normAx,value)
end
fprintf('\n\nA= eye(2,2). Special Case when the singular values are repeated:\n\n')
A = eye(2,2);
x = EllipsePointCVL(A,.5);
normx = norm(x);
normAx = norm(A*x);
value = lambda;
fprintf('%10.6f   %10.6f    %10.6f   \n',normx,normAx,value)


