ToolPythonIntermediate Ray tracing a black hole with solve_ivp: the shadow of M87* ax.plot(z_out, 12 / np.cos(alpha12) - np.tan(alpha12) * z_out, color=MUTED, ls="--", lw=1) # outgoing
ToolPythonIntermediate The Sievert integral with scipy.integrate.quad: dose beside a shielded pipe x_el = np.tan(u) # pipe element seen at angle u