In this paper, we have considered supersymmetry algebra for pseudoharmonic potentials with an arbitrary radial angular momentum. We have obtained the bound states and wave functions by the associated Laguerre differential equation. We have factorized Schrödinger equation for pseudoharmonic potential in terms of first order differential equations, where these operators can make corresponding partner Hamiltonian. Therefore, by introducing superpotential, we obtained the pseudo-harmonic potential which gives rise to existence of shape invariance condition. In that case, we have calculated supercharge algebras. Finally, we have computed energy eigenvalues for several diatomic molecules.