Script para modelo de crecimiento poblacional logistico en plantas

 import matplotlib.pyplot as plt

def logistic_growth(K, r, N0, timesteps):

    """

    Modelo de crecimiento logístico para una población de plantas.

    

    Parámetros:

    - K: Capacidad de carga del ambiente

    - r: Tasa de crecimiento intrínseco

    - N0: Población inicial

    - timesteps: Número de pasos de tiempo para la simulación

    

    Retorna:

    - Una lista con la población de plantas en cada paso de tiempo

    """

    N = N0  # Población inicial

    population = [N0]  # Lista para almacenar la población en cada paso de tiempo

    

    for _ in range(1, timesteps):

        delta_N = r * N * (1 - N / K)  # Cambio en la población

        N += delta_N  # Actualizar la población

        population.append(N)  # Almacenar la nueva población

    

    return population

# Parámetros del modelo

K = 1000  # Capacidad de carga del ambiente

r = 0.1  # Tasa de crecimiento intrínseco

N0 = 50  # Población inicial

timesteps = 100  # Número de pasos de tiempo

# Simulación del crecimiento de la población

population = logistic_growth(K, r, N0, timesteps)

# Visualización de los resultados

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

plt.plot(population, label='Población de plantas')

plt.axhline(y=K, color='r', linestyle='--', label='Capacidad de carga (K)')

plt.xlabel('Tiempo')

plt.ylabel('Población')

plt.title('Crecimiento de una Población de Plantas')

plt.legend()

plt.show()