intone.m

改进等高线拟合代码

% INTONE  : Index determination for CUBVAL                           [f]

function i = intone(x,n,v,istart)

i=istart;
if  i >= n,
   i = 1;
elseif i< 1,
    i = 1;
end
if v < x(1) | v>x(n),
  error('v is buiten bereik')
end
if v < x(i),
  l = i;
  i = 1;
else
  if v <= x(i+1),
    return,
  else
    l = n;
  end
end

while l> i+1,
  k = floor((i + l)/ 2);
  if (v <= x(k)),
     l=k;
  else
     i=k;
  end
end  

return