[pyar] Diversión para el finde: Acortando código

[Ariel Hernán Curiale]

El tema de numpy es la precisión (float64), pero con sympy hacemos algo similar a lo que hace tu decimal. En fin, algo con precisión super extendida sería

import numpy as np
from  sympy import *
print(N((((np.array(range(200)) +1)/ np.array([factorial(i) for i in range(200)])).sum() * N(.5,200)), 200))


Saludos
 ___________________________
| Ariel Hernán Curiale, Ph.D
| Research assistant at CONICET  
| Medical Physics Department 
| Centro Atómico Bariloche - FCEN,
| Univ. Nacional de Cuyo 
| Padre Jorge Contreras 1300
| 5500 Mendoza, Argentina
| Phone: +54 261 423 6003
| Web: www.curiale.com.ar <http://www.curiale.com.ar/>
|___________________________

> On Jun 17, 2017, at 08:42, Facundo Batista <facundobatista en gmail.com> wrote:
> 
> 2017-06-16 21:15 GMT-03:00 Ariel Hernán Curiale <curiale en gmail.com>:
> 
>> import numpy as np
>> from scipy.special import factorial
>> x=np.r_[0:150]
>> e=0.5* ((x+1)/factorial(x)).sum()
>> print(e)
> 
> ¡Muy bueno!
> 
> 
> -- 
> .    Facundo
> 
> Blog: http://www.taniquetil.com.ar/plog/
> PyAr: http://www.python.org/ar/
> Twitter: @facundobatista
> _______________________________________________
> Lista de Correo de PyAr - Python Argentina - pyar en python.org.ar
> Sitio web: http://www.python.org.ar/
> 
> Para administrar la lista (o desuscribirse) entrar a http://listas.python.org.ar/listinfo/pyar
> 
> La lista de PyAr esta Hosteada en USLA - Usuarios de Software Libre de Argentina - http://www.usla.org.ar

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://listas.python.org.ar/pipermail/pyar/attachments/20170617/87f46960/attachment.html>