// build the P matrix elements from Taylor series formulas
final BigFraction[] pI = pData[i];
final int factor = -(i + 1);
int aj = factor;
for (int j = 0; j < pI.length; ++j) {
pI[j] = new BigFraction(aj * (j + 2));
aj *= factor;
}
}
return new Array2DRowFieldMatrix<BigFraction>(pData, false);