Olá
Boa noite
Estou tentando criar um código que me retorne os divisores comuns entre 2 variáveis que recebem números inteiros. E depois mostre os divisores comuns.

Para me retornar os divisores estou usando esse código para teste:

import numpy as np

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

divisores(30) 


Porém eu preciso que essa função me retorne os divisores de 2 valores e não apenas um. Além de comparar os divisores comuns entre eles, pois não consegui pensar em uma forma de fazer. =)