import numpy as np

def divisores(num):
    n = np.arange(1,num)
    d = num % n
    zeros = d == 0
    return n[zeros]

def divisores_comuns(num1, num2):
    divs1 = set(divisores(num1))
    divs2 = set(divisores(num2))
    return divs1 & divs2


print(divisores_comuns(15,20))