"""confronto distribuzione binomiale p=q=1/2 - curva di Gauss """
from graftarta2 import Graf
import math, random, fatt
from graphics import *

g=Graf("Distr. binom. e curva di Gauss",600,300)
g.coordinate(-10,-0.1,40,0.4,True)
g.nascondi_tarta()

def rettangolo(altezza,base):
        for i in range(2):
            g.avanti(altezza)
            g.destra(90)
            g.avanti(base)
            g.destra(90)

def graf_isto(b,k):
	for n in b:
		rettangolo(n,k)
                g.destra(90)
		g.avanti(k)
		g.destra(-90)
        g.asPos(0,0)

def grafico_gauss(num):
	if num % 2 == 1:
		media = (num+1)/2
	else:
		media = num/2+0.5
	sigma = math.sqrt(num)/2
	x = 0
	while x<=40:
		y = math.exp(-(x-media)**2/(2*sigma**2))/(sigma*math.sqrt(2*math.pi))
		g.plot(x,y,nome)
		x = x+0.02
	
#numero casuale per la scelta di un colore casuale
rr=random.randint
i = 'w'
while i<>'f':
	x = input('num lanci ')
	b = fatt.distr_bin(x)
	nome= str(color_rgb(rr(0,255),rr(0,255),rr(0,255)))
	#color_rgb modulo graphics ritorna il colore
	print nome
        g.colorePenna(nome)
        graf_isto(b,1)
        grafico_gauss(x)
        i = raw_input('f e invio per finire ')
g.close()