Ulam spiral

[Damien]

Posted by Gilles Carpentier on 22 July 2013, 5:38 p.m.

Code:


EXPORT Ulam
BEGIN
 LOCAL a,b,xy:={160,120}, n:=1, m:=.9, d:={{1,0}, {0,1}, {-1,0}, {0,-1}};
 RECT();
 WHILE n<100000 DO
  FOR a FROM 1 TO 4 DO
   m:=m+.5;
   FOR b FROM 1 TO m DO
    IF isprime(n) THEN PIXON_P(xy,127);END;
    xy:=xy+d(a);
    n:=n+1;
   END;
  END;
 END;
 FREEZE;
END;

another one !

Regards,

Damien.