The integration is still accurate for a larger step cap
Further speed up the code by using larger step cap
Bug fix: missed a factor of 1/4pi in L_j_ff()
Add references to the emissivities
Setup and cleanup IDL environment for plotting to EPS files
Add the IDL scripts for interstellar scattering