程序与数学：应用泰勒展开式计算自然对数

import math
  
#计算泰勒展开式的最大前n项
max_n = 100000
  
# 计算ln(x)函数的泰勒展开式第n项的值
def taylor(n,x):
    value = (1/n)*(((x-1)/(x+1))**n)
    return value
     
  
# 返回ln(x)泰勒展开式当前项数
def get_n():
    n=1
    while True:
        # 若当前项数大于max_n，计算结束
        if n > max_n:
            break
        yield n
        n=n+2
         
# 计算ln(x)函数
def ln(x):
    sum = 0
    for n in get_n():
        value = taylor(n,x)
        sum += value
        # 若当前项的值<1e-10,计算结束
        if value <= 1e-10:
            break
    return 2 * sum
     
  
# 程序入口
if __name__ == '__main__':
  
    x = float(input("请输入一个正实数:\n"))
    print("ln函数计算值：%.3f" % ln(x))
    print("math.log函数计算值：%.3f" % math.log(x))

%{
绘图：
（1）绘制ln(x)在区间[1,10]内的曲线
（2）绘制ln(x)在x=3邻域内的泰勒展开式
目的：观察泰勒展开式对函数的拟合
%}
  
% 绘制ln(x)函数曲线
% 区间[1,10]内创建100个点
x1 = linspace(1,10,100);
% 计算f(x)=log(x)函数的y坐标
y1 = log(x1);
% 绘制ln(x)函数曲线
plot(x1,y1)
hold on
  
% 绘制ln(x)在x=3邻域内的泰勒展开式曲线
% 定义符号变量x,y,f
syms x y f;
% 定义自然对数函数
% ln(x)在matlab为log(x)
y = log(x);
% 应用matlab的taylor函数获取log(3)泰勒展开式前10项
f = taylor(y,x,3,'order',10);
% log(3)邻域内内创建100个点
t = linspace(2,4,100);
% 使用t替换符号变量
f1=subs(f,'x',t);
% 绘制泰勒展开式图像
plot(t,f1,'.-r')