How to get first n derivatives all at once efficiently?

[orebas]

Feature request: Based on the package design, it seems like it is probably possible to calculate the first n derivatives of a function at once (with one call).

Right now to get that out probably requires using internal details. It would be nice to have an API to return the value and the first n derivatives at once (in a tuple or vector). It would also be nice to have an example in the docs.