lambdaEnd = lambdaStart + (nObs-1)*lambdaDelta
lambda = (lambdaStart:lambdaDelta:lambdaEnd)'
s=spectra(:,6)
loglog(lambda,s,".-")
xlabel("Wavelength")
ylabel("Intensity")
lambdaEnd = lambdaStart + (nObs-1)*lambdaDelta
lambda = (lambdaStart:lambdaDelta:lambdaEnd)'
s=spectra(:,6)
loglog(lambda,s,".-")
xlabel("Wavelength")
ylabel("Intensity")