ToolPythonIntermediate Crank-Nicolson for the Schrödinger equation: tunneling through a barrier s = np.where(E < V0, np.sinh(q * a), np.sin(q * a))
ToolPythonIntermediate Numba from the ground up: a Monte Carlo Ising model a hundred times faster m[T < TC] = (1 - np.sinh(2 / T[T < TC]) ** -4) ** (1 / 8)