".times..function. ##EQU00001## where k is a spring constant connecting the nodes i and j, x.sub.ij=x.sub.j-x.sub.i, and x.sub.ij* is a vector that is calculated at each time step but regarded as constant when differentiating, wherein at least" . . . .