Have added option of using b=aRs*cos(i) instead of the
inclination directly in the transit.ma02_aRs() routine.