Dear mapleprimes ODE enthusiasts,

I am analyzing a stiff 3D system characterizes by fast-slow dynamics. I believe I can show analytically that, for some parameter values, in a neighborhood of some critical points (where the system is stationary), the dynamic system has a local center manifold that is center-stable. I have simulated the system for different parameter values.I am interested in one particular trajectory. With a random set of initial values, the system...

dsolve event handling looks like a great tool, but I can't understand how it works and it's not sufficiently popular to google my way to examples online. There are two things I'd like to do: 1) interrupt the computation once some variable leaves a certain range, 2) interrupt the computation if the solver gets stuck into an infinite loop. Both are documented in the help, but I don't understand it. Suggestions welcome, thanks!

http://www.maplesoft.com/support/help/Maple/view...

I am computing eigenvalues of a 3X3 square Matrix that contains symbolic elements and many zeros. There are 2 zero eigenvalues and 1 non-zero. I'd like to request Maple to order the eigenvalues in a certain way, with the non-zero eigenvalue first. Oddly, the command LinearAlgebra:-Eigenvectors(A); orders them in a seemingly random way, with the non-zero eigenvalue either first or last, never in the middle. The LinearAlgebra:-Eigenvalues(A); commands seems...

There was a recent exchange I was a part of on Maple classic versus standard, etc. It disappeared from the list of "recent discussions" so I tried to track it by other means. A keyword search did not immediately return a helpful link, so I thought I would click on my profile and look for my recent comments from the lists there. But I could only find a list of "questions," "answers" and "posts". My contribution there was merely a "comment" and I can't seem to find any list about those...

I have this non-linear ODE in q(p):

ode := diff(q(p), p) = (q(p)-a*p)/(p-a*q(p))*(p^(-b)-1)/q(p)^(-b);

for some real parameters a and b, it may be solved in exact form for the special value a=1.

Is there any hope of getting a closed-form solution or of expressing the solution in a neat way for the case 0<a<1? even some specific value of a would be of interest, e.g. a=1/2.

dsolve(eval(ode,{a=1}));