x, y, z, t = var('x y z t') P = desolve_system_rk4([(9.0)*(y - x),x*(28.0 - z)- y, x*y - (8.0/3.0)*z], [x, y, z], ics = [0, 1.0, 1.0, 1.0], ivar = t, end_points = 40, step = 0.01) Q = [[j, k, l] for i, j, k, l in P] # colorindo o gráfico p = 0 for i in range(1,4): p += line( Q[i* 10^3 : (i + 1)*10^3 + 1 ], thickness = 2, color = hue(1./3*i)) p.show(spin = 2)
3D rendering not yet implemented
sage: x, y, z, t = var('x y z t') sage: P = desolve_system_rk4([(9.0)*(y - x),x*(28.0 - z)- y, x*y - (8.0/3.0)*z], [x, y, z], ics = [0, 1.0, 1.0, 1.0], ivar = t, end_points = 40, step = 0.01) sage: Q = [[j, k, l] for i, j, k, l in P] sage: # colorindo o gráfico sage: p = 0 sage: for i in range(1,4): ....: p += line( Q[i* 10^3 : (i + 1)*10^3 + 1 ], thickness = 3, color = hue(1./3*i)) sage: p.show() ︠f888bd6d-da89-4f14-984d-0a0337216563︠