/* 
   Ergodic Theory and Dynamical Systems (1998), 18 : pp 739-758
    1998 Cambridge University Press
    Describing quadratic Cremer point polynomials by parabolic perturbations
   DAN ERIK KRARUP SØRENSEN 

  Maxima CAS code 

*/


kill(all);
remvalue(all);

/* ------- functions --------------- */

GivePoints(i):= points([[1/(2^i) ,0],[1/(2^i),1]])$

compile(all);

/* ----- constans and variables     --- */

iMax:100;
p:[];   /* list with points object */
pp:[]; /* list of points objects  */

/* --------- main -------------*/

pp:[points([[0,0],[1,0]]), points([[0,0],[0,1]] )]$

for i:0 thru iMax step 1 do
(
 p:GivePoints(i),
 pp:cons(p,pp)
)$

/* ---------- draw image  ----------- */
load(draw);

draw2d(
 title     = "Topologist\\47s comb", /* Syntax for postscript enhanced option */
 xrange        = [-0.2,1.2],
 yrange        = [-0.2,1.2],
 file_name = "comb_p",
 terminal      = screen,
 point_type = filled_circle,
 point_size = 0.5,
 points_joined = true, /* draws segments between points  */
 color         = "blue",
 pp )$