Added script to test the parameter estimation part of the code - it is almost working, but some issues still require sorting out