jacobi
in nr2.h
extern void jacobi(float** a, int n, float d[], float** v, int* nrot);