Как высчитать распределение интенсивности света в фокальной плоскости линзы?

Нужно смоделировать распределение интенсивности сфокусированного сферической линзой (конечно со сферической аберрацией) гауссового (или хотя бы просто параллельного) пучка на экране вблизи фокальной плоскости фокусирующей линзы. С чего начать копать?
У меня мысли такие, что нужно рассчитывать интенсивность в каждой точке экрана чтобы создать "картинку" распределения интенсивности. Тогда как это сделать?
Заранее спасибо! :)