Professor Puntje schreef: ↑di 18 mei 2021, 10:43
\( \frac{\mathrm{d}x^2}{\mathrm{d}t^2} = \frac{(1 - \frac{r_s}{r})^2}{ \frac{r_s}{r} \frac{x^2}{r^2} + 1 - \frac{r_s}{r} } \cdot c^2 \,\,\,\,\,\,\,\, (6) \)
Die formule leidde ik ergens helemaal in het begin van dit topic af nog
voordat ik Huygens' principe toepaste. Hier een Python-progje en bijbehorend plotje van die functie:
Code: Selecteer alles
import numpy as np
import matplotlib.pyplot as plt
def aanw(x):
r0 = 7e8
rs = 2.95e3
c = 3e8
r = np.sqrt(x**2 + r0**2)
A = (1 - rs/r)**2
B = (rs/r)*((x**2)/(r**2)) + 1 - rs/r
C = c**2
return (A/B)*C
fig, ax= plt.subplots(figsize=(15, 15))
x=np.linspace(-1e10,1e10, 1000)
y=aanw(x)
ax.plot(x,y)
Zodat Huygens' principe vrijuit gaat - de twee pieken moeten ergens anders door komen...