public static (MathObject quotient, MathObject remainder) PolynomialDivision(MathObject u, MathObject v, MathObject x)
{
var q = (MathObject)0;
var r = u;
var m = r.DegreeGpe(new List <MathObject> {
x
});
var n = v.DegreeGpe(new List <MathObject> {
x
});
var lcv = v.LeadingCoefficientGpe(x);
while (m >= n && r != 0)
{
var lcr = r.LeadingCoefficientGpe(x);
var s = lcr / lcv;
q = q + s * (x ^ (m - n));
r = ((r - (lcr * (x ^ m))) - (v - lcv * (x ^ n)) * s * (x ^ (m - n))).AlgebraicExpand();
m = r.DegreeGpe(new List <MathObject> {
x
});
}
return(q, r);
}
public static MathObject LeadingCoefficientGpe(this MathObject u, MathObject x) =>
u.CoefficientGpe(x, u.DegreeGpe(new List <MathObject>()
{
x
}));