Eksempel 15.5 fjærarbeid.py

# Program fjærarbeid.py

# Definerer størrelser
def F(x): # fjærkraft [N]
  k = 300 # fjærkonstant [N/m]
  return k*x

x1 = 0 # x startverdi [m]
x2 = 0.05 # x sluttverdi [m]
x = x1 # aktuell x-verdi [m]
n = 200 # antall intervaller
dx = (x2-x1)/n # x-intervall [m]
W = 0 # W arbeid startverdi [J]

# Summerer opp ΔW-rektanglene
while x <= x2: 
  W = W + F(x)*dx # summering av arbeid [J]
  x = x + dx # neste rektangelposisjon [m]

# Skriver ut arbeidet W
print(f'W = {W:.4f} J') # utskrift med 4 desimaler

W = 0.3731 J