function [x,y,iter]=bsegno(a,b,c,d,numpoints,maxiter,epsilon)

x=linspace(a,b,numpoints);
y=linspace(c,d,numpoints);
nx=length(x);
ny=length(y);
 for i=1:nx
  for j=1:ny
	z=x(i)+y(j)*1i;
	[s,n]=segno(z,maxiter,epsilon);
	iter(j,i)=n;
  end
 end
end
