[…] The method of determining an optimal step size, is a little complicated. Therefore, the following implementation of Newton's method always takes a full step. That is, it is not the damped method. Because of this, the output of this function will not always match the output of the NLPNRA subroutine. However, this function implements the main ideas behind Newton's method in higher dimensions. Because it is tedious to write modules that evaluate the gradient and Hessian of the objective function, this implementation uses the NLPFDD function in SAS IML to obtain finite-difference derivatives. …