%function [pinds,FV] = convhullinds(P)
function pinds = convhullinds(P)
% pinds are the corners of conv(P)

pinds = [];
FV=[];

for i=1:size(P,2)
  Pi = [P(:,1:i-1) P(:,i+1:end)];
  q = P(:,i);
  %[b,fv] = inconvhull(q,Pi);
  [b,x] = inconvhull(q,Pi);
  if ~b
    pinds(end+1) = i;
  end
  %FV(end+1)=fv;
end

return