Polynomial long division: Difference between revisions

Line 306:
double[] d = mul(0.0, dN - dD) ~ D;
double mult = q[dN - dD] = N[$-1] / d[$-1];
foreachd[] (ref*= elmult; d)
N[] el *-= multd[];
double[] newN;
foreach (i; 0 .. min(N.length, d.length))
newN ~= N[i] - d[i];
N = newN;
dN = degree(N);
}
Anonymous user