Facebook Twitter Google +1     Admin

Se muestran los artículos pertenecientes a Agosto de 2009.

El Atractor de Lorenz con SCILAB

por Jorge Gonzalez Casanovas el 04/08/2009 21:05, en Quantitative Art

20090823211436-atractor-mariposa.jpg

 

// ODE definition

function ydot=lorenz(t,y)

x=y(1);

a=[-10,10,0;28,-1,-x;0,x,-8/3];

ydot=a*y

endfunction

function j=jacobian(t,y)

x=y(1);yy=y(2);z=y(3);

j=[-10,10,0;28-z,-1,-x;-yy,x,-8/3]

endfunction

/// Integration

y0=[-3;-6;12];t0=0;step=0.01;t1=10;

instants=t0:step:t1;

y=ode(y0,t0,instants,lorenz,jacobian);

//Visualization

clf()

a=gca();a.title.text=’Lorenz differential equation’;

a.title.font_size=4;

param3d(y(1,:),y(2,:),y(3,:))

//Después de copiar este código pulsar ENTER dos veces

 

Más...

Test de homogeneidad de frecuencias

por Jorge Gonzalez Casanovas el 11/08/2009 10:22, en sin tema

Si se quiere comparar las frecuencias de las 4 estaciones del año, con respecto a conteos de cualquier variable aleatoria. El código R sería el siguiente:

.Table <- matrix(c(24,32,10,9))
 colnames(.Table) <- c('1', '2', '3', '4')
 .Table  # Counts
 .Test <- chisq.test(.Table, correct=FALSE)
 .Test

Da como resultado:

    Chi-squared test for given probabilities

data:  .Table
X-squared = 19.9867, df = 3, p-value = 0.0001708

Más...



Tissue Modelling

por Jorge Gonzalez Casanovas el 16/08/2009 20:25, en sin tema

20090816202614-liver-tissue-modelling.png

Más...

por Jorge Gonzalez Casanovas el 23/08/2009 20:34, en sin tema

20090823205212-tumoral.jpg
Más...


Conjunto de Julia (Generado con XAOS)

por Jorge Gonzalez Casanovas el 23/08/2009 20:45, en Quantitative Art

20090823204523-fractal.png

Una familia muy notable de conjuntos de Julia se obtienen a partir de funciones cuadráticas simples: fc(z) = z2 + c, donde c es un número complejo. El conjunto de Julia que se obtiene a partir de esta función se denota Jc.

Más...

El embudo de Gaston Julia

por Jorge Gonzalez Casanovas el 24/08/2009 21:21, en Quantitative Art

20090824212150-julia.jpg

#ejemplo de Python para modelar el conjunto de Julia con
#variable externa generaciones. Lo que va dentro del
#FOR tiene que ir indentado, es decir tabulador o
#espacios
def crecimiento(generaciones):
    N=0
    C=complex(-1.3,0.006)
    
    for tiempo in range(generaciones):
        N=N**2+C
        print abs(N)  

crecimiento(100)

Más...