The natural factor formulation of the stiffness for the matrix displacement method