function y = eigvec(x,s)
% Calcul de pre-vecteurs propres normalises en norme L2.
% Pour s dans le spectre, on obtient les vecteurs propres.

% Declaration des parametres globaux
global a
global b

y = (x<=b) .* sin(pi*s*x/a) * cos(pi*s*(b-1/2)) +...
	(x>b & x<1-b) .* cos(pi*s*(x-1/2)) * sin(pi*s*b/a) +...
	(x>=1-b) .* sin(pi*s*(1-x)/a) * cos(pi*s*(b-1/2)) ;
	
normL2 =  cos(pi*s*(b-1/2))^2 * (b - sin(2*pi*s*b/a) /(2*pi*s/a)) +...
		  sin(pi*s*b/a)^2 * ( (1/2-b) - sin(pi*s*(2*b-1)) / (2*pi*s)) ;
normL2 = sqrt(normL2) ;

y = y / normL2 ;
