function [nber_primes,vec_primes ] = allprimes(ub)
%Provides the number of primes and a list of the primes that are <=ub
nber_primes=0;
vec_primes=[];
for i=2:ub
    if isprime_myfun(i)==1
        nber_primes=nber_primes+1;
        vec_primes=[vec_primes; i];
    end
end

