Flux point method

Implements #50 .
We removed this before as the code just calls into scipy - is it still worth including this?