Cálculo y python

Buenas!

Les cuento, para todos aquellos que no saben, que este lindo semestre estoy cursando cálculo numérico (y métodos numéricos / computacional). Entonces decidí que sería interesante aprender algo de python aplicado a cálculo. Razones:

  • No hay que andar declarando cosas.
  • Tiene un nombre chévere.
  • Su módulo math me convenció bastante.
  • En lignux muchas cosas funcionan bajo él.
  • Para hacerme el copado… (jaja! todos ustedes que hacen sus pseudoprogramas en excel).

En fin, aqui tenemos el primer código. Implementa bisección para una función determinada. Tiene parada por número de iteraciones y por |f(x)<e

¿Qué le falta?

Bien, algunos controles que te soliciten nuevos puntos si las ordenadas no tienen signo diferente. También, y lo más importante, algo para poder cambiar las funciones dinámicamente. Es decir, que las solicite al comienzo.

Just show the code…

#!/usr/bin/python
# -*- coding: iso8859-15 -*-
# Método de bisección
from math import *
#se importa math para dar soporte a funciones como sqrt y otras hierbas
def funx(x):
	return x**3-3*x-5 #función introducida
xa = float(raw_input("Introduzca margen izquierdo: "))
xb = float(raw_input("Introduzca margen derecho: "))
ya = funx(xa)
yb = funx(xb)
print "las ordenadas respectivas son:", ya,"y",yb
#criterios de parada
err = float(raw_input("introduzca e/ |f(r)|< e: "))
u = float(raw_input("Introduzca el número de iteraciones: "))
k = 0
yi = 100
while   (err < abs(yi)) and (k < u):
	ya = funx(xa)
	yb = funx(xb)
	xi = xb-(funx(xb)*(xb-xa)/(funx(xb)-funx(xa)))
	yi = funx(xi)
	if yi * ya < 0:
		xb = xi
	else:
		xa = xi
	k += 1

print "r es aproximado a %5.20f" % xi ,"y se obtuvo en la iteración número %d" % k
Consola corriendo el script

Consola corriendo el script

Puede descargarse desde acá Biseccion, listo para ejecutar.
¿Recomendaciones?

Espero seguir con mi progreso. Voy a ir presentando los scripts conforme transcurra la cátedra.

./.

Anuncios

About fclad

Estudiante de ingeniería, con todo lo que conlleva!
Esta entrada fue publicada en Comp, Labore, Proyectos y etiquetada , . Guarda el enlace permanente.

3 Responses to Cálculo y python

  1. Natu dice:

    jaja con el solo hecho de ser alumno de la facu de ingenieria ya sos copado(ni creidos eran ah?)bueno fer suerte con esto, espero ver como sigue.

  2. Sr009 dice:

    si bien el nombre es chevere tyambien es medio gay ya veremos su uso cuando lo aprenda a usar
    sigue asi

  3. Bubux dice:

    Es bien chevere el nombre. Me tenes que ensañar a usar python. Así te doy mas ayuda con el planet y esas fumadas. Nos vemos Fer, en algun asado supongo que espero organizar pronto

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s