Afstandskvadratlov

import numpy as np

import matplotlib.pyplot as plt

# Funktion til at beregne intensitet baseret på afstand (i AU)

def intensitet(kilde_intensitet, afstand):

    return kilde_intensitet / (afstand ** 2)

# Konstant for solens intensitet ved 1 AU (solkonstanten ved Jorden)

solens_intensitet = 1361  # Watt per kvadratmeter på Jorden (1 AU)

# Afstande fra Solen til planeterne (i astronomiske enheder, AU)

afstande_planeter = {

    'Merkur': 0.39,  # AU

    'Venus': 0.72,   # AU

    'Jorden': 1.00,  # AU

    'Mars': 1.52,    # AU

    'Jupiter': 5.20  # AU

}

# Afstande i AU for grafen

afstande_au = np.linspace(0.36, 6, 500)  # Fra 0.1 AU til 6 AU

# Beregn intensiteterne ved forskellige afstande fra Solen

intensiteter = intensitet(solens_intensitet, afstande_au)

# Plot resultaterne

plt.figure(figsize=(10, 6))

plt.plot(afstande_au, intensiteter, label='Solens stråling ift. afstand', color='orange')

# Marker planeterne med prikker og deres afstande fra Solen

for planet, afstand_au in afstande_planeter.items():

    intensitet_planet = intensitet(solens_intensitet, afstand_au)

    plt.scatter(afstand_au, intensitet_planet, label=f'{planet}', s=100)

    plt.text(afstand_au, intensitet_planet, f'  {planet}', fontsize=12, verticalalignment='bottom')

# Tilføj titler og labels

plt.title('Afstandskvadratloven: Solens stråling som funktion af afstand', fontsize=14)

plt.xlabel('Afstand fra Solen (AU)', fontsize=12)

plt.ylabel('Stråling (W/m²)', fontsize=12)

plt.grid(True)

plt.legend()

plt.show()