Matlab 数值计算泰勒公式

      1.利用n阶泰勒多项式展开来估算函数f(x)=e^x在指定点x处的近似值。要求绝对误差在最大阶数MAXN以内达到给定精度EPS。

    taylorcal函数：

function taylorcal()
eps=0.00001;
MAXN=20;
x=input('存储输入的浮点数/向量');
f1=zeros(1,length(x));
for j=1:length(x)
     xzs=1;
     xjc=1;
     xsum=1;
     accexp=exp(x(j));
for i=1:MAXN
    xzs=xzs.*x(j);
    xjc=xjc.*i;
    xsum=xsum+xzs./xjc;
    if abs(xsum-accexp)<eps
        f1(j)=xsum;
        break
    end
    f1(j)=-1;
end
end
disp([x' f1'])
end

截图：

2.测试程序，分别查看当x=-6,5时的泰勒展开的前20项的值。（泰勒公式过程分析）

  taylorcaltest函数

function taylorcaltest()
eps=0.00001;
MAXN=20;
x=input('存储输入的浮点数/向量');
flp=zeros(1,length(x));
for j=1:length(x)
    xzs=1;
    xjc=1;
    xsum=1;
    accexp=exp(x(j));
 for i=1:MAXN
     xzs=xzs.*x(j);
     xjc=xjc.*i;
     term(i)=xzs./xjc;
     xsum=xsum+xzs./xjc;
     flp(i)=xsum;
 if abs(xsum-accexp)<eps
       fl(j)=xsum;
  break
 end
      fl(j)=-1;
 end
end
disp([[1:20],flp,term])
end
         

截图：