function [ y ] = isprime_myfun( a )
%The function takes as input a number and outputs 0 or 1 depending on if it
%is prime or not
y=1;
for i=2:ceil(sqrt(a))
    if mod(a,i)==0
        y=0;
        break
    end
end
if a==2
    y=1;
end

