function diagstamov (zc)
% Dominic Grenier janvier 2016

Gammac= (zc-1)/(zc+1);

d=0:0.025:1; % longueur d'onde lambda=1
t=0:0.05:1; % période T=1
[D,T]=meshgrid(d,t);

Vd=exp(j*2*pi*D)+Gammac*exp(-j*2*pi*D);
Id=exp(j*2*pi*D)-Gammac*exp(-j*2*pi*D);

vdt=abs(Vd).*cos(2*pi*T+phase(Vd));
idt=abs(Id).*cos(2*pi*T+phase(Id));

disp('GEL-19881 présente');
pause(1);
disp('un film de Dominic Grenier');

h=figure(1);
clf;
M=moviein(length(t),h);
for i=1:length(t);
    subplot(2,1,1), plot(d,vdt(i,:),'r',d,abs(Vd(1,:)),'r--',d,-abs(Vd(1,:)),'r:')
    axis([0,1,-2,2]); set(gca,'Xdir','reverse'); grid on;
    daspect([0.1 1 1]);
    xlabel('d[m]'); ylabel('v(d,t)[V]');
    subplot(2,1,2), plot(d,idt(i,:),'b',d,abs(Id(1,:)),'b--',d,-abs(Id(1,:)),'b:')
    axis([0,1,-2,2]); set(gca,'Xdir','reverse'); grid on;
    daspect([0.1 1 1]);
    xlabel('d[m]'); ylabel('i(d,t)[A] (\times Z_o)');
    M(:,i)=getframe(h);
end
disp('Tapez sur une touche');
pause
close(h);
disp('Sortez le pop-corn');

figure(1)
clf;
axis off;
disp('période de 1 seconde, longueur d''onde de 1 mètre');
movie(M,10,length(t)-1,[-70 -50 0 0])
disp('FIN :-(');