function b=primzahl(m) n=floor(sqrt(m)); b=1; for i=2:n a=mod(m,i); if a==0 b=0; break end end end