ConceptPythonIntermediate Finite volumes: why a conservation law is solved by bookkeeping the fluxes path_eff = np.trapezoid((u0 / speed(x_fine[path]))**2, x_fine[path])
ConceptPythonIntermediate Stiffness: why an explicit solver crawls on a reaction that has long settled predicted = np.trapezoid(1 / bound, s.t)