bpoly = BPoly.from_derivatives(x, yyp, orders=3)