p, u, dx1 = self._init_poly_coefs(dx, dydx, n, p, D)