import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt

def model(y,t):
    k = 0.3
    dydt = -k * y
    return dydt

initial_conditions = [-10, -8, -6, -4, -2, 0, 2, 4, 6, 8, 10]

t = np.linspace(0,20)

for y0 in initial_conditions:
    y = odeint(model,y0,t)
    plt.plot(t, y, label=f'y(0) = {y0}')

plt.xlabel('time')
plt.ylabel('y(t)')
plt.grid()
plt.legend(loc='upper right')
plt.show()
