Assistance Request with the Diffuse Approximation Model (SRDA)

[ImanKafian]

Hi,
I am trying to use the SRDA class to estimate the reflectance value using the diffuse approximation method. To this end, I have two questions:

1. Could you please introduce any literature support for the PyXopto-specific implementation of the diffuse approximation? An implementation of the diffuse approximation which is frequently seen in the literature is just the "_Rf" method of the SRDA class (shown below). I would like to know how the other terms are contributing to the estimated values of the reflectance in SRDA.

2. Could you please explain how we can use the SRDA class to estimate the reflectance for a given set of (mua, musp, r, nmed, noutside)-values? My intention is to use it for estimating the reflectance in a fiberoptic configuration with an acceptance angle of 60 degrees and a separator distance of 2 mm . As shown below, I am getting a value of 11.755 which is exceeding the maximum possible value for reflectance (=1).

I couldn't understand what is the unit of the "r" argument in the call dunder method and assumed it should be in "m".