diophant(8)

NAME
diophant - diophantine (Bezout) equation

CALLING SEQUENCE
[x,err]=diophant(p1p2,b)

PARAMETERS

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

EXAMPLE
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)