function y = FourIntegrand(x,j,cl,f,h,m)
% x variable d'espace dans [0,1]
% j rang du coeff de Fourier
% cl condition aux limites en 0 et 1 ˆ choisir entre
%     'dd' pour Dirichlet
%     'nd' pour Neumann-Dirichlet
%     'nn' pour Neumann
% f type de la fonction que l'on decompose
% h parametre (demi-largeur du support)
% m milieu du support
% C'est la fonction que l'on va integrer pour calculer le coeff de Fourier.
% Cette fonction sera integree entre m-h et m+h
%    'a' est la fonction indicatrice de l'intervalle [m-h, m+h] (paire)
%    'b' est la fonction chapeau sur [m-h, m+h] (paire)
%    'c' est une periode de la fonction sin sur [m-h, m+h] (impaire)
%    '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)
if f == 'a'
	y1 = 1 ;
elseif f == 'b'
    y1 = max (0*x , 1 - abs( (x-m)/h ) ) ;
elseif f == 'c'
    y1 = sin( pi*(x-m)/h + pi )  ;
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
if cl == 'dd'
	y = y1 .* sqrt(2) .* sin(pi*j*x) ;
elseif cl == 'nd'
	y = y1 .* sqrt(2) .* cos(pi*(j-1/2)*x) ;
elseif cl == 'nn'
    if j == 0
		y = y1 ;
	else
		y = y1 .* sqrt(2) .* cos(pi*j*x) ;
	end
end
