solve (1+q)^2 * (1-q^-2) * UnitStep[q^2-1] = 28/9, for q