  Based on a quantitative version of the inverse function theorem and an
appropriate saddle-point formulation we derive a quasi-optimal error estimate
for the finite element approximation of harmonic maps into spheres with a nodal
discretization of the unit-length constraint. The estimate holds under natural
regularity requirements and appropriate geometric stability conditions on
solutions. Extensions to other target manifolds including boundaries of
ellipsoids are discussed.
