#!/usr/bin/python3

from galpy.potential import (
    NFWPotential,
    HernquistPotential,
    MiyamotoNagaiPotential
)

mp= MiyamotoNagaiPotential(a=0.5,b=0.0375,normalize=.6)
np= NFWPotential(a=4.5,normalize=.35)
hp= HernquistPotential(a=0.6/8,normalize=0.05)
