|
# normalized lower incomplete beta integral ibeta(a,b, x)
save_encoding = GPVAL_ENCODING
set encoding utf8
set title "Incomplete beta integral {Γ(a+b)/(Γ(a)Γ(b))} {/*2 ∫@_{/*0.5 0}^{/*.5 x}} {/:Italic t^{a-1}(1-t)^{b-1}dt}\n" offset 0,-1
set xrange [0:1]
set yrange [0:1]
set xlabel "x"
set ylabel "ibeta(a,b, x)"
set key outside
plot ibeta(9,1,x) lw 2 title "a=9 b=1", ibeta(7,2,x) lw 2 title "a=7 b=2",\
ibeta(5,5,x) lw 2 title "a=5 b=5", ibeta(2,7,x) lw 2 title "a=2 b=7",\
ibeta(1,9,x) lw 2 lt 6 title "a=1 b=9"
Click here for minimal script to generate this plot
|