function [vals vecs] = matlabproject1 (A, c, iter)

vecs    = zeros(size(A));
n       = max(size(A));

for k=1:n
    intermat = inv(A - (c(k)*eye(n)));
    vecs(:, k) = rand(n,1);
 
    for j =1:iter
        vecs(:,k) = intermat*vecs(:,k);
        vecs(:,k) = vecs(:,k)/max(abs(vecs(:,k)));
        
    end
   
    test = A*vecs(:,k);
    place = find(test == max(test));
    vals(k) = test(place)/vecs(place,k);
  

end