Lower/Upper dekompozicija

Resujemo sistem enacb velikosti N. Matrika h predstavlja levo stran sistema enacb. Pri fortranu so lahko matrike le staticne in zato parameter NEL opisuje velikost matrike h, ki je obicajno vecja kot pa sistem enacb (n < NEL). INDX je celoctevilcni vektor permutacji v matriki A in se prenasa naprej tako kot parameter D v podprogram LUBKSB, kateri zahteva se desno stran sistema enacb v vektorju B. Po izracunu se rezultat nahaja v vektorju U. Potek izracuna sistema linearnih enacb bi potekal nekako takole:

     parameter (NEL=100)
     real h(NEL,NEL), g(NEL, NEL), q(NEL), u(NEL), d
     integer indx(NEL)
     ...
     ...
     n = elements
     ...
     ...
     call ludcmp(h, n, NEL, indx, d)
     call lubksb(h, n, NEL, indx, u)
     ...