| Hosted by CoCalc | Download
def f(beta): return gamma(5/beta)*gamma(1/beta)/(gamma(3/beta)^2)
f(1)
6
%time plot(f, .5, 2, plot_points=50, gridlines='minor')
CPU time: 0.41 s, Wall time: 0.43 s
def pdf(b): x = var('x') return (1/(2*b)) * exp(-abs(x)/b) plot(pdf(1), -10,10)