How to do symbolic differentiation in Python with SymPy
How to perform numerical integration using the Trapezoidal, Midpoint or Simpsons Rule in Python