The system stiffness matrix K is square since the vectors R and r have the same size.
Subsequently, the members' characteristic forces may be found from Eq.(1) where \mathbf{q}^m can be found from r by compatibility consideration.