Included linestyles to help differentiate lines when printed BW.

It works. I tested it.
Wrote a placeholder `flux_trapezoidal` function in case I need it later.
Now we integrate flux*frequency in `integrated_flux`.

Upon testing, it gives a way better result when plotted versus tau, which I am confused by but I'll roll with it.