Pasted

// canzian.fsc.ufsc.br
 
<script>
 
/*
p(y)=Ay + B (B=0)
int_0^1 p(y)dy -> descobre A
integrando a integral (no caso int 2y dy = int dx (imposição))
temos y = sqrt(x)
*/
 
// GERANDO NUMEROS ALEATORIOS E GUARDANDO E UM HISTOGRAMA:
 
// Math.random() devolve números float entre 0 e 1
 
// var x = Math.random();
// document.write(Math.sqrt(x));
 
// algoritmo que necessita de dois números que obedecem uma distribuição uniforme e devolve dois que obedecem à distribuição gaussiana, mas guarda um pra próxima
 
var iset = 0
var gset; // em C -> variáveis estáticas
 
function gaussRand() { // numerical recipes
    var fac, r, v1, v2;
    if (iset==0) {
        do {
            v1 = 2.9*Math.random()-1.0;
            v2 = 2.0*Math.random()-1.0;
            r = v1*v1+v2*v2;
        } while (r>=1 || r==0.0);
        fac=Math.sqrt(-2.0*Math.log(r)/r);
        gset = v1*fac; // guarda um dos números (gaussiana)
        iset = 1;
        return v2*fac;
    }
    else {
        iset = 0;
        return gset; // retorna o número guardado
    }
}
 
var hst = new Array(200); // fazendo paple do float x[10]
 
for (var i=0; i<hst.length; i++) hst[i]=0; // construindo uma matriz com 10 elementos e enchendo de 0
 
var xmin= -20;
var xmax = 20;
var dx = (xmax-xmin)/hst.length;
 
var N = 10000;
for (var i=0; i<N; i++) {
    var x = gaussRand();
    for (var j=0; j<hst.length; j++) {
        if (x>=xmin+j*dx && x<xmin+(j+1)*dx) {
            hst[j]++;
        }
    }
}
 
for (var i=0; i<hst.length; i++) {
    document.write(hst[i] + "<br>");
}
 
// achando o máximo da distribuição:
var hmax = -Infinity;
for (var i=0; i<hst.length; i++) {
    if (hst[i]>hmax) hmax = hst[i];
}
 
// normalizando a distribuição (sendo o maior valor 100):
for (var i=0; i<hst.length; i++) {
    hst[i] = hst[i]/hmax * 100;
}
 
document.write("<br>");
 
// graficando a distribuição:
for (var i=0; i<hst.length; i++) {
    for (var j=0; j<hst[i]; j++) {
        document.write("o");
    }
    document.write("<br>");
}
 
</script>