НОД двух многочленов

[KonraDgeo]

Доброе время суток
Столкнулся с проблемой во время написания функции по нахождению НОД двух многочленов. Мой вариант ниже:

	def gcd(self,pol):
		spol = self
		while spol.deg>1:
			spol,pol=pol,spol%pol
		return spol

Подразумевается, что spol - больший (по степени) многочлен
Вот как перегружен остаток от деления:

def __mod__(self, pol):
		new_list = list(self.coeffs)
		mod = Polynom(new_list)
		div = Polynom([0.] * (self.deg - pol.deg + 1))

		while mod.deg >= pol.deg:
			k = mod.deg - pol.deg
			a = mod.coeffs.pop()
			mod.deg -= 1
			div.coeffs[k] = a / pol.coeffs[-1]

			for i in range(k, mod.deg + 1):
				mod.coeffs[i] -= a * pol.coeffs[i - k]
			mod.remove_trailing_zeros()
		return mod

Заранее спасибо

[Alex11223]

Так проблема-то в чем?

[KonraDgeo]

Считает неправильно. Не могу понять, где ошибка

[KonraDgeo]

Все, я во всем наконец-то разобрался
Если кому интересно, внизу решение:

def __mod__(self, pol):
		mod = Polynom(list(self.coeffs))
		div = Polynom([0.]*(self.deg-pol.deg+1))
		
		while mod.deg>=pol.deg:
			k = mod.deg-pol.deg
			a = mod.coeffs[-1]/pol.coeffs[-1]
			div.coeffs[-k]=a
	
			for i in range(pol.deg,-1,-1):
				mod.coeffs[i+k]-=a*pol.coeffs[i]
				
			mod.remove_trailing_zeros()
		return mod

def gcd(self,pol):
		spol = self
		while pol.deg!=0:
			buff = pol
			spol,pol=pol,spol%pol
			if pol.coeffs == []:
				pol = buff
				break
		return pol