Programmet gjør en beregning på et elektron.

from pylab import *

def f(m,c,v):
    return m*c**2/(sqrt(1 - (v/c)**2)) - m*c**2

m = 9.11e-31
c = 3.00e8
dv = 0.01*c
v = 0
v_liste = [v]
f_liste = [0]

while v < 0.99*c:
    v = v + dv
    v_liste.append(v)
    f_liste.append(f(m,c,v))
plot(v_liste, f_liste)
show()

q)

Hva viser grafen programmet skriver ut?

Fasit

q)

B

LøsningsforslagKI-generert

q)

Funksjonen returnerer

$$\frac{m{c}^2}{\sqrt{1-{\left(\frac{v}{c}\right)}^2}}-m{c}^2=\gamma m{c}^2-m{c}^2=E-E_0=E_k$$

altså totalenergien minus hvileenergien, som er den kinetiske energien. Programmet plotter f_liste mot v_liste, altså kinetisk energi som funksjon av fart. Totalenergien (D) ville vært $\gamma m{c}^2$ uten å trekke fra $m{c}^2$, og bevegelsesmengden (C) er $\gamma mv$.