File:Greatestprimefactor,to1000.png

#!/usr/bin/python
import numpy as np
import matplotlib.pyplot as plt
import fractions
from mpmath import *
from pylab import rcParams
rcParams['figure.figsize'] = 10, 10

def gpf(n):
        i = 2
        while i * i <= n:
                if n % i:
                        i += 1
                else:
                        n //= i
        return n

x=np.arange(0,1000,1);
f=np.vectorize(gpf)
y=f(x)

fig, ax = plt.subplots()

plt.xlabel(r'$n$')
plt.ylabel(r'$\mathrm{gpf}(n)$')

plt.scatter(x,y,linewidth=0.01,color='Black')
plt.xlim(0,1000)
plt.ylim(0,1000)
plt.savefig('greatestprimefactor,to1000.png',bbox_inches='tight',pad_inches=0.15)