function[ em ] = emirp( A ) index=1; for(i=1:A) if(isprime(i)) Q=str2num(fliplr(num2str(i))); if(isprime(Q)&&Q~=i) em(index)=i; index=index+1; end end end