function v = prim_e (a)
v=1;
if (a<2) v=1;
endif;
for i = 2:sqrt(a)
if (a==floor(a/i)*i) v=0;
 endif
endfor

endfunction
