CoCalc Shared Filescode / william-crazy-distributions.sagewsOpen in CoCalc with one click!
Authors: Maarten Derickx, JB Nam, James Rickards, William A. Stein, Jan Vonk
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

Laplace

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)