function dplanck(wav,temp)
	implicit double precision (a-h,o-z)
	parameter (one=1.d0,c1=1.192d-5,c2=1.439565d0)
	dum = c2*wav/temp
	dplanck = dum/(one-dexp(-dum))
	return
        end