def density_profile(self,ixaxis='<STR_LIT>',ifig=None,colour=None,label=None,fname=None):
pT=self._classTest()<EOL>if pT is '<STR_LIT>':<EOL><INDENT>x = self.get(ixaxis)<EOL>if ixaxis is '<STR_LIT>':<EOL><INDENT>x = x*ast.rsun_cm<EOL><DEDENT>y = self.get('<STR_LIT>')<EOL><DEDENT>elif pT is '<STR_LIT>':<EOL><INDENT>if fname is None:<EOL><INDENT>raise IOError("<STR_LIT>")<EOL><DEDENT>x = self.se.get(fname,ixaxis)<EOL>y = np.log10(self.se.get(fname,'<STR_LIT>'))<EOL><DEDENT>else:<EOL><INDENT>raise IOError("<STR_LIT>")<EOL><DEDENT>if ixaxis is '<STR_LIT>':<EOL><INDENT>x = np.log10(x)<EOL>xlab='<STR_LIT>'<EOL><DEDENT>else:<EOL><INDENT>xlab='<STR_LIT>'<EOL><DEDENT>if ifig is not None:<EOL><INDENT>pl.figure(ifig)<EOL><DEDENT>if label is not None:<EOL><INDENT>if colour is not None:<EOL><INDENT>pl.plot(x,y,color=colour,label=label)<EOL><DEDENT>else:<EOL><INDENT>pl.plot(x,y,label=label)<EOL><DEDENT>pl.legend(loc='<STR_LIT>').draw_frame(False)<EOL><DEDENT>else:<EOL><INDENT>if colour is not None:<EOL><INDENT>pl.plot(x,y,color=colour)<EOL><DEDENT>else:<EOL><INDENT>pl.plot(x,y)<EOL><DEDENT><DEDENT>pl.xlabel(xlab)<EOL>pl.ylabel('<STR_LIT>')<EOL>
Plot density as a function of either mass coordiate or radius.

Parameters
----------
ixaxis : string
    'mass' or 'radius'
    The default value is 'mass'
ifig : integer or string
    The figure label
    The default value is None
colour : string
    What colour the line should be
    The default value is None
label : string
    Label for the line
    The default value is None
fname : integer
    What cycle to plot from (if SE output)
    The default value is None