import astropy.coordinates as coord
coord.galactocentric_frame_defaults.set('v4.0')
pal5_c = coord.SkyCoord(ra=229.018*u.degree, dec=-0.124*u.degree,
                        distance=22.9*u.kpc,
                        pm_ra_cosdec=-2.296*u.mas/u.yr,
                        pm_dec=-2.257*u.mas/u.yr,
                        radial_velocity=-58.7*u.km/u.s)
rep = pal5_c.transform_to(coord.Galactocentric).data
pal5_w0 = gd.PhaseSpacePosition(rep)
pal5_mass = 2.5e4 * u.Msun
pal5_pot = gp.PlummerPotential(m=pal5_mass, b=4*u.pc, units=galactic)
mw = gp.MilkyWayPotential()
gen_pal5 = ms.MockStreamGenerator(df, mw, progenitor_potential=pal5_pot)
pal5_stream, _ = gen_pal5.run(pal5_w0, pal5_mass,
                              dt=-1 * u.Myr, n_steps=4000)
pal5_stream_c = pal5_stream.to_coord_frame(coord.ICRS)

fig, ax = plt.subplots(1, 1, figsize=(6, 4))
ax.scatter(pal5_stream_c.ra.degree, pal5_stream_c.dec.degree,
           marker='o', s=4, color='k', lw=0, alpha=0.1)
ax.set_xlim(250, 215)
ax.set_ylim(-10, 10)
ax.set_xlabel('RA [deg]')
ax.set_ylabel('Dec [deg]')
fig.tight_layout()