Solving nonlinear differential third order equation using Python

Question

I would like to solve a nonlinear third order differential equation using Python.

In my case it is :

d^3f/dx^3 = (1-f)/(f^3)

I wrote the following progr