8   cc=<a href="curvature.html" class="code" title="function curv=curvature(x,y,n,t)">curvature</a>(rho,eta,n,ak);
0019 <span class="keyword">end</span>
0020 mm=min(find(cc==max(cc)));
0021 <span class="keyword">if</span> (mm==1)|(mm==length(cc)-1),
0022     mm=1;flag=0;
0023     <span class="keyword">while</span>((flag==0)&amp;(mm+1&lt;length(rho))),
0024         mm=mm+1;
0025         <span class="keyword">if</span>((cc(mm)&gt;cc(mm-1))&amp;(cc(mm)&gt;cc(mm+1))) flag=1; <span class="keyword">end</span>
002