fluxratio2beta.pro 205 Bytes
FUNCTION fluxratio2beta,fluxratio,lambda1,lambda2,Tdust

hck=1.44e4
beta=(alog(fluxratio)+alog((exp(hck/lambda1/Tdust)-1.)/(exp(hck/lambda2/Tdust)-1.)))/ $
      alog(lambda2/lambda1)-3.

RETURN,beta

END