function [f] = flux(gam,r,u,e)
%
% calcul du flux des equations d'Euler
%
p = (gam-1.)*r.*e ; 
f(1,:) = r.*u ; 
f(2,:) = r.*u.^2+p ;
f(3,:) = (r.*(0.5*u.^2+e)+p).*u ;
   



