Pythagorean Triple Function

pythagoreantriple.py

def ppt(limit):
    triples = []
    for m in range(2,int(math.sqrt(limit))+1):
        for n in range(1,m):
            if (m+n) % 2 == 1 and math.gcd(m,n) == 1:
                a = m**2 + n**2
                b = m**2 - n**2
                c = 2*m*n
             
                p = max(a,b,c)
                
                for k in range(1,int(limit/p)+1):
                    triples.append([k*b,k*c,k*a])
    return triples