def f(beta): return gamma(5/beta)*gamma(1/beta)/(gamma(3/beta)^2)
f(1)
%time plot(f, .5, 2, plot_points=50, gridlines='minor')
https://en.wikipedia.org/wiki/Laplace_distribution
def pdf(b): x = var('x') return (1/(2*b)) * exp(-abs(x)/b) plot(pdf(1), -10,10)