def plotprofMulti(self, ini, end, delta, what_specie, xlim1, xlim2,<EOL>ylim1, ylim2, symbol=None):
plotType=self._classTest()<EOL>if plotType=='<STR_LIT>':<EOL><INDENT>for i in range(ini,end+<NUM_LIT:1>,delta):<EOL><INDENT>step = int(i)<EOL>if symbol==None:<EOL><INDENT>symbol_dummy = '<STR_LIT:->'<EOL>for j in range(len(what_specie)):<EOL><INDENT>self.plot_prof_1(step,what_specie[j],xlim1,xlim2,ylim1,ylim2,symbol_dummy)<EOL><DEDENT><DEDENT>else:<EOL><INDENT>for j in range(len(what_specie)):<EOL><INDENT>symbol_dummy = symbol[j]<EOL>self.plot_prof_1(step,what_specie[j],xlim1,xlim2,ylim1,ylim2,symbol_dummy)<EOL><DEDENT><DEDENT>filename = str('<STR_LIT>' % step)+'<STR_LIT>'<EOL>pl.savefig(filename, dpi=<NUM_LIT>)<EOL>print('<STR_LIT>', filename)<EOL>pl.clf()<EOL><DEDENT><DEDENT>else:<EOL><INDENT>print('<STR_LIT>'+str(self.__class__))<EOL>return<EOL><DEDENT>
create a movie with mass fractions vs mass coordinate between
xlim1 and xlim2, ylim1 and ylim2. Only works with instances of
se.

Parameters
----------
ini : integer
    Initial model i.e. cycle.
end : integer
    Final model i.e. cycle.
delta : integer
    Sparsity factor of the frames.
what_specie : list
    Array with species in the plot.
xlim1, xlim2 : integer or float
    Mass coordinate range.
ylim1, ylim2 : integer or float
    Mass fraction coordinate range.
symbol : list, optional
    Array indicating which symbol you want to use. Must be of
    the same len of what_specie array.  The default is None.