function [i,mu,x]=puissinv(N)
%
%  Methode de la puissance ou methode de la puissance inverse.
%
%  N        : Dimension de la matrice
%
% construction de A
%
diag = 2*speye(N,N);
subd=-sparse(2:N,1:N-1,1,N,N);
A=diag+subd+subd';
%
% Vecteur initial
x=ones(N,1);
mu=dot(x,A*x)/dot(x,x);
for i=1:1000
  % A completer pour la methode de la puissance ou pour la methode de la 
  % puissance inverse
  x=A???x; 
  munew=dot(x,A*x)/dot(x,x);
  if (abs(munew-mu)/munew<0.0001)
    break
  end 
  mu=munew
end
x=x/norm(x); 
end

