donnees=get(gca,'CurrentPoint');
xo=donnees(1,1);
yo=donnees(1,2);
set(obs,'XData',xo,'YData',yo)

for li=1:nbli
    for co=1:nbco
        xi=x(li,co);
        yi=y(li,co);
        dx=xo+(xi-xo).*(mur-yo)./(yi-yo);
        dx=s.*(dx-xo)+xo;
        dz1=mur+(oh.*(mur-yo))./(yi-yo);
        dz2=((h-oh).*(mur-yo))./(yi-yo);
        set(piq(li,co),'XData',[dx dx],...
		'YData',[dz1 dz1+dz2]);
    end
end
axis([0 nbco 0 nbli+mur+1])
drawnow
