function p=prime(n)
  s=0
  for i=2:sqrt(n)
    if rem(n,i)==0
      s=s+1
      break
    endif
  endfor
  if s>0
    p='not prime'
    else p='prime'
  endif
endfunction