import gala.dynamics as gd

# when using dimensionless units, we don't need to specify units for the
# initial conditions
w0 = gd.PhaseSpacePosition(pos=[1.,0,0.25],
                           vel=[0.,0.3,0.])

# by default this uses Leapfrog integration
orbit = potential.integrate_orbit(w0, dt=0.1, n_steps=10000)

fig = orbit.plot(marker=',', linestyle='none', alpha=0.5)