  We consider the issue of deriving superintegrable systems with position
dependent mass (PDM) in two dimensions from certain known superintegrable
systems using the recently introduced method of master symmetries and complex
factorization by M. Ranada \cite{Rana1,Rana2,Rana3,Rana4}. We introduce a
noncanonical transformation to map the Hamiltonian of the PDM systems to that
of ordinary unit mass systems. We observe a duality between these systems. We
also study Tsiganov's method \cite{Tsiganov1,Tsiganov2,Tsiganov3,Tsiganov4} to
derive polynomial integrals of motion using addition theorems for the
action-angle variables using famous Chebyshev's theorem on binomial
differentials. We compare Tsiganov's method of generating an additional
integral of motion with that of Ranada's master symmetry method.
