##// END OF EJS Templates
use better integration boundaries...
use better integration boundaries Because this example intends to illustrate the difference between trapezoid rule and high-order integration, choose an appropriate integration interval where such a difference exists. Also implement the trapezoid rule directly for pedagogical purposes, instead of relying on a canned version from scipy.

File last commit:

r12797:d1ebf45a
r19878:9f0fa954
Show More
ipython_config.py
20 lines | 559 B | text/x-python | PythonLexer
c = get_config()
app = c.InteractiveShellApp
# This can be used at any point in a config file to load a sub config
# and merge it into the current one.
load_subconfig('ipython_config.py', profile='default')
lines = """
from __future__ import division
from sympy import *
x, y, z, t = symbols('x y z t')
k, m, n = symbols('k m n', integer=True)
f, g, h = symbols('f g h', cls=Function)
"""
app.exec_lines.append(lines)
# Load the sympy_printing extension to enable nice printing of sympy expr's.
app.extensions.append('sympy.interactive.ipythonprinting')