0
私はShijun Liaoの "非線形微分方程式のホモトピー解析法"というタイトルの本を読んでいます。第13章金融のアプリケーション：アメリカン・プット・オプションには、この章のmathematicaコードに関する付録があります。しかし、私は代わりにMathematica butMatlabを持っていないので、Matlabのコードを書き直したいと思います。ここMatlabの任意の関数のテイラー展開を構築するには？
は、Mathematicaコードの最初の部分である：
<<Calculus`Pade`;
<<Graphics`Graphics`;
(* Define approx[f] for Taylor expansion of f *)
approx[f_] := Module[{temp},
temp[0] = Series[f, {t, 0, OrderTaylor}]//Normal;
temp[1] = temp[0] /. t^(n_.)*Derivative[j_][DiracDelta][0] -> 0;
temp[2] = temp[1] /. t^(n_.)*DiracDelta[0] -> 0;
temp[3] = temp[2] /. DiracDelta[0] -> 0;
temp[4] = temp[3] /. Derivative[j_][DiracDelta][0] -> 0;
temp[5] = N[temp[4],60]//Expand;
If[KeyCutOff == 1, temp[5] = temp[5]//Chop];
temp[5]
];
Iシリーズ[F {X、X 0、N}]は点約fに対するベキ級数展開を生成することのWolfram参照サイトから知りますx = x0を（x-x0）^ nとする。したがって、私はfの拡張のべき級数を構築する必要があります。しかし、限り私が知っているように、私たちは、MATLABでfを最初に、例えばf = cos（x）などを定義する必要がありますので、私の質問は、Matlabの任意の関数のテイラー展開を構築する方法です？
ありがとうございます。