Trying to write a script that solves a series
显示 更早的评论
The series goes 1-1/3!+1/5!-1/7!+.....20 terms
This is my first guess on how to solve it. Don't have software available to me at the moment. This is a start but I have a feeling something isn't right here.
a0=1/factorial(0); %first term
an=1/factorial(38); %final term
n=0; Initialize
sumf=0; %initialize sum of series
sign=1; %initialize sign
for tn=a0: -1 / factorial(n+2):an ; %step size negative due to decreasing values.
sumf = (sumf + tn) . * sign;
sign= (sign .* -1); %update sign
n=n+2; %The n is update within the factorial increment
end
any problems with this? can you update an increment?
回答（1 个）
A few problems, mostly regarding numerics
- A factorial can only be represented accurately in double precision up to 21. Same applies for the inverse
- Your different terms will span many orders of magnitude. Adding terms like that kills accuracy.
If you want accurate results and have the symbolic math toolbox you could look into vpa().
类别
在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
