;perturb.pro

i=[.0095,.01758,.00563,.08564,.03905]
tens=FLTARR(100)
fens=tens
FOR j=0,99 DO BEGIN
    di=0.001*RANDOMN(s,5)
    i2b,i+di,w,b,db
    bbfit,w,b,db,t,f,/noplot
    tens(j)=t
    fens(j)=f*1E10
    ENDFOR
PLOT,tens,fens,PSYM=1,XTITLE='Temperature (K)' $
    ,YTITLE='Attenuation Factor * 1E10'
END