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

numeric_result = (sym)

     5     4    3    2
    x     x    x    x
    ─── + ── + ── + ── + x + 1
    120   24   6    2
计算10阶泰勒展开，结果如下
>> [text_result, numeric_result] = func39(10)
text_result =
       10        9       8       7     6     5     4    3    2
      x         x       x       x     x     x     x    x    x
    ─────── + ────── + ───── + ──── + ─── + ─── + ── + ── + ── + x + 1
    3628800   362880   40320   5040   720   120   24   6    2

numeric_result = (sym)

       10        9       8       7     6     5     4    3    2
      x         x       x       x     x     x     x    x    x
    ─────── + ────── + ───── + ──── + ─── + ─── + ── + ── + ── + x + 1
    3628800   362880   40320   5040   720   120   24   6    2
友链