What is the error in this code.

%expansion of maclaurine series
function result =  expansion(n,x)

c = [1];
d = [1:n];
a = [x.^d];
for i=1:n
    vec = [vec,1./prod(1:i)];
end
e = sum(a.*vec);
result = [1+sum];
end

Answer 1

[James Tursa]

Did you mean:

result = 1 + e;