Page 64, Section 10.3. lsolve should be replaced with utsolve. The
utsolve2.f
file has the corrected program.