HP Forums
HP Prime: CAS taylor - Printable Version

+- HP Forums (https://archived.hpcalc.org/museumforum)
+-- Forum: HP Museum Forums (https://archived.hpcalc.org/museumforum/forum-1.html)
+--- Forum: Old HP Forum Archives (https://archived.hpcalc.org/museumforum/forum-2.html)
+--- Thread: HP Prime: CAS taylor (/thread-258367.html)



HP Prime: CAS taylor - Alberto Candel - 12-10-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?


Re: HP Prime: CAS taylor - parisse - 12-11-2013

The third parameter is the order used for series expansion, not necessarily the final order.


Re: HP Prime: CAS taylor - Alberto Candel - 12-11-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)


Re: HP Prime: CAS taylor - parisse - 12-12-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...


Re: HP Prime: CAS taylor - Alberto Candel - 12-13-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"



Re: HP Prime: CAS taylor - Alberto Candel - 12-13-2013

But it looks like "truncate" is in Prime, and does more or less the same.