function AK = svdapprox(A,k)
tol=0.001
r=rank(A,tol)
if k>r           % return error message if k exceeds rank of matrix
        fprintf('Value of k must not exceed that of r' );
end
[U, S, V]=svd(A)             %calculate the svd of a