A generalized numerical method for the fit of regression curves to experimental data is given. It is an ordered iteration procedure that has enabled representation of data by non-linear forms of function containing many coefficients, including integral and non-linear differential equations that themselves require numerical solution. It requires only that the function is entered into the computer program and unlike existing techniques requires no further analysis.