coupling_angles.m  - calculates parameters of the beam to be coupled into the disk via prism
beam_trace.m  - calculates beams propagation and their branching on dielectric interfaces (faces) for arbitrary set of beams and faces.
example_beam_trace.m - example use of above routine.

rest of the code: helper functions