//////////////////////////////
// un schéma centré qui ne marche pas
// équation de transport linéaire
//////////////////////////////
function y=u0(x)
    y=zeros(x);
    if (x<0) then
       y=1 
    end    
endfunction
////////////////////
function y=f(x)
   y=c*x
endfunction
////////////////////
c=1;
N=100;
h=10/N;
x=-5:h:5;
U=zeros(1,N+1);
for i=1:N+1
    U(i)=u0(x(i))
end
///////////////////
Dt=1/100;Nt=500;
////////////////////
// schéma centré
///////////////////
figure(0)
for i=1:Nt
  V=U
  drawlater
  clf()
  plot2d(x,[0,U(2:N),0],rect=[-5,0,5,4]);  
  for j=2:N-1
    U(j)=V(j)-Dt/(2*h)*(f(V(j+1))-f(V(j-1)))  
  end
  drawnow 
end  
//////////////////










