sp;&nbsp;&nbsp;&nbsp;&nbsp;model=mmgauss(sets.X,options,init_model);<br>
&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>
<br>
&nbsp;&nbsp;&nbsp;mi=&nbsp;model.Mean;<br>
&nbsp;&nbsp;&nbsp;sigma&nbsp;=&nbsp;model.Cov;<br>
&nbsp;&nbsp;&nbsp;solution&nbsp;=&nbsp;model.exitflag;<br>
&nbsp;&nbsp;&nbsp;minp&nbsp;=&nbsp;model.lower_bound(<span class=keyword>end</span>);<br>
&nbsp;&nbsp;&nbsp;topp&nbsp;=&nbsp;model.upper_bound(<span class=keyword>end</span>);<br>
&nbsp;&nbsp;&nbsp;h.ellipse.N&nbsp;=&nbsp;model.