Неверный вывод графика

Question

Я пытаюсь построить следующее в gnuplot в Windows,

    reset

# png
set terminal wxt size 650,550 enhanced font 'Verdana,10'
set output 'example_plot.png'

set border linewidth 1

#unset key

set xtics auto
set ytics auto

set title 'example plot of 3-resonance K-matrix'
set xlabel 'E(GeV)'
set ylabel 'sigma(1/GeV)'
set xrange [3:6]
set yrange [0:0.15]

gee1 = 1
gDD1 = 1
m1 = 3.9
me = 0.000511
mf = 1.86965

p(x) = sqrt(1-(4*mf*mf/(x*x)))
q(x) = sqrt(1-(4*me*me/(x*x)))

rho(x) = (1/(16*3.14159))*sqrt(1-((4*mf*mf)/(x*x)))

P(x) = ((gee1*gDD1)/(x*x-m1*m1))
K(x) = rho(x)*((gDD1*gDD1)/(x*x-m1*m1)

Fsq(x) = (P(x)*P(x))/(1+(K(x)*K(x)))

f(x) = (3/(64*3.14*x*x))*(p(x)/q(x))*Fsq(x)
plot Fsq(x) title "T^2" with lines linestyle 1
#plot f(x) title "g(s)" with lines linestyle 1

, но получаю неправильный график.Пик должен быть около 31300, но вместо этого мой пик достигает 200.

Я действительно не вижу, что не так с моим кодом.Может ли кто-нибудь помочь?

Спасибо, Кристиан введите описание изображения здесь

Answer 1

Пик в вашей функции настолько резок, что, если вы не оцените его непосредственно по координате х, соответствующей максимуму, график не уловит истинный максимум. Чем более точно вы сэмплируете функцию построения графика (например, set sample 5000), тем больше вероятность того, что один из сэмплов будет близок к истинному максимуму, но вы не можете гарантировать его точное попадание.