ConceptPythonIntermediate Symplectic integrators: why leapfrog's energy error is bounded and RK4's drifts print(f"angle per step: {np.arccos(1 - h**2 / 2) / h:.6f} h")
ToolPythonIntermediate findiff.PDE with mixed boundary conditions: seepage under a dam harr = H / np.pi * np.arccos(2 * x_base / B)gap = np.abs(hb - H / np.pi * np.arccos(2 * xb / B)).max()