% Skript som forsøker å finne høgast mogleg primtal

% Øvre grense
Nmax=1000;

% Initerar Nprimes - talet på primtal under Nmax
Nprimes=0;

for n=2:Nmax
  divisor=2;                    % Divisor er kandidat til faktor
  % Sjekkar om n er deleleg med alle tal <sqrt(n)
  while abs(floor(n/divisor)-n/divisor)>1e-13 & divisor<=floor(sqrt(n))
    divisor=divisor+1;
  end
  if divisor>sqrt(n)  
    disp([num2str(n), ' er eit primtal'])
    Nprimes=Nprimes+1;          % Oppdaterar talet på primtal
  end
end

% Bestemmer talet på primtal
disp(['Der er ',num2str(Nprimes),...
    ' primtal mellom 2 og ',num2str(Nmax),'.'])
