¿Se puede resolver el modelo de Lorenz usando Mathematica?

Sí.

NDSolve[{x'[t] == ​​​​-10(x[t] - y[t]), y'[t] == ​​-x[t] z[t ] + 28 x [t] - y[ t], z'[t] == ​​​​x[t] y[ t] - (8/3)z[t], x[0] == z[0 ] == 0, y [0] == 1}, {x, y, z}, {t, 0, 200}, MaxSteps -> Infinito]

ParametricPlot3D[Evaluar[{x[t ], y[t] , z[t]} /. %], {t, 0, 200}, PlotPoints -> 50000]

Puedes ajustar los parámetros.

Consulta mi blog para más detalles.

/%CD%F2%B6%BE%BF%F1%C4%A7/blog/item/75344d138e806c2bdc540107.html