用Octave求泰勒展开(3)
广告
{{v.name}}
常见麦克劳林展开式(3)
\({\rm cos}(x)=x-\frac{x^2}{2!}+\cdots+\frac{(-1)^{n} }{(2n)!}x^{2n}+o(x^{2n})\)
求\({\rm cos}(x)\)的泰勒展开.
程序代码如下
function [text_result, numeric_result] = func41(order)
    pkg load symbolic;
    x = sym('x');
    question = cos(x);
    d = taylor(question, 'order', order + 1);
    text_result = ["\n", disp(d)];
    numeric_result = eval(d);
endfunction
计算5阶泰勒展开,结果如下
>> [text_result, numeric_result] = func41(5)
text_result =
     4    2
    x    x
    ── - ── + 1
    24   2

numeric_result = (sym)

     4    2
    x    x
    ── - ── + 1
    24   2
计算10阶泰勒展开,结果如下
>> [text_result, numeric_result] = func41(10)
text_result =
         10       8      6     4    2
        x        x      x     x    x
    - ─────── + ───── - ─── + ── - ── + 1
      3628800   40320   720   24   2

numeric_result = (sym)

         10       8      6     4    2
        x        x      x     x    x
    - ─────── + ───── - ─── + ── - ── + 1
      3628800   40320   720   24   2
友链