> |
with(DEtools): with(plots): |

> |
Lorenz := diff(x(t), t) = sigma*(y(t)-x(t)), diff(y(t), t) = r*x(t)-y(t)-x(t)*z(t), diff(z(t), t) = x(t)*y(t)-b*z(t); |

Lorenz := diff(x(t), t) = sigma*(y(t)-x(t)), diff(y(t), t) = r*x(t)-y(t)-x(t)*z(t), diff(z(t), t) = x(t)*y(t)-b*z(t) |

> |
params := sigma = 10, r = 21, b = 8/3; |

params := sigma = 10, r = 21, b = 8/3 |

> |
init1 := x(0) = 4, y(0) = 5, z(0) = 6; |

init1 := x(0) = 4, y(0) = 5, z(0) = 6 |

> |
Lorenz1 := [op(subs(params, [Lorenz])), init1]; |

Lorenz1 := [diff(x(t), t) = 10*y(t)-10*x(t), diff(y(t), t) = 21*x(t)-y(t)-x(t)*z(t), diff(z(t), t) = x(t)*y(t)-(8/3)*z(t), x(0) = 4, y(0) = 5, z(0) = 6] |

> |
LS1 := dsolve(Lorenz1, numeric): |

> |
odeplot(LS1, [x(t), z(t)], 0 .. 30, numpoints = 2400); |

> |

> |