(* Copyright (c) Dror Bar-Natan 1994 *)

f[x_,y_] := Sin[20(x^2+y^2)]/(1+5(x^2+y^2))  (* the function drawn *)

Show[Graphics[Table[ 	    (* loop 500 times *)
  x=Random[Real,{-1,-0.6}]; (* generate a random number between -1 and -0.6 *)
  y=Random[Real,{-1,1}];    (* generate a random number between -1 and 1 *)
  Table[		    (* loop 6 times *)
    x0=x; x = x + 0.4 + f[x,y]/15;  (* compute coordinates of next pixel *)
    Point[{x0,y}],	    (* put a point *)
  {6}],			    (* end inner loop *)
{500}]]]		    (* end outer loop *)
