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

numeric_result = (sym)

     5    4    3    2
    x    x    x    x
    ── + ── + ── + ── + x
    5    4    3    2
计算10阶泰勒展开，结果如下
>> [text_result, numeric_result] = func45(10)
text_result =
     10    9    8    7    6    5    4    3    2
    x     x    x    x    x    x    x    x    x
    ─── + ── + ── + ── + ── + ── + ── + ── + ── + x
    10    9    8    7    6    5    4    3    2

numeric_result = (sym)

     10    9    8    7    6    5    4    3    2
    x     x    x    x    x    x    x    x    x
    ─── + ── + ── + ── + ── + ── + ── + ── + ── + x
    10    9    8    7    6    5    4    3    2
友链