A numerical method for the kinematic analysis of planar higher pairs in rolling contact