Flux point method

Implements #50 .
Put the `spectral_index` option at the end together with the `model` option.