function r = isprime(number)
  p = factor(number);
  r = (size(p, '*') == 1);
endfunction

