def angles(nmodes,iseed,iy,iv):

    import numpy as np
    import math
    from randf_1 import randf_1

    ang=np.zeros(nmodes+1)
    teta=np.zeros(nmodes+1)

    fi,iseed,iy,iv=randf_1(nmodes,0.,2.*math.pi,iseed,iy,iv)
    psi,iseed,iy,iv=randf_1(nmodes,0.,2.*math.pi,iseed,iy,iv)
    alfa,iseed,iy,iv=randf_1(nmodes,0.,2.*math.pi,iseed,iy,iv)
    ang,iseed,iy,iv=randf_1(nmodes,0.,1.,iseed,iy,iv)
    for m in range(1,nmodes+1):
      teta[m]=np.arccos(1.-ang[m]/0.5) 

    return fi,psi,alfa,teta,iy,iv,iseed