#!/usr/bin/python

from pylab import *

def gcd(a,b):
    r = a%b
    while r!=0:
        a,b,r = b,r,b%r
    return b

def Coprime(x,y):
    return gcd(x,y)==1

X,Y = [],[]
for i in range(1,51):
    for j in range(1,51):
        if Coprime(i,j):
            X.append(i)
            Y.append(j)
plot(X,Y,'b.')
show()