%!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

(lenses.inc) run

180 180 scale
0.005 setlinewidth
2 2 translate

%%%%%%%%%%%%%%%%%%%%%%%%%%Lens (R=1 and R= inf)%%%%%
newpath
0 0 1 140 220 arc 
-1.55 0 1 -40 40 arc
closepath 
gsave
0.8 0.8 1 setrgbcolor
fill
grestore 
stroke



0 setlinewidth

newpath
-3 0 moveto
3 0 lineto
0 0 0 setrgbcolor
stroke


0 0 0 setrgbcolor
/y 0.53 def
/dy 0.53 4 div def

9 { 
/p0 [-4 y] def
/v0 [1 0] def

/p1 p0 v0 [0 0] 1 hit def
/n2 [0 0] p1 gradient normalized def
/v1 n2 v0 1.5 refraction normalized def

/p2 p1 v1 [-1.55 0] 1 hit def

/n3 [-1.55 0] p2 gradient normalized def
/v2 n3 v1 1 1.5 div refraction normalized def
/p3 p2 v2 [1 0 -10 3 div] line-intersection def


newpath
p0 aload pop moveto
p1 aload pop lineto
p2 aload pop lineto
p3 aload pop lineto

stroke

/y y dy sub def

} repeat
0 0 0 setrgbcolor
-1.4 -1 moveto
/Times-Roman findfont 0.14 scalefont setfont
(Spherical Aberration) show
showpage