from random import gauss

# gauss(mu, sigma)

def angles_M1(n, alpha, gamma, sigma):
    A = []
    A.append( gauss(0, (alpha**2.0)/16.0) )
    for i in range(1,n+1):
        A.append( gauss(-(gamma**4.0)*A[i-1]/2.0, (alpha**2.0)/16.0 + (sigma**2.0)*A[i-1]/alpha) )
    return A

def distantes_M1(n, delta):
    A = []
    for i in range(n):
        A.append( gauss(delta, (delta**2.0)/64) )
    return A

if __name__=="__main__":
    A = angles_M1(15, 22.5, 0.04, 0.2)
    for a in A:
        print(a)