function [w] = soltr(w,nx,dx,dt,a,n) // // solution exacte de l'equation de transport // for i=1:nx xi = (i-1)*dx ; w(i) = sin((xi-a*n*dt)*2*pi) ; end endfunction