Solution to the intrinsic method for the evaluation of beam shape coefficients in spheroidal coordinates