Flux point method

Implements #50 .
Put the `numpy` import to the other numpy import at the top.