#============================================================================== import numpy as np import matplotlib import matplotlib.pyplot as plt #============================================================================== # Programa mostrando modelo discreto de crescimento logístico # de duas populações com taxas de crescimento distintas # -------------------------------------------------------- n=121 # número de passos lamb=0.1 # taxa média mensal de crescimento da população (coelhos) alfa=0.08 # taxa média mensal de crescimento da população (lebres) x=np.zeros(n) y=np.zeros(n) x[0] = 30 # população inicial de coelhos y[0] = 30 # população inicial de lebres mx = 2000 # saturação de coelhos my = 2500 # saturação de lebres # # segue a integração do modelo pelo número de passos for k in range(0,n-1): x[k+1]=x[k]+lamb*x[k]*(1-x[k]/mx) y[k+1]=y[k]+alfa*y[k]*(1-y[k]/my) # # Alguns gráficos da evolução das populações # ------------------------------------------ plt.axis([0,40,0,1000]) plt.plot(x,'r-') plt.plot(y,'g-') plt.title('Evolução da população no tempo') plt.ylabel('número de animais') plt.xlabel('meses') plt.plot([5,8],[500,500],'r-') plt.plot([5,8],[400,400],'g-') plt.text(9,500,r'coelhos',color='r') plt.text(9,400,r'lebres',color='g') plt.savefig('coelhos_e_lebres_30.png') plt.show() plt.axis([0,n-1,0,3000]) plt.title('Evolução da população no tempo') plt.ylabel('número de animais') plt.xlabel('meses') plt.plot(x,'r-') plt.plot(y,'g-') plt.plot([5,8],[2000,2000],'r-') plt.plot([5,8],[1500,1500],'g-') plt.text(9,2000,r'coelhos',color='r') plt.text(9,1500,r'lebres',color='g') plt.savefig('coelhos_e_lebres_60.png') plt.show() #