Я пишу функцию, которая вычисляет ряд Тейлора любой функции.
syms x
y=cos(x);
y0=0;
a=0;
for i=0:25
diff(y,i); %%Gives the derivative formula
y0=y0+diff(y,i)*((x-a)^i)/factorial(i); %%sums every new element of the series
end
x=0:0.1:2*pi;
res = subs(y0,x);
plot(x,res,x,cos(x))
Это код Matlab.
Моя проблема в том, что он отображает cos(2x)
вместо cos(x)
, аналогично он отображает ln(2x)
вместо ln(x)
и т. Д.
Я проверил факториалы, и они отображаются чтобы быть правильным. В чем может быть проблема, я испортил сериал или допустил ошибку Matlab?