f(x) = 5*besj1(x)*exp(-x/5)
set samples 20
set style fill solid 0.5
set xrange [0:10]
set yrange [-1:2.5]
unset key
seed = rand(time(1))
set term pngcairo enh font "Times, 10" size 480,360
set out "plot_boxxyerrorbars.png"
set label "plot with boxxyerrorbars" at sc 0.4, 0.85 font "Helvetica,18"
set xlabel "Ｘ軸 IPAゴシック" font "IPAGothic,12"
set ylabel "Ｙ軸 IPA明朝" font "IPAMincho,12"
plot '+' u 1:(f($1)):((rand(0)+0.3)/5):(rand(0)/2.0) w boxxyerrorbars