diophant

Arguments

p1p2

polynomial vector p1p2 = [p1 p2]

b

polynomial

x

polynomial vector [x1; x2]

Description

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.

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)