Algunas veces queremos obtener de manera cíclica resultados de algún problema o alguna función que estamos trabajando, en esta ocasión veremos como usar el bucle for en Python. Para ejemplificar esto usaremos la famosa y siempre bien ponderada función sin(x)
from numpy import *
for x in range(1, 10):
print repr(math.sin(x)).rjust(2), repr(math.sin(x*x)).rjust(8)
print repr(math.sin(x*x*x)).rjust(11)
math.sin es la manera que se escribe en Python las funciones matemáticas.
rjust() justifica a la derecha dando un grosor en especifico.
Así que vayan y abran su editor de código favorito y comiencen a probar lo escrito anteriormente.
Technorati Tags: for,, Python,, seno
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
Si lo que queremos comprobar si el valor de erfc obtenido a partir de la función incluida en Python podemos resolver las integrales de manera algebraica y después realizar la implementación utilizando parte del código de arriba.
So happy hacking!
Technorati Tags: Integrales,, Python,, erfc,, science
En mi clase de Biologia Computacional comenzaremos a jugar con erfc, aka función de error complementario 
Para la clase siguiente debo hablar sobre como se realiza la gráfica de esta función usando MATLAB, aprovechando que ando entusiasmado escribiendo código y más código de Python. Así que sin más preámbulos les presento la versión utilizando las utilidades que nos ofrece el paquete scipy.
from numpy import *
import scipy
import pylab
from scipy import special
loco = scipy.special.erfc(1.2)
print 'El resultado es', loco
La salida de este programa nos dice que:
El resultado es 0.0896860217704
En otro post pondré como resolver este mismo problema usando las definiciones formales, es decir usando integrales y raíces así como su respectiva representación gráfica. Código que les aseguro no abarcara más de 30 lineas.
Asimismo con este post inauguro una nueva categoría en este blog, obviamente llamada Python.
P.D. Si te pirateas el código, por favor deja un comentario.
Technorati Tags: erc,, matlab,, código,, code,, python,, scipy
