Recently, the simplified spherical harmonics equations (SPN) model has attracted much attention in modeling the light propagation in small tissue geometries at visible and near-infrared wavelengths. In this paper, we report an efficient numerical met