  We discuss supersymmetric biorthogonal systems, with emphasis given to the
periodic solutions that occur at spectral singularities of PT symmetric models.
For these periodic solutions, the dual functions are associated polynomials
that obey inhomogeneous equations. We construct in detail some explicit
examples for the supersymmetric pairs of potentials V_{+/-}(z) = -U(z)^2 +/-
z(d/(dz))U(z) where U(z) = \sum_{k>0}u_{k}z^{k}. In particular, we consider the
cases generated by U(z) = z and z/(1-z). We also briefly consider the effects
of magnetic vector potentials on the partition functions of these systems.
