;pro plot_adc_lin, ps = ps
;restore, 'LOG10-09-09B.sav'
; loads date, and data


; Upward slope

ndata = n_elements(data)
err = sqrt(data + 1.) ; assume Poisson errors
index = findgen(ndata)


; Calculate the starting position
slopeStart = ( min(data) - max(data) ) / ndata
interceptStart = min(data)

start = [slopeStart, interceptStart]

result = MPFITFUN('MYLIN', index, data, err, $
                  start, yfit = yfit, bestnorm = best, dof = dof, /quiet, $
                 perr = dp)

residuals = data - yfit

end