Определить интенсивность излучения на расстоянии X от источника, расположенного между двумя плоскостями Z=0 и Z=1 в точке X=0, Y=0, Z=0.5. Плоскости отражают излучение по закону Ламберта с коэффициентом отражения R=0.5. Задача трехмерная. Найти

Задача Соболева В.В заключается в определении интенсивности излучения на расстоянии X от источника, расположенного между двумя плоскостями Z=0 и Z=1 в точке X=0, Y=0, Z=0.5. Плоскости отражают излучение по закону Ламберта с коэффициентом отражения R=0.5. Задача трехмерная. Найти распределение интенсивности на прямой X =[ -8; 8] при Y=0 и Z=0.5 . При столкновении фотона с плоскостями Z=0 и Z=1 происходит отражение фотона с коэффициентом R . Принять во внимание вес фотона ph_weight(n) = ph_weight(n)*R (фотоны не поглощаются). Рассчитывается вклад фотона из точки его столкновения с плоскостью с учетом веса, квадрата расстояния от этой точки до каждой точки, в которых ищется решение (учет изменения телесного угла), и квадрата косинуса между соединяющими их прямыми (т.к. и источник и приемник ламбертовские получается зависимость четвертой степени): E_loc = E_loc + ph_weight(n)./(vecnorm(r_loc-r)).^4 *r(3);
Реализуется ламбертовское отражение от плоскости, вес фотона домножается на R, находится координата следующего столкновения и т.д., пока вес фотона не станет меньше 10-4 и пока не закончатся фотоны. Облученность, создаваемая в результате многократных отражений, равна E_loc/N/(4*pi). Реши эту задачу методом Монте Карло и предложи код на Python.