x = np.linspace(-2.0, 2.0, 100) y1 = 1 / (4 * 2) * x**2 # Parabola with focal length f = 2 y2 = 1 / (4 * 6) * x**2 # Parabola with focal length f = 6 plt.figure ...