gin&lt;4, t=reshape(1:nn,size(x)); <span class="keyword">end</span>
0053 <span class="keyword">if</span> nargin&lt;3, n=0; <span class="keyword">end</span>
0054 vor=n;nach=n;
0055 <span class="keyword">if</span> n==0,
0056     ys=(y(3:nn)-y(1:nn-2))./(x(3:nn)-x(1:nn-2));
0057     dydx=diff(y)./diff(x);
0058     yss=diff(dydx)./(x(3:nn)-x(1:nn-2))*2;
0059     curv=yss./((1+ys.^2).^1.5);
0060     curv=[curv(1);curv(:);curv(length(curv))];
0061 <span class="keyword">else</span>
0062   <span class="keyword">if<