Fixes for willard chandler surfaces

The original implementation I did of this was numerically unstable because
a minimisation rather than a root finding algorithm was being used.  In this
new implementation I have also changed it so that you can specify an upper
and lower bound in your search for a surface