function Ondes1DTrans(tinit,tfin,dt,K,f,h)
% tinit le temps initial o� l'on regarde 
% (les conditions initiales sont toujours en t = 0).
% tfin le temps final d'observation
% dt l'increment de temps
% K la frequence de coupure
% Les conditions aux limites sont Dirichlet
% Les conditions initiales sont paires par rapport au milieu 0.5
% Voir la doc pour les autres parametres dans FourIntegrandTrans.m
% Onde1DTrans montre la propagation des ondes dans un milieu anisotrope 1D
% (decrit dans FourIntegrandTrans.m)
% Appelle spec.m qui calcule le spectre et determine les param�tres a et b
% du materiau qui sont declares en global.
% Appelle eigvec.m qui calcule les vecteurs propres.
% Pour l'utiliser, faire d'abord tourner spec.m en definissant a et b
% ˆ l'intŽrieur, et vŽrifier que ses param�tres internes sont convenablement
% r�glŽs pour obtenir le nombre dŽsirŽ de valeurs propres.

% Declaration des parametres globaux
global a
global b

% Preparation des proprietes de la figure 1
figure(1)
WHITEBG(1,'k')
set(1,'Renderer','zbuffer')

% Calcul du spectre
Spectre = spec ;
K = min(K,length(Spectre)) 
% Calcul des coefficients de Fourier
m = 1/2 ;
for j = 1:K
	FC(j) = quad8('FourIntegrandTrans',max(0,m-h),min(1,m+h),[],[],...
	Spectre(j),f,h) ;
end
% Formule de representation
x = [0:0.001:1] ;
for tt = tinit:dt:tfin
	z = 0*x ;
	for j = 1:K
		z = z + FC(j) * eigvec(x,Spectre(j)) * cos(Spectre(j)*pi*tt) ;
	end
% 	Graphe
	zmin = min(z,0) ;
	zmax = max(z,0) ;
	for i=1:length(x)
		if zmin(i) == 0, zmin(i) = NaN ; end
		if zmax(i) == 0, zmax(i) = NaN ; end
	end	
	plot(x,zmin,'c',x,zmax,'b') ;
	title (['t  =  ',num2str(tt)]) ;
	axis([0 1 -1 1])
	pause(0.05)
end
