Hello,
I am trying to program the function exp(-x) using the power series and for loop and plotting it against the actual exp(-x) to check what the error is, however, for some reason the program does not show the correct graph for the power series expansion. My code is below:
x1 = 0:0.1:5;y = exp(-x1);figureplot(x,y,'k') % plots the true function e^-x in black
hold on;ex1 = ones(1,51);ex2= 0; for n = 1:5 ex1 = ex1 + ((-ex1).^n).*((x1.^n)/(factorial(n))) end plot(x1,ex1,'r');
The issue is with the expression for ex1, but I can't figure out why it isn't working as I am sure that I have used the correct formulae. Thanks.
Best Answer