CoCalc Shared Filescode / wstein-fitting.sagewsOpen in CoCalc with one click!
Authors: Harald Schilly, William A. Stein, Nicholas Triantafillou
%auto %attach compute_lambda.sage
plot(log(x/5)^2, 0, 5)
t.plot_histogram(bins=400)
def rmt_pdf(b): x = var('x') return b*abs(log(abs(x))^2)
t = get_drew_data('11a', 3) @interact def f(b=.1): show(plot(rmt_pdf(.1), -2,2, color='red', thickness=1) + t.plot_histogram(bins=400), ymax=10)
Interact: please open in CoCalc
report(t)
kurtosis: 10.5467804636 mean: 0.00426001527055 sd: 1.77988021702