File:Mobiusplot,on0to100.png

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

def numberOfFactors(n):
        counter = 0
        while (n > 1):
                i=2
                while (i <= n):
                        if(n%i ==0):
                                n /= i
                                counter += 1
                                break
                        i += 1
        return counter

def isDivisibleBySquare(n):
        i=2
        while(i**2 <= n):
                if (n%(i**2)==0):
                        return 1
                i += 1
        return 0

def isValid(n):
        if (n <= 0):
                return 0
        elif (n-int(n) != 0):
                return 0
        return 1

def mobius(n):
        if not isValid(n):
                return 0
        elif (n==1):
                return 1
        elif(isDivisibleBySquare(n)):
                return 0
        else:
                return (-1)**numberOfFactors(n)
fig, ax = plt.subplots()

x=np.arange(0,100,1)
f=np.vectorize(mobius)
y=f(x)

plt.xlabel(r'$n$')
plt.ylabel(r'$\mu(n)$')