We implement a Gauss-Newton method for fitting a specified nonlinear model to observed values by using iterative linear-regression steps to reduce the sum-of-squares error between the model and observations. A brief example of fitting observed laser energy measurements to a Gaussian model is provided.