%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)

report(t)