Flux point method

Implements #50 .
Can I add this as a separate unit test, or as a second part to the existing test_compute_differential_flux_points() unit test?