Taylor series error estimate built-in?

Matlab Symbolics has the Function Taylor to derive the series expansion, to a given order, and about a given point. But amazingly, the function does not generate the corresponding error estimate. The user can devise a work-around by getting the series for the next highest order, and evaluating the n+1th derivative . . . . It would be better to have the function estimate the error automatically.