matrix
in nr2.h
extern float** matrix(long nrl, long nrh, long ncl, long nch);