DoubleMatrix2D V = alge.mult(B4, W); // V <-- B4*W
// T = B2.copy();
W = B4.copy();
W.assign (F.mult(cPade[4])); // W <-- W*cPade[4]
W.assign (T, F.plusMult(cPade[2])); // W <-- W + T*cPade[2]
W.assign (B0, F.plusMult(cPade[0])); // W <-- W + B0*cPade[0]
V.assign (W, F.plus); // V <-- V + W
W = V.copy();
W.assign(U, F.plus); // W = V + U, Padé numerator