updated CompareNbodyToAnalytic.py for more general comparisons with variable input data
Added global variable for designating the maximum j-value to include when computing the 0/1F-term Laplace coefficients
Numerous improvements to printed output plus fixed issue with coordinate transformation
modified plotting methods and fixed the N-body coordinate transformation method for case when input is given in the MCMC parameter form