function [I,II]=findI(A,b,x);

eps=10^(-6);

[m,n]=size(A);
II=zeros(1,m);
I=[];
for i=1:m
  if abs(b(i)-A(i,:)*x')<eps
   II(i)=1;
   I=[I i];
  end;
end;

return;


