Higher Order Asymptotics Unleashed: Software Design for Nonlinear Heteroscedastic Regression