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

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)