function Sp = spec ;
% Calcule le spectre de l'operateur div alpha^2 grad  
% avec alpha = a sur [0,b] \cup [1-b,1] et alpha = 1 sur [b, 1-b]
% pour les conditions de Dirichlet et la parite paire
% Appelle la fonction specfun.m
% On doit y inserer les valeurs de a et b
% On peut etre amene a regler l'Žchantillonage en frequence (s) pour 
% la detection du spectre.
% La mŽthode consiste en une dŽtection grossi�re de la localisation des
% zŽros d'une fonction, suivie de l'utilisation de la routine fzero.m

% Declaration des parametres globaux
global a
global b

% Initialisation des param�tres a et b
a = 0.5 
b = 0.33 ;

% Reglage de l'Žchantillonage en s pour la localisation grossi�re des zŽros
s = [0.1:0.1:100] ;

y = specfun(s) ;
z = sign(y) ;
plot(s,y,s,0*s,'c',s,z,'r')

k = 0 ;
for i = 2:length(s)
	if z(i-1)*z(i) == -1
		k = k + 1 ;
		Sp(k) = fzero('specfun',s(i)) ;
	end
	if z(i-1) == 0
		k = k + 1 ;
		Sp(k) = s(i-1) ;
	end
end
