vladomiro
2f9bb6d4e4
modificados: .gitignore nuevos archivos: Mossbauer.jpg modificados: notas5.pdf modificados: notas5.tex modificados: pres7.pdf modificados: pres7.tex
22 lines
323 B
Python
Executable file
22 lines
323 B
Python
Executable file
import numpy as np
|
|
from random import random
|
|
|
|
|
|
#Definimos la funcion
|
|
def f(x):
|
|
return np.sin(1/(x*(2-x)))**2
|
|
|
|
#Definimos las inicializaciones
|
|
N= 1000
|
|
count = 0
|
|
|
|
#ciclo principal
|
|
for i in range(N):
|
|
x=2*random()
|
|
y= random()
|
|
if y<f(x):
|
|
count+=1
|
|
|
|
#Evaluamos la integral
|
|
I=2*count/N
|
|
print(I)
|