Higher-order asymptotics unleashed: Software design for nonlinear heteroscedastic regression