"comment">% identifiy species index</span>
[labels, ~, inds]= unique(str.species);

colors = [1 0 0; 0 0.8 0; 0 0 1];

<span class="keyword">for</span> i = 1:3
    pts_i = pts(inds == i, :);
    drawPoint(pts_i, <span class="string">'marker'</span>, <span class="string">'x'</span>, <span class="string">'color'</span>, colors(i,:), <span class="string">'linewidth'</span>, 2);

    elli = equivalentEllipse(pts_i);
    drawEllipse(elli, <span class="string">'color'</span>, colors(i,:), <span class="string">'li