ludcmp
in nr2.h
extern void ludcmp(float** a, int n, int* indx, float* d);