Provides functions for fast matrix inversion, creation of empirical CDF from sample data including
handling of asymmetric errors, and fitting to a funtion using chi-squared. The fitting procedure return
the full covariance matrix describing the fitted parameters.