Ecuaciones Diferenciales con SAGE

 

Insc√≠bete a nuestro curso y lleva tus habilidades de Ecuaciones Diferenciales al siguiente nivel resolviendo un problema real utilizando el lenguaje de programaci√≥n del futuro: SAGE



\begin{eqnarray} \frac{dX}{d t} &=& \sigma(Y-X) \label{eq:X}\\ \frac{dY}{d t} &=& -XZ + rX -Y \label{eq:Y}\\ \frac{dZ}{d t} &=& XY - bZ \label{eq:Z} \end{eqnarray}

 def f(state, t):
x, y, z = state # unpack the state vector
  return sigma (y - x), x (rho - z) - y, x y - beta z # derivatives

  state0 = [1.0, 1.0, 1.0] t = np.arange(0.0, 40.0, 0.01)
  states = odeint(f, state0, t)

  fig = plt.figure()
  ax = fig.gca(projection='3d')
  ax.plot(states[:,0], states[:,1], states[:,2])
  plt.show()