Dfp h = step.divide(field.newDfp(100000));
for (int j = 0; j < 20; ++j) {
Dfp x = field.newDfp(j).multiply(step);
Dfp[] y = interpolator.value(x);
Dfp[] yP = interpolator.value(x.add(h));
Dfp[] yM = interpolator.value(x.subtract(h));
Assert.assertEquals(p.length, y.length);
for (int k = 0; k < p.length; ++k) {
Assert.assertEquals(p[k].value(x.getReal()),
y[k].getReal(),
1.0e-8 * FastMath.abs(p[k].value(x.getReal())));