Blogia
casanovas

Quantitative Art

El embudo de Gaston Julia

El embudo de Gaston Julia

#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)

Conjunto de Julia (Generado con XAOS)

Conjunto de Julia (Generado con XAOS)

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.

El Atractor de Lorenz con SCILAB

El Atractor de Lorenz con SCILAB

 

// 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