fermi_catalog = FermiCatalog()
flux_points = fermi_catalog.flux_points('2FGL XXX')


flux_points = read_flux_points('hess_crab.txt')
flux_points.write('hess_crab.fits')

compute_differential_flux_point(table, )

compute_integral_flux_point(table)

'ENERGY', 'DIFF_FLUX', 'DIFF_FLUX_ERR'

'ENERGY_MIN', 'ENERGY_MAX', 'INT_FLUX', 'INT_FLUX_ERR'


_compute_x_lafferty(f, x_min, x_max)

ebounds = []
int_flux = []
spectral_index = 2.3

flux_points = FluxPoints(ebounds, int_flux)
flux_points.calc_diff_flux(xmethod='lafferty', ymethod='powerlaw', spectral_index=spectral_index)
return flux_points




from astropy.table import Table

t = Table()
t.meta['spectral_index'] = 3.2
t['flux'] = [1, 2, 3]

print(t.meta)