En estos dias he andado probando algunas cosas usando Python, la mayoría de las cosas con las que he estado jugando consisten básicamente en reescribir algunos programas que tenia escritos en C y también en realizar los programas que debería escribir en Matlab para mi clase de Biología Computacional.
Para mi sorpresa gracias al post anterior donde hable sobre la función de error complementario he recibido bastantes visitas, francamente no creí que hubiese mucha gente interesada sobre el tema y con esto ahora pondré la manera de como resolver numéricamente las integrales y a partir de ahí pueden construir su propio erfc.
- import scipy
- from scipy.integrate import quad #para resolver las integrales
- def intgrnd1(t):
- return exp(-t**2) #definir la funcion a integrar
- result, err = quad(intgrnd1,0.95,inf) #0.95 e inf son los limites inferior y superior
- print 'El resultado de la integral es ', result
Me ha venido genial tu post. Sigue así.
gracias :)