Changes from BHU workshop . changed variable names in list intro . pendulum.pychanged to pendulum_plot.py . Also an extra space added on slide 33 session1 day1
import math
def aliquot(n):
sum = 1
i = 2
while i * i < n:
if n % i == 0:
sum += i + (n / i)
i += 1
if i*i == n: sum += i
return sum
amicable = []
for n in range(1000, 10000):
m = aliquot(n)
if m > n and aliquot(m) == n:
amicable.append((m, n))
print amicable
# please please please profile this.