FUNCTION bb_specf,t h=float(6.62606957e-34) ;physical constants c=float(299792458) k=float(1.3806488e-23) b=2*h*c^2 d=h*c/(k*t) spectrum=fltarr(1001,2) for i=0,1000 do begin spectrum(i,0)=i+1000 w=double(i+1000)*1e-10 spectrum(i,1)=b/((exp(d/w)-1)*w^5) endfor return,spectrum end