matlab - 我如何在不影响我的（xa）术语的情况下获得“a”处的函数的导数任何建议？

Question

我需要这个函数来输出一个以 a 为中心的 P n 次泰勒多项式，这是我拥有的代码，它工作得几乎完美，但我不确定如何在不影响我的 (xa) 项的情况下在“a”处获得函数的导数...从我看到的 diff(f,k) 计算 f 的第 k 个导数，但不能插入 a。看起来我的代码将需要另一个 matlab 函数来完成任何建议？

function [P] = mytaylor(f,a,n) 
  f = sym(f);
  syms x;
  terms = 0;
  for k = 0:n
    fk = diff(f,k);
    terms = terms + fk*(x-a)^(k)/factorial(k);   
  end
  P = terms;
end

Answer 1

您需要使用该函数在处subs评估您的符号表达式。fka