matlab关于函数的有关问题

matlab关于函数的问题
底下的代码关于求矩阵中素数的为什么不能运行？求指点。。
function [res fctr]= is_prime(n)
res = ones(size(n));
fctr = ones(size(n));
for K=1:length(n(:))
  cur_n = n(K);
  [res(K) fctr(K)] = is_a_prime(cur_n);
end
function [res fctr] = is_a_prime(n)
res = true;
for I=2:floor(sqrt(n))
  if mod(n, I) == 0
  fctr = I; % 找到一个因子
  res = false;
  return;
  end
end
A=[2:100]
[a,b]=is_prime(A)


[解决办法]
错误提示是什么呢，麻烦贴上来。而且函数写的貌似有点问题。
给一个自带的判断素数的。

C/C++ code

isprime(V):判断V中元素是否为质数，是返回1，不是返回0比如“ v=[1 2 4 5];isprime(v)