;HW2: Peak Flux Program
fits_read,'M7.fits',flux
x=findgen(n_elements(flux))*0.5+500
peakw=x[where(flux eq max(flux))]
print,"Peak flux: "+strcompress(max(flux))+" at "+strcompress(peakw[0])+" angstroms."
peakfciv=max(flux[where(x ge 1500 AND x le 1600)])
peakwciv=x[where(flux eq peakfciv)]
print,"Peak of CIV: "+strcompress(peakfciv)+" at "+strcompress(peakwciv[0])+" angstroms."
end