['illusion', 'measurement', 'earth', 'meditation', 'fusion energy']
Vull que generis un codi en Python que simuli la trajectòria d'un pèndol simple amb fricció de l'aire. Mostra la posició angular del pèndol en funció del temps en una gràfica.
```python
import matplotlib.pyplot as plt
import numpy as np

# Constants del pèndol
g = 9.81  # Acceleració deguda a la gravetat (m/s^2)
l = 1.0  # Longitud del pèndol (m)
b = 0.1  # Coeficient d'amortiment

# Condicions inicials
theta_0 = 0.1  # Angle inicial (rad)
omega_0 = 0.0  # Velocitat angular inicial (rad/s)

# Paràmetres de temps
t_0 = 0.0  # Temps inicial (s)
t_f = 10.0  # Temps final (s)
dt = 0.01  # Pas de temps (s)

# Arrays per emmagatzemar els resultats
temps = np.arange(t_0, t_f, dt)
theta = np.zeros(len(temps))
omega = np.zeros(len(temps))

# Condicions inicials en els arrays
theta[0] = theta_0
omega[0] = omega_0

# Mètode d'Euler per a la integració numèrica
for i in range(1, len(temps)):
    # Calcula les derivades
    dtheta_dt = omega[i - 1]
    domega_dt = - (g / l) * np.sin(theta[i - 1]) - (b / l) * omega[i - 1]
    
    # Actualitza els valors
    theta[i] = theta[i - 1] + dtheta_dt * dt
    omega[i] = omega[i - 1] + domega_dt * dt

# Crea la gràfica
plt.plot(temps, theta)
plt.xlabel('Temps (s)')
plt.ylabel('Angle (rad)')
plt.title('Trajectòria del pèndol simple amb fricció')
plt.show()
```