4828.sagews

#Atractor de Lorenz


h=0.01;        
k=2000        

sigma=10;       
rho=28;
beta=8/3;

x=1; y=1; z=1;         

puntos=[]
for i in range(k):
    x,y,z=x+h*( sigma*(y-x) ), y+h*( x*(rho - z) - y ), z+h*( x*y - beta*z )
    puntos.append((x,y,z))
point3d(puntos)