;[mi,sigma,solution,minp,topp,h.ellipse.N,h.ellipse.t]=...</span><br>
<span class=comment>%%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;mmln(sets.X,epsilon,iter,h.ellipse.t,h.ellipse.N);</span><br>
<br>
&nbsp;&nbsp;&nbsp;options.tmax&nbsp;=&nbsp;iter+h.ellipse.t;<br>
&nbsp;&nbsp;&nbsp;options.eps&nbsp;=&nbsp;epsilon;<br>
&nbsp;&nbsp;&nbsp;<span class=keyword>if</span>&nbsp;h.ellipse.t&nbsp;==&nbsp;0,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;model=mmgauss(sets.X,options);<br>
&nbsp;&nbsp;&nbsp;<span class=keyword>else</spa