vector package



On 06.08.2012 20:22, Mario Rodriguez wrote:
> El 06/08/12 18:40, Adam Majewski escribi?:
>>
>> Thx. Like this ?
>>
>>
>>
>> pp:[points([[1,0],[1,1]]), points([[1/2,0],[1/2,1]] )]$
>
> Yes. This is exactly what I meant.
>
> --
> Mario
  Here is full code using different definition :

/* ------ code -------- */
/*

    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

*/


kill(all);
remvalue(all);


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

/* compute   */
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 )$