function p1
    verify_order_of_convergence('Method 1', [4.0e-2, 9.1e-4, 4.8e-7], 2)
    verify_order_of_convergence('Method 2', [3.7e-4, 1.2e-15, 1.5e-60], 3.9)
    verify_order_of_convergence('Method 3', [4.3e-3, 1.8e-8, 1.4e-24], 3)
end

function verify_order_of_convergence(label, seq, order)
    disp(label)
    alpha = (log(seq(2))-log(seq(3))) / (log(seq(1))-log(seq(2)))
    lambda = exp(log(seq(2))-alpha*log(seq(1)))
end
