Anonymous user
Polynomial long division: Difference between revisions
→{{header|D}}
Line 290:
poly.popBack();
return (cast(int)poly.length) - 1;
}
Line 302 ⟶ 298:
throw new Exception("ZeroDivisionError");
if (dN >= dD) {
q =
while (dN >= dD) {
double mult = q[dN - dD] = N[$-1] / d[$-1];
d[] *= mult;
Line 312 ⟶ 308:
r = N;
} else {
q = [0.0];
r = N;
}
Line 322 ⟶ 318:
double[] D = [-3.0, 1.0, 0.0, 0.0];
auto qr = polyDiv(N, D);
writefln(" %s / %s = %s remainder %s", N, D, qr._0, qr._1);
}</lang>
| |||