prime =[]; for i=500:599 isDivided = false; for j=2:i-1 if rem(i,j) == 0 isDivided = true; end end if isDivided == false prime = [prime, i]; end end