Added new anisotropic density function to allow for smooth truncation (i.e. n = Real), velocity is still done numerically, but the code can be used to do fast fitting with potonly=true