salut ,les liens ce sont utile mais malheureusement on n'a pas encore fait des cours concernant ceci et moi je cherche des trucs trés détaillé.voilà le programme en scilab que j'ai fait mais je n'ai pas pu le terminer :
lg = 10. ; // lg = demie longueur du domaine
dx = 0.05 ; // dx = pas d'espace
nx = lg/dx ; // nx = demi nombre de mailles
cfl = 0.4 ; // cfl
dt = dx*dx*cfl ; // dt = pas de temps
nt = 500 ; // nt = nombre de pas de temps effectues
//
// initialisation
//
x=zeros(1,2*nx+1) ;
u0=zeros(1,2*nx+1) ;
for i=1:2*nx+1
x(i) = (i-nx-1)*dx ;
u0(i) = max(0.,1.-x(i)**2) ;
end
u=u0 ;
up=u0 ;
um=u0 ;
uexacte=u0 ;
tics=[4,10,4,10];
plotframe([-lg,-0.1,lg,1.1],tics);
plot2d(x,u0,1,"000")
xtitle ('donnee initiale' ,' ',' ') ;
halt() ;
////////////////////////////////////////////////////////////////
// schema explicite: cfl=0.4
////////////////////////////////////////////////////////////////
//
for n=1:nt
//
up = sh('+1',u) ;
um = sh('-1',u) ;
u=u + dt/(dx*dx)*(up+um-2*u) ;
plotframe([-lg,-0.1,lg,1.1],tics);
plot2d(x,u,[1,1],"100","schema explicite")
plot2d(x,u0,[2,2],"100","donnee initiale")
xtitle('schema explicite, cfl=0.4',' ',' ');
xpause(1000) ;
end
avec sh une fct qui décale un tableau d'element de +1 ou -1 .
ici j'ai traité le cas de shéma explicite , est ce que vous pouvez m'aider à le terminer ?merci et by.


