import matplotlib.pyplot as plt

def pop_cigognes(a,n):
    u=0.3
    for i in range (n):
        u=a*u*(1-u)
    return(u)

for a in range(1,5):
    x=[n for n in range(50)]
    y=[pop_cigognes(a,n) for n in range(50)]
    plt.scatter(x,y)
    plt.show()
    plt.figure()

print(round(pop_cigognes(1,5),3))
print(round(pop_cigognes(3,10),3))

print(round(pop_cigognes(1,5),3))
print(round(pop_cigognes(3,5),3))

print(round(pop_cigognes(1,8),3))
print(round(pop_cigognes(3,8),3))

print(round(pop_cigognes(1,10),3))
print(round(pop_cigognes(3,10),3))

print(round(pop_cigognes(1,15),3))
print(round(pop_cigognes(3,15),3))