load_pha("3c273.pi")
subtract()
notice(0.4, 6.0)
set_source(xsphabs.abs1 * powlaw1d.p1)
abs1.nh = 0.07
freeze(abs1)
guess(p1)
fit()

from sherpa_contrib.chart import *

plot_chart_spectrum()
plot_chart_spectrum(elow=1.0, ehigh=8.0)
save_chart_spectrum("source_flux_chart.dat", elow=1.0, ehigh=8.0)

quit