Этот код не дает мне соз (х) в Matlab

Question

Я изменил код греха (х) из книги Отто, чтобы получить соз (х), и я сделал следующее:

v=0.:0.25:0.75; 
cosx=zeros(size(v)); 
N=10; 
range=0:N; 
ints=2*range; 

for n=range 
    cosx=cosx + (-1)^n*v.^ints(n)/factorial(ints(n)); 
end 

Но что это дает 0 всюду.

Может кто-нибудь помочь мне исправить это, пожалуйста? Большое спасибо

Answer 1

Вы должны изменить две вещи в коде, чтобы получить правильное расширение серии Taylor. Вот код. Обратите внимание на изменения.

v=0:0.25:0.75; 
cosx=zeros(1,size(v,2)); %Now, it's a vector 
N=10; 
range=0:N; 
ints=2*range; 

for n=range 
    cosx=cosx + (-1)^n*v.^ints(n+1)/factorial(ints(n+1)); % replace n by n+1,   
                 % MATLAB indexing starts from 1. 
end