This paper presents a relative coordinate formulation for constrained mechanical systems. Relative coordinates are defined along degrees of freedom of a joint. Graph theoretic analyses are performed to identify topological paths in mechanical systems. Cut constraints are generated to handle closed loop systems. Equations of motion are derived in the Cartesian space and transformed to the joint space. Relative generalized coordinates are corrected to satisfy the cut constraints by a parametrizatiom method.