code: def plotDensity(self,t=0.,
                      rmin=0.,rmax=1.5,nrs=21,zmin=-0.5,zmax=0.5,nzs=21,
                      phi=None,xy=False,
                      ncontours=21,savefilename=None,aspect=None,log=False,
                      justcontours=False):
          
          return plotDensities(self,rmin=rmin,rmax=rmax,nrs=nrs,
                               zmin=zmin,zmax=zmax,nzs=nzs,phi=phi,xy=xy,t=t,
                               ncontours=ncontours,savefilename=savefilename,
                               justcontours=justcontours,
                               aspect=aspect,log=log)
docs: NAME:
  
             plotDensity
  
          PURPOSE:
  
             plot the density of this potential
  
          INPUT:
  
             t= time to plot potential at
  
             rmin= minimum R (can be Quantity) [xmin if xy]
  
             rmax= maximum R (can be Quantity) [ymax if xy]
  
             nrs= grid in R
  
             zmin= minimum z (can be Quantity) [ymin if xy]
  
             zmax= maximum z (can be Quantity) [ymax if xy]
  
             nzs= grid in z
  
             phi= (None) azimuth to use for non-axisymmetric potentials
  
             xy= (False) if True, plot the density in X-Y
  
             ncontours= number of contours
  
             justcontours= (False) if True, just plot contours
  
             savefilename= save to or restore from this savefile (pickle)
  
             log= if True, plot the log density
  
          OUTPUT:
  
             plot to output device
  
          HISTORY:
  
             2014-01-05 - Written - Bovy (IAS)