Solving a Boundary Value Problem with the Finite Difference Method - Python

Assume f is a function of x and satisfies the following differential equation:
f''' + f = 1
With conditions f(0) = 0 and f(pi/2) = 0.
Here's how to solve this with the finite difference method.
Python code here
trinket.io/glo...