function [power_rec, lvc] = powrec(A,x0,maxit,tol)
locnit=2;
ok=0;
power_rec = [0 Inf Inf NaN 0];
nit=0;
lvc=x0;
while (locnit == 2 & ~ok & nit < maxit ),
	if nit == 0 	
		[lvl,lvc,ler,locnit,pvl,pvc,per,ok]= power(A,lvc,2,tol);
	else
		[lvl,lvc,ler,locnit,pvl,pvc,per,ok]= power(A,lvc,2,tol,lvl,ler);
	end
	if locnit==0,
		if nit == 0 ,
			power_rec = [nit+locnit lvl ler 0 ok];
		else 
			power_rec = [ power_rec ; [nit+locnit lvl ler 0 ok]];
		end
	elseif  locnit==1,
		if nit > 0,
		power_rec = [power_rec ; [nit+locnit lvl ler ler/power_rec(nit+1,3) ok]];
		else
			power_rec = [power_rec ; [nit+locnit lvl ler ler/per ok]];
		end
	else 
		if nit > 0, 
		power_rec = [power_rec ; 
				[nit+1 pvl per per/power_rec(nit+1,3) 0]; 
				[nit+locnit lvl ler ler/per ok]];
		else
%		power_rec, nit, pvl, per,
		power_rec = [power_rec ; 
				[nit+1 pvl per per/power_rec(nit+1,3) 0]; 
				[nit+locnit lvl ler ler/per ok]];
		end
	end
	nit = nit+locnit;
end;
return

