Hello again. When using the subroutine ludcmp I am getting error with the matrix indices, in particular this statement is giving me the problem,


if (dum.ge.aamax) then
imax = i
aamax = dum
endif

if this statement is not met imax is never defined causing the first index of a to be too high for the following statement

dum = a(imax, k)
a(imax,k) = a(j,k)

any thoughts as in how to avoid this error?

Thanks in advance.