Software

Buenas noches continuando con la marathon de publicaciones en mi blog jeje es broma muestro ahora uno de mis tantos codigos que he hecho en momentos de inspiracion y como muchos me dicen en momentos de vagancia =)

1. #Programa: Integrales por metodo de Simpson
2. #Programador: William Mendez
3. #Fecha: 14/01/2012
4. #Objetivo: """Crear un programa que permita resolver las integrales por medio de la regla de simpson"""
5. #Lenguaje: Python 2.7.2
7. def funcion(x):
8.    integral = 1/(1+x)#"""aqui ingrese la funcion que desea integrar"""
9.    return integral
11. a = float (raw_input("Ingrese limite inferior: "))
12. b = float (raw_input("Ingrese limite Superior: "))
13. n = int( raw_input("Numero par de intervalos: "))
14. h = (b-a)/(2*n)
15. n1 = 2*n-1
16. n2 = 2*n-1
17. j = 1
18. si = 0
19. while j<=n1:
20.   x = a + h * j
21.   si += funcion(x)
22.   j+=2
24. j = 2
25. sp = 0
26. while j<= n2:
27.    x = a + h * j
28.    sp += funcion(x)
29.    j+=2
30. is = h*(funcion(a)+4*si + 2*sp+funcion(b))/3
31. print"La integral es: "+is+"\n"
32. print"pulse cualquier tecla para salir...."
33. raw_input()