from matplotlib import pyplot as plt
from math import *
import time
def suite_euler(n,N):       # à définir
    a=1
    for i in range (n):
        a=a+(a**2)*(1/(3*N))
    return(a)


def liste_euler(a,H):
    F=[]
    for i in range(a):
        F.append(suite_euler(i,H))
    return(F)

def liste_euler_abs(a,H):
    K=[]
    for i in range(a):
        K.append(i/H)
    return(K)

K=[0.01*i for i in range(300)]                      # ces trois ligne servent a comparer avec la vrai solution
C=[-3/(i*0.01-3) for  i in range(300)]
plt.plot(K,C)
a=time.clock()
B=liste_euler(1650,550)
print(time.clock()-a)
plt.plot(liste_euler_abs(1650,550),B)   #en abcisse on va jusqua a/H
plt.show()