ToolPythonBeginner Units with pint: catch the mistake that lost the Mars Climate Orbiter E_run = cumulative_trapezoid(q, t, initial=0)E_run = Q_(cumulative_trapezoid(q.m_as("W"), t.m_as("h"), initial=0), "W*h").to("kWh")
ToolPythonIntermediate findiff.PDE with mixed boundary conditions: seepage under a dam psi = cumulative_trapezoid(-k * dhdx, y, axis=1, initial=0) # flow between bedrock and node