Resujemo sistem enacb velikosti N. Matrika A predstavlja levo stran sistema enacb. Pri fortranu so lahko matrike le staticne in zato parameter NP opisuje velikost matrike A, ki je obicajno vecja kot pa sistem enacb (N < NP). 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 B.
Potek izracuna sistema linearnih enacb bi potekal nekako takole:
parameter (NP=100)
real a(NP,NP), dQ(NP), d
integer indx(NP)
...
...
n = pipes
...
...
call ludcmp(a, n, NP, indx, d)
call lubksb(a, n, NP, indx, dQ)
...