function [y,nu,iter]=powmeth(A,nmax,toll);
[n,n]=size(A);
y=rand(n,1);
y=y/norm(y);
nuold=y'*A*y
   y=A*y;
   y=y/norm(y);
   nu=y'*A*y;
   err=abs(nu-nuold);
iter=1;   
while err>toll & iter<nmax
   iter=iter+1;
   y=A*y;
   y=y/norm(y);
   nuold=nu;
   nu=y'*A*y;
   err=abs(nu-nuold);
end;


