Examples
s = poly(0, 's'); p1 = (s+3)^2; p2 = (1+s); x1 = s; x2 = (2+s); [x,err] = diophant([p1 p2], p1*x1 + p2*x2); p1*x1 + p2*x2 - p1*x(1) - p2*x(2) |  |  |
diophantine (Bezout) equation
[x,err]=diophant(p1p2,b)
polynomial vector p1p2 = [p1 p2]
polynomial
polynomial vector [x1; x2]
diophant solves the bezout equation:
p1*x1+p2*x2=b with p1p2 a polynomial vector.
If the equation is not solvable, err = ||p1*x1 + p2*x2 - b|| / ||b||
else err = 0.
s = poly(0, 's'); p1 = (s+3)^2; p2 = (1+s); x1 = s; x2 = (2+s); [x,err] = diophant([p1 p2], p1*x1 + p2*x2); p1*x1 + p2*x2 - p1*x(1) - p2*x(2) | | |