def make(N_part=100, rmax=30, rmin=2, omega=10, seed=42):
    random = np.random.RandomState(seed=seed)
    r = (((rmax - rmin) * random.random_sample(size=N_part)) + rmin)
    phi0 = ((2 * np.pi) * random.random_sample(size=N_part))
    mass = (100000000.0 * np.ones_like(r))
    x = (r * np.cos(phi0))
    y = (r * np.sin(phi0))
    z = ((1 * random.random_sample(size=N_part)) - 0.5)
    xdot = ((((- 1) * omega) * r) * np.sin(phi0))
    ydot = ((omega * r) * np.cos(phi0))
    zdot = np.zeros_like(xdot)
    pos = np.array([x, y, z]).T
    vel = np.array([xdot, ydot, zdot]).T
    return (mass, pos, vel)
def make(N_part=100, rmax=30, rmin=2, omega=10, seed=42):
    random = np.random.RandomState(seed=seed)
    r = (((rmax - rmin) * random.random_sample(size=N_part)) + rmin)
    phi0 = ((2 * np.pi) * random.random_sample(size=N_part))
    mass = (100000000.0 * np.ones_like(r))
    x = (r * np.cos(phi0))
    y = (r * np.sin(phi0))
    z = ((1 * random.random_sample(size=N_part)) - 0.5)
    xdot = ((((- 1) * omega) * r) * np.sin(phi0))
    ydot = ((omega * r) * np.cos(phi0))
    zdot = np.zeros_like(xdot)
    pos = np.array([x, y, z]).T
    vel = np.array([xdot, ydot, zdot]).T
    return (mass, pos, vel)