function costf = f_TS_no_offset(x0, spectra, psiT, psiS, wavel);


    costf = sum(abs(spectra - psiT.*x0(1) - psiS.*x0(2)));

return