//////////////////////////////
// un schéma  qui marche (avec terme diffusion)
// équation de transport non linéaire
//////////////////////////////
function y=u0(x)
    y=zeros(x);
    y=exp(-10*(x-1/2)^2)
//    end    
endfunction
////////////////////
function y=f(x)
   y=c*x^2
endfunction
////////////////////
c=1/2;
N=400;
h=1/(N);
x=0:h:1;
U=zeros(1,N+3);
for i=2:N+2
    U(i)=u0(x(i-1))
end
///////////////////
Dt=1/401;Nt=600;
////
// schéma diffusif figure(1)
///////////////////
for i=1:Nt
  V=U
  drawlater
  clf() 
  plot2d(x,V(2:N+2),rect=[0,0,1.1,1.1]); 
  for j=2:N+2
    U(j)=(V(j-1)+V(j+1))/2-Dt/(2*h)*(f(V(j+1))-f(V(j-1)))  
  end
  U=[U(N+2),U,U(2)]
  drawnow 
end









