Построение функции черного тела (кривая Планка) с использованием gnuplot
Я пытаюсь использовать gnuplot для построения стандартной кривой черного тела (функция Планка) и приближения Вены, чтобы соответствовать набору заданных звездных данных.
Данные отображаются точно в (произвольном) диапазоне от 0 до 1,4 единиц. Тем не менее, независимо от того, что я делаю, кривые Планка и Вены выходят почти на 0 везде - все точки находятся в диапазоне 1e-27.
Возможно, я неправильно определяю функции или gnuplot устанавливает нулевые значения без моего ведома?
Вот код gnuplot, который я использую:
set title "Planck and Wien functions (blackbody curve)" font "Helvetica,14"; \
set out "ps3g0ivplot.eps"
set xrange ["0":"1.1e4"]; \
#set yrange ["0":"1.4"]; \
set ylabel "Arbitrary flux" font "Helvetica,14"; \
set xlabel "Wavelength (Angstrom)" font "Helvetica,14"; \
set mxtics; \
set mytics; \
#constants
h = 6.62606876e-34 # Planck's constant
c = 2.99792458e8 # speed of light
k = 1.3806503e-23 # Boltzmann's constant
Planck(x,t) = \
(2*h*c*c/x**5)*(exp(-1*h*c/(x*k*t))-1)
Wien(x,t) = \
(2*h*c**2/x**5)*exp((-1*h*c)/(x*k*t))
#set temperature, t
t = 6400 #in Kelvin, to match given star in ps3g0iv.dat
#print Planck(4500, t)
#print Wien(4500, t)
plot \
Wien(x,t) lt 2 lw 3 title "Wien" , \
Planck(x,t) lt 1 lw 3 title "Planck", \
#"ps3g0iv.dat" lt 9 pt 1 pointsize 1 title "ps3g0iv.dat" with lines
quit