%  divisors2
%   divisors2(N) returns a vector containing 
%   the divisors of N  which are  ¾ sqrt(N)
%   See also PRIMES, ISPRIME, FACTOR
%   MB
function f = divisors2(n)
if prod(size(n))~=1, error('N must be a scalar.'); end
if (n < 0) | (floor(n) ~= n), error('N must be a positive integer.'); end
f=1; sqrtn=floor(sqrt(n));
if (n==1),  return; end
d=factor(n);  ld=length(d);
if ld==1,  f=[1 n]; return,  end
for  count=1:ld
	f1=f*d(count);
	f1(f1>sqrtn)='';
	f=unique([f f1]);
end
