11.3 Matematisk programmering

def f(x):
    return x**3 - 2*x - 5

# Søk mellom 1 og 3 med liten steglengd
x = 1.0
steg = 0.001

while x <= 3.0:
    if abs(f(x)) < 0.01:
        print(f"Mogleg løysing: x = {x:.3f}, f(x) = {f(x):.4f}")
    x += steg

Utskrift:

Mogleg løysing: x = 2.094, f(x) = -0.0053
Mogleg løysing: x = 2.095, f(x) = 0.0079

Løysinga er omtrent $x \approx 2{,}094$.

Vi kan sjekka: $2{,}094^3 - 2 \cdot 2{,}094 - 5 = 9{,}19 - 4{,}19 - 5 = -0{,}005 \approx 0$ \checkmark

# Enkel tekstgraf av f(x) = x^2
for x in range(-5, 6):
    y = x**2
    stolpar = "*" * y
    print(f"x = {x:3d} | {stolpar} ({y})")

Utskrift:

x =  -5 | ************************* (25)
x =  -4 | **************** (16)
x =  -3 | ********* (9)
x =  -2 | **** (4)
x =  -1 | * (1)
x =   0 |  (0)
x =   1 | * (1)
x =   2 | **** (4)
x =   3 | ********* (9)
x =   4 | **************** (16)
x =   5 | ************************* (25)

Vi kan tydeleg sjå parabelforma: grafen er symmetrisk om $x = 0$ og stig raskt for store $|x|$.

For ein finare graf kan du bruka matplotlib:

import matplotlib.pyplot as plt

x_verdiar = []
y_verdiar = []

for i in range(-50, 51):
    x = i / 10
    x_verdiar.append(x)
    y_verdiar.append(x**2)

plt.plot(x_verdiar, y_verdiar)
plt.xlabel("x")
plt.ylabel("f(x)")
plt.title("f(x) = x²")
plt.grid(True)
plt.show()

import random

antall_forsok = 10000
antall_syv = 0

for i in range(antall_forsok):
    terning1 = random.randint(1, 6)
    terning2 = random.randint(1, 6)
    if terning1 + terning2 == 7:
        antall_syv += 1

sannsynlegheit = antall_syv / antall_forsok
print(f"Antall gongar sum = 7: {antall_syv}")
print(f"Estimert sannsynlegheit: {sannsynlegheit:.4f}")

Typisk utskrift:

Antall gongar sum = 7: 1674
Estimert sannsynlegheit: 0.1674

Eksakt sannsynlegheit: Det er 6 gunstige utfall for sum 7: $(1,6), (2,5), (3,4), (4,3), (5,2), (6,1)$. Totalt $6 \times 6 = 36$ utfall.

$$P(\text{sum} = 7) = \frac{6}{36} = \frac{1}{6} \approx 0{,}1667$$

Simuleringen gjev eit resultat nær den eksakte verdien! \checkmark

import math

tal = int(input("Skriv inn eit tal: "))

if tal < 2:
    print(f"{tal} er ikkje eit primtal")
else:
    er_primtal = True
    for i in range(2, int(math.sqrt(tal)) + 1):
        if tal % i == 0:
            er_primtal = False
            break

    if er_primtal:
        print(f"{tal} er eit primtal")
    else:
        print(f"{tal} er ikkje eit primtal")

Triks: Vi treng berre å sjekke divisorar opp til $\sqrt{n}$, fordi viss $n = a \cdot b$ og begge er større enn $\sqrt{n}$, ville $a \cdot b > n$.

For $n = 97$: $\sqrt{97} \approx 9{,}8$, så vi sjekkar 2, 3, 4, 5, 6, 7, 8, 9. Ingen deler 97, altså er 97 eit primtal.

karakterar = [4, 5, 3, 6, 5, 4, 5, 3, 4, 5]

# Gjennomsnitt
gjennomsnitt = sum(karakterar) / len(karakterar)

# Sorterer for median
sortert = sorted(karakterar)
n = len(sortert)
if n % 2 == 1:
    median = sortert[n // 2]
else:
    median = (sortert[n // 2 - 1] + sortert[n // 2]) / 2

# Variasjonsbreidde
variasjonsbreidde = max(karakterar) - min(karakterar)

# Typetal (enkel metode)
typetal = max(set(karakterar), key=karakterar.count)

print(f"Karakterar: {karakterar}")
print(f"Sortert: {sortert}")
print(f"Gjennomsnitt: {gjennomsnitt:.1f}")
print(f"Median: {median}")
print(f"Typetal: {typetal}")
print(f"Variasjonsbreidde: {variasjonsbreidde}")

Utskrift:

Karakterar: [4, 5, 3, 6, 5, 4, 5, 3, 4, 5]
Sortert: [3, 3, 4, 4, 4, 5, 5, 5, 5, 6]
Gjennomsnitt: 4.4
Median: 4.5
Typetal: 5
Variasjonsbreidde: 3

Sjekk: Gjennomsnitt $\displaystyle = \frac{4+5+3+6+5+4+5+3+4+5}{10} = \frac{44}{10} = 4{,}4$ \checkmark

Median: Dei to midtarste er 4 og 5, median $\displaystyle = \frac{4+5}{2} = 4{,}5$ \checkmark