I have both NR in C and in C++. I took the example xsimplx.c from my old floppy and rewrote it in C++. However the returned values in izrov and iposv are different in the C version and the C++ version. Is there a bug in the support routines (I already applied the update to simp2.cc) or am I missing something. I have attached (I hope) my implementation of xsimplx.cc. I should add that I use the extention hh for C++ header files which explains the nr.hh

I should add that I add methods to VecNR and MatNR to print a vector or matrix. The methods take an ostream and a width. This was added so that I can easily debug routines.