Added looptools functionality
Added debug msg
Modded b0 to produce more accurate answers - to the 10^-5 level when x->1
Various debuggings to do with the function b1
Provided a simple case where my implementation of b1 screws up
Added bisection routine and fixed predMZsq initialisation issue
Added stuff to print out bisected solutions for a certain parameter point
Added different functionality to your bisections: different mu signs etc
Added functionality to m0scan: can change sign of mu now
Added negative mu option to m0scan