! table of angles of incidence and exitance for each surface for a specific ray
! Georg Nadorff    4-apr-2002

input "Hx: ", hx
input "Hy: ", hy
input "Px: ", px
input "Py: ", py

print
format 5.2
print "Ray Hx Hy Px Py:", hx, hy, px, py
print "     surf      aoi        aor"
print "==============================="
format 10.4
!trace chief ray for entire system:
raytrace hx, hy, px, py


for s = 1, nsur(), 1


normal = 57.29577951*atan(rany(s)/ranz(s))
slope  = 57.29577951*atan(raym(s-1)/rayn(s-1))
aoi = (slope - normal)

slope  = 57.29577951*atan(raym(s)/rayn(s))
aor = (slope - normal)

print s, aoi, aor

next