cipher

 25 June 12:59   % Aqueduct with amplification - endure adapted 11/03/99

    % Reproduit avec laimable permission de Prof. L. Mongeau

    clear all;

    clc;

    close all;

    Nsteps=30;

    f=86

    rho=1.15;

    c=344.;

    rhoc=rho

    l=1; %length of duct

    a1=0.0762/2;

    S1=pi

    a2=0.1524/2;

    S2=pi

    m=S2/S1;

    w=2

    k=w/c;

    t=linspace(2

    alpha1=2.93e-5/a1.

    alpha2=2.93e-5/a2.

    khat=k-i

    kl=khat

    U0=0.01;

    R=(1-m)/(m+1);

    A=rho

    B=R

    C=A+B;

    x=linspace(0,2,100); % computational domain

    y=linspace(-0.2,0.2,32);

    for m=1:50

    for k=8:24

    pt(k,m)=A

    end

    end

    for m=51:100

    for k=1:32

    pt(k,m)=C

    end

    end

    pmag=abs(pt);

    

    pphase=angle(pt);

    clim=max(max(pmag));

    V=[-clim clim]; % sets blush scheme

    pt=pmag.

    % blue: minimum abrogating pressure

    figure(1)

    pt2=real(pt);

    H=pcolor(x,y,pt2);

    axis according % to annihilate distortion

    axis([0 2 -.2 .2])

    M=moviein(Nsteps); % for Matlab cine animation

    for k=1:Nsteps

    pt2=real(pt

    H=pcolor(x,y,pt2);

    axis equal

    xlabel(x);

    ylabel(y);

    title(pressure)

    shading interp;

    caxis(V);

    axis([0 2 -.2 .2])

    M(:,k)=getframe;

    % to abundance the images for web announcement (bitmaps)

    % filnme=strcat(spher,int2str(k));

    % eval([print -dbitmap filnme -f]);

    end

    movie(M,20)