Finding Taylor Coefficients Using Octave

[Gregg Anderson]

Sirs,
In MATLAB one can find the 15th-order Taylor polynomial of f(x) = x^3*tan(x^2) (assuming you have the symbolic module installed) by typing at the command line the following:

       >> taylor(x^3*tan(x^2),16)       ->ans = x^5+1/3*x^9+2/15*x^13

MATLAB is using n! = gamma(n+1).

How can I do this using Octave?  I'm aware of Octave packages, having installed the following:
1. gsl
2. integration
3. odepkg
4. plot
5. specfun
6. symbolic
Also, any documentation on how to actually use these packages?
Thanks.
Gregg  Anderson