Added a function to get the phase interference pattern from beams on a surface and return it as a numpy array