Motivated by the experimental measurement of the decay rate, Γ, and the longitudinal polarization, PL, in the Cabibbo favored decay Ds+→ϕρ+, we have studied theoretical prediction within the context of factorization approximation invoking several form factors models. We were able to obtain agreement with experiment for both Γ and PL by using experimentally measured values of the form factors A1Dsϕ(0), A2Dsϕ(0) and VDsϕ(0) in the semi-leptonic decay Ds+→ϕl+νl. We have also included in our calculation the effect of the final state interaction (fsi) by working with the partial waves amplitudes S, P and D. Numerical calculation shows that the decay amplitude is dominated by S wave, and that the polarization is sensitive to the interference between S and D waves. The range of the phase difference δSD=δS−δD accommodated by experimental error in PL is large.