# Matplotlib compatibility patch for Pyodide
import matplotlib
if not hasattr(matplotlib.RcParams, "_get"):
    matplotlib.RcParams._get = dict.get

… using differential equations

\[ -GI \cfrac{{d}^2\varphi_x}{{d}x^2} = -q_{t,z}\]

This is not treated in the books.