In [6]:

    

import matplotlib.pyplot as plt
import numpy as np
from aleatorio import *

def AreaCirculo(radio):
    return np.pi*radio

radio=3
muestras = 1000
t=np.linspace(0,AreaCirculo(radio),100);
x=radio*np.cos(t);
y=radio*np.sin(t);

plt.plot(x,y,'-',2);

#Con el metodo MidSquare
#w = MidSquare("23",muestras)
#z = MidSquare("15",muestras)

#Con el metodo Congruencial
#w= Cm(205,muestras)
#z = Cm(45,muestras)

w= np.random.uniform(-radio, radio, size=1000)
z= np.random.uniform(-radio, radio, size=1000)
for i in range(muestras):
    if np.sqrt(pow(w[i],2)+pow(z[i],2)) < radio:
        plt.plot(w[i],z[i],'bo', alpha=.5)
    else:
        
        plt.scatter(w[i],z[i], marker="o",c='r')

plt.show()


    

In [ ]: