Flux point method

Implements #50 .
Put imports from astropy and gammapy at the top.