Posts: 72
Threads: 16
Joined: Sep 2013
Hi,
In CAS, for example, if I input
taylor(LN(1+x),x,3)
the output is
x -(1/2)*x^2+(1/3)*x^3+x^4*order_size(x)
as expected from the documentation. But
taylor((x+x^2+2*x^3)/(1+2*x-x^2+6*x^3), x,3)
outputs
x-x^2+5*x^3-17*x^4+x^5*order_size(x)
instead of
x-x^2+5*x^3+x^4*order_size(x)
Why the extra term?
Posts: 38
Threads: 0
Joined: Aug 2013
The third parameter is the order used for series expansion, not necessarily the final order.
Posts: 72
Threads: 16
Joined: Sep 2013
Thank you. So if I want to obtain the Taylor polynomial of degree n of a function f(x), what is the correct command to use? (I want the result to be a polynomial of degree less than or equal to n, and I do not want the extra "order_size" term)
Posts: 38
Threads: 0
Joined: Aug 2013
rem(convert(series(expression,x=limit_point,order),polynom),x^(order+1)) should do that: convert(.,polynom) will remove the order_size remainder term, and rem will remove all monomials of degree>=order+1. You have no warranty to get the full Taylor expansion however, since the remainder term has been erased...
Posts: 72
Threads: 16
Joined: Sep 2013
Thanks. But it looks that such command is not available in the HP prime. The example in page 157 of the "Symbolic algebra and Mathematics with Xcas"
convert(series(sin(x),x=0,6),polynom)
outputs
"Error: Unmatch control word"
Posts: 72
Threads: 16
Joined: Sep 2013
But it looks like "truncate" is in Prime, and does more or less the same.