function r=my_sum (n); x_values=0:.1:20; taylor_sin=zeros(length(x_values),1); for i=1:length(x_values) for j=1:n x(j)=[(-1)^(j-1)*(x_values(i))^(2*j-1)]/factorial(2*j-1); end; taylor_sin(i)=sum(x); end; length(taylor_sin) length(sin(x_values)) subplot(2,1,1); plot(x_values,taylor_sin) hold on; plot(x_values,sin(x_values)) subplot(2,1,2); plot(x_values,sin(x_values)-taylor_sin)