// La fonction times(path p, real x=a) bien pratique pour connaitre le 
// "time" des points d'intersection de p avec la droite d'équation x=a.

import geometry;
size(8cm,0);
// Pour matérialiser le repère :
draw((1,0)--(0,0)--(0,1),Arrows());
// Tracé d'un chemin(path) p et placement des points le définissant.
path p=(-.5,0)..(2,0)..(2,1)..(0,2)..(-1.5,1)..(-0.5,1.5)..cycle;
draw("p",p);
dot(p,2bp+black);
for(int k; k<length(p); ++k) label(string(k),point(p,k),unit(accel(p,k)));

real a=-1, b=-.5, c=2;
real[] ta=times(p,a),
       tb=times(p,b),
       tc=times(p,c);
for(int k=0; k<ta.length; ++k) dot(point(p,ta[k]),3bp+red);
for(int k=0; k<tb.length; ++k) dot(point(p,tb[k]),3bp+blue);
for(int k=0; k<tc.length; ++k) dot(point(p,tc[k]),3bp+green);
draw(Label(format("$x=%f$",a),BeginPoint,S),line(1,0,-a),dashed+red);
draw(Label(format("$x=%f$",b),EndPoint,N),line(1,0,-b),dashed+blue);
draw(Label(format("$x=%f$",c),BeginPoint,S),line(1,0,-c),dashed+green);

addMargins(0,1cm);