Source : Symetriquedtesecaxe.mp

input geometriesyr16;
vues:=20;
pair A,B,E,F,C,D,M[],N[];
A=u*(1,4);
B=u*(8,5);
C=u*(5,1);
D=u*(6,6);
figure(0,0,11u,8u);
path dd,de;
dd=droite(A,B);
de=1.1[B,A]--1.4[A,B];
  path d;
  d=droite(C,D);
  for k=0 upto vues:
    M[k]=point(k*length de/vues) of de;
    N[k]=symetrie(M[k],C,D);
  endfor;
fin;
marque_p:="plein";
_nfig:=0;
for k=0 upto vues:
  figure(0,0,11u,8u);
  trace feuillet withcolor blanc;
  trace dd;
  trace d;
  label.rt(btex $(d)$ etex,C);
   nomme.top(A);
  nomme.top(B);
  drawoptions(withcolor vert);
  pointe(A,B);
  drawoptions();
  for j=0 upto (k-1):
    pointe(M[j],N[j]);
  endfor;
  drawoptions(withcolor rouge);
  pointe(M[k],N[k]);
  drawoptions();
  trace perpendiculaire(C,D,M[k]) dashed evenly;
  trace codeperp(M[k],iso(M[k],N[k]),C,8);
  trace coupdecompas(iso(M[k],N[k]),M[k],10);
  trace coupdecompas(iso(M[k],N[k]),N[k],10);
  trace codesegments(M[k],iso(M[k],N[k]),iso(M[k],N[k]),N[k],2);
fin;
endfor;
end