calculus - Using Taylor Polynomials Programmatically in Maple

I am trying to use a Taylor polynomial programmatically in Maple, but the following does not seem to work...

T[6]:=taylor(sin(x),x=Pi/4,6);convert(T[6], polynom, x);
  convert(T[6], polynom, x);
end proc;

All of the following also do not work:

  • f:=convert(T[6], polynom);
  • f:=convert(T[6], polynom, x);
  • f:=x->convert(T[6], polynom);
  • f:=x->convert(T[6], polynom, x);.

Is there a way of doing this without copying and pasting the output of convert into the definition of f?

3 Answers

  1. Albert- Reply


    It might also be natural to define T as a function.

  2. Andrew- Reply


    If I understood you correctly, this accomplishes what you want:

    f := proc(z)
        local p :: polynom;
        p := convert(T[6], polynom); 
        return subs(x = z, p)
    end proc
  3. Andy- Reply


    Several earlier answers involving procedures and subs will do the entire taylor series derivation, as well as the conversion to polynom, for each and every input. That is highly inefficient.

    You only need to produce the taylor result, and convert to polynom, once. With that result in hand you can then create an operator (with which to act on as many inputs as you wish, merely by evaluating the polynomial at the point but without having to recompute the whole taylor answer).

    Below is a way to create a procedure f with which to evaluate at any given point for the argument x. It computes the (truncated) taylor series and converts to polynom just once.

    > f:=unapply(convert(taylor(sin(x),x=Pi/4,6),polynom),x):

Leave a Reply

Your email address will not be published. Required fields are marked *

You can use these HTML tags and attributes <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>