This contribution deals with the solution of geometrically nonlinear elastic problems solved by the least-squares mixed finite element method (LSFEM). The degrees of freedom (displacements and stresses) will be approximated using suitable spaces, namely W 1,p with p > 4 and H(div, Ω). In order to define the stress response of the material, different hyperelastic free energy functions will be presented. The residual forms R i of the balance of momentum and the constitutive equation build a system of differential equations of first order. Choosing suitable weighting operators and applying L 2 -norms lead to a leastsquares functional F (P ,u). The interpolation of the unknowns is accomplished using a standard polynomial interpolation for the displacements and vector-valued Raviart-Thomas functions for the approximation of the stresses. The formulations presented will be compared considering a uni-axial spatial tension test.