Monte Carlo for Spherical Harmonics
The BDRF is defined as
U angle the ray enters
V angle the ray leaves
p probablity density function
R mean attenuation of all the rays incident with U and scattered in direction V