result, error = quad(integrand, 0, 4) print(result)
: NR is written for procedural/compiled languages; naive Python loops are too slow. numerical recipes python pdf