I've incorporated the 4D polynomial fitting routines from EAM into the CVS tree and added fairly simple tests. I changed from GaussJordan to LUD since the former was overflowing.