  Hereby we complete the proof of integrability of the Lax systems, based on
pseudo-Riemannian coset manifolds G/H^{*}, we recently presented in a previous
paper [arXiv:0903.2559]. Supergravity spherically symmetric black hole
solutions have been shown to correspond to geodesics in such manifolds and, in
our previous paper, we presented the proof of Liouville integrability of such
differential systems, their integration algorithm and we also discussed the
orbit structure of their moduli space in terms of conserved hamiltonians. There
is a singular cuspidal locus in this moduli space which needs a separate
construction. This locus contains the orbits of Nilpotent Lax operators
corresponding to extremal Black Holes. Here we intrinsically characterize such
a locus in terms of the hamiltonians and we present the complete integration
algorithm for the Nilpotent Lax operators. The algorithm is finite, requires no
limit procedure and it is solely defined in terms of the initial data. For the
SL(3;R)/SO(1,2) coset we give an exhaustive classification of all orbits,
regular and singular, so providing general solutions for this case. Finally we
show that our integration algorithm can be generalized to generic
non-diagonalizable (in particular nilpotent) Lax matrices not necessarily
associated with symmetric spaces.
