.Alternate to using for loop or symsum for the summation ∑(const)^n/(n*n!) ?
7 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Dear all,
Is there a more computationally efficient way compared to using for loop or symsum (from Symbolic math toolbox) to compute:
∑(const)^n/(n*n!)
const is some constant value, n is the range of limit varying from 1 to infinity (or some high value like 200 for approximating the sum).
-- Thanks, Ram.
2 Commenti
Sean de Wolski
il 26 Giu 2013
Modificato: Sean de Wolski
il 26 Giu 2013
Why not symsum? You're going to need it for factorial greater than 170 anyway:
factorial(171)
Risposta accettata
Roger Stafford
il 26 Giu 2013
Your sum is equal to the integral
int('(exp(x)-1)/x','x',0,const)
so you could do numerical integration of this rather than summing the infinite series. That integrand is actually well-behaved in the vicinity of x = 0, but computing it might give you some problems, so you could substitute a Taylor series approximation very near x = 0.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Calculus in Help Center e File Exchange
Prodotti
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!