A new numerical scheme has been introduced for viscoelastic fluid flows. First, it was studied the viscoelastic behavior with particular emphasis on equations in wich the extra stress tensor is not explicity dependent of the velocity field. Two of the most difficult to solve constitutive models have been considered, the Maxwell-B and Oldroyd-B. The numerical scheme is based upon a mixed finite element approximation satisfying the Babuska-Brezzi condition for the three-fields. The adopted approximation was designed in order to optimize computer efforts. To solve the non-linear system of equation it was chosen the incremental Newton method. The SUPG (Streamline Upwind Petrov-Galerkin) technique was chosen to stabilize the hyperbolic characteristic of the constitutive law. Various numerical tests confirming the good performance of the numerical scheme are introduced.