1 2 3 4 5 6 7 8 9
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
adding librairies...