function y = FourIntegrandTrans(x,s,f,h)
% x variable d'espace dans [0,1]
% s parametre spectral
% C'est Dirichlet pour des fonctions paires pour div A^2 grad
% avec A = a sur [0,b] \cup [1-b,1] et A = 1 sur [b, 1-b].
% f type de la fonction que l'on decompose
% h parametre (demi-largeur du support)
% m milieu du support
% 'a' est la fonction indicatrice de l'intervalle [m-h, m+h] (paire)
% 'b' est la fonction chapeau sur [m-h, m+h] (paire)
% 'd' est le carre de la precedente (paire)
% 'e' est une demi-periode de la fonction cos sur [m-h, m+h] (paire)
% 'f' est le carre de la precedente (paire)
% FourIntegrandTrans.m
% est la fonction que l'on va integrer pour calculer le coeff de Fourier.
% Cette fonction sera integree entre m-h et m+h
% Elle appelle la fonction eigvec.m qui calcule les vecteurs propres.

% Declaration des parametres globaux
global a
global b

m = 1/2 ;

if f == 'a'
	y1 = 1 ;
elseif f == 'b'
    y1 = max (0*x , 1 - abs( (x-m)/h ) ) ;
elseif f == 'd'
    y1 = sin( pi*(x-m)/h + pi ).^2  ;
elseif f == 'e'
    y1 = cos( pi*(x-m)/(2*h) )  ;
elseif f == 'f'
    y1 = cos( pi*(x-m)/(2*h) ).^2  ;
end
y = y1 .* eigvec(x,s) ;
