Flux point method

Implements #50 .
Looks like you introduced an indentation problem here and in a bunch of other places ... please review the diff on Github and fix those cases.