File:Onde cisaillement impulsion 1d 30.gif

 clear;
 
 Nx=35; // corde
 Ny=27;
 // Ny=5; // test
 Nxaff = 30; // corde
 Nyaff = 26;
 // Nyaff=4; // test
 L=3; // largeur (loi empirique)
 
 points_initiaux_x = zeros(Nx,Ny); // initialisation
 points_initiaux_y = zeros(Nx,Ny);
 for i=1:Nx
    for j=1:Ny
       points_initiaux_x(i,j) = i;
       points_initiaux_y(i,j) = j;
    end
 end
 points_finaux_y = points_initiaux_y;
 
 // fonction
 
 deff('[e]=profil(r)','e=r*exp(-(r/L)^2)'); // profil empirique
 
 // programme principal
 
 xmin = (Nx-Nxaff)/2+1;
 xmax = Nx-xmin;
 ymin = (Ny-Nyaff)/2+1;
 ymax = Ny-ymin;
 unite=0;dizaine=0; // pour le nom de fichier
 for r=-1:0.5:Nx+1 // position du front d'onde
    intensite = 2; // variation avec la distance
    for i=1:Nx
       for j=1:Ny
          // e = deplacement du point
          e = intensite*profil(i-r);
          points_finaux_y(i,j) = points_initiaux_y(i,j) + e;
       end
    end
    // dessin
    clf;
    plot2d(0,0,-1,"010"," ",[xmin,ymin,xmax,ymax]) 
    for i=1:Nx
       for j=1:Ny
          xpoly(points_initiaux_x(i,:),points_finaux_y(i,:));
          xpoly(points_initiaux_x(:,j),points_finaux_y(:,j));
       end
    end
    nom='test\test'+code2str(dizaine)+code2str(unite)+'.gif';
    xs2gif(0,nom,0);
    unite = unite+1;
    if unite == 10 then
       dizaine = dizaine + 1;
       unite = 0 ;
    end
 end