Items tagged with chaos chaos Tagged Items Feed

I am going to use maple in connection with my project concerning non-linear dynamic and chaos, but i have never used the program before, so I hope you guys can help me with that,

I have a model of non-linear dynamic system as following:

dA/dr = k + m * Y(r) - A(r) * B(r) * B(r) - A(r)

dB/dr=1/q * [A(r) * B(r) * B(r) + A(r) - B(r)]

dY/dr=1/s * [B(r) - Y(r)] 

the values of the parameters are:
k= 10
q= 5 * 10^(-3)
s=2 * 10^(-2)

In this post I will describe a little about the OU course MS325: Computer Algebra, Chaos and Simulations, which I took last year.

MS325 is a level 3 OU applied mathematics course, which means, roughly that it is pitched at the level of a final year mathematics undergraduate. It is split into three components: Computer Algebra, which teaches the use of Maple and Maple programming; Chaos, which teaches dynamical systems, deterministic chaos and fractals, with an emphasis...

I spent this past week preparing a Webex presentation to a client who was interested in using Maple for a physics course in chaos. Of the two texts selected for the course, I had one on my own bookshelf. So I scanned Steven Strogatz' text Nonlinear Dynamics and Chaos (Addison Wesley, 1994) for topics that would profit from investigation with Maple.

Here is the Rossler system, one of the simplest examples of 3 dimensional deterministic chaos (under certain conditions according to "params"). Thanks to Doug and Joe for various assists. Comments and critiques most welcome !

restart;
interface(displayprecision=10):
ross_x:=diff(x(t),t)=-y(t)-z(t):
ross_y:=diff(y(t),t)=x(t)+a*y(t):
ross_z:=diff(z(t),t)=b+x(t)*z(t)-c*z(t):
rossler_sys:=ross_x,ross_y,ross_z;

#Find fixed points:
sol:=solve({rhs(ross_x...