%!PS-Adobe-2.0
%%BoundingBox: 0 0 200 200

(lenses.inc) run

20 20 translate

1 setlinecap
1 setlinejoin

72 dup scale
1 72 div setlinewidth

gsave
0.7 setgray
newpath
0 0 2 0 90 arc
stroke
grestore

/N 10 def
/h 0 def
/dh 0.99 N div def

/D { 4 dict begin
/h exch def
h 1.33 div asin 4 mul h asin 2 mul sub
end } def

newpath
N 1 add {
0 0 moveto
/d h D def
d cos 2 mul d sin 2 mul lineto
-4 0 rlineto
/h h dh add def
} repeat
stroke

   showpage