Why repost?

http://www.mapleprimes.com/questions/145196-How-To-Compare-Or-Perform-Phase-Portrait--Diffthetat

For help with the syntax you could look at my answer to your question about the damped driven pendulum.

http://www.mapleprimes.com/questions/145198-The-Damped-Driven-Pendulum

I don't get anything I can see.

Could you give us the equations in Maple syntax in text form we can copy and paste or as an uploaded worksheet?

I don't see your equations. Could you give us them in text form or as an uploaded worksheet?

"but throughout mathematics teachings I've always learned cos(45) as 1/sqrt(2) as I'm sure

the rest of you all have as well"

Not this one!

I don't know anything about CUDA, but from what I can read in the help pages this technology accelerates floating point matrix computations only.

One of the eigenvalues is zero, since you make sure that det(K.M) = det(K)*det(M) = 0. The corresponding eigenvectors are computed fast:

restart;
M := Matrix(3, symbol = m, shape = symmetric);
det := LinearAlgebra:-Determinant(M);
M[1, 2] := solve(det, m[1, 2])[1];
K := Matrix(3, symbol = k, shape = symmetric);

LinearAlgebra:-LinearSolve(K.M, <0,0,0> ,free=t);

@dhonkabulo I have posted my answer to this as an 'answer' and not a 'comment'.  The system happens to be so simple that no numerical computation is necessary. That makes the solution simpler. I have done that in my newly posted answer.

@dhonkabulo I have posted my answer to this as an 'answer' and not a 'comment'.  The system happens to be so simple that no numerical computation is necessary. That makes the solution simpler. I have done that in my newly posted answer.

However, with that approach also

w();

and

w(77,bb,pp); #or whatever

would result in 99.

However, with that approach also

w();

and

w(77,bb,pp); #or whatever

would result in 99.

Could you give us text in MaplePrimes or upload a worksheet instead of a picture?

@Markiyan Hirnyk Nonlinearity has nothing to do with. The jacobian is the result of linearizing about the equilibrium point. There would be no point in linearizing if the system was linear.

I resent the tone of your comment, as I have done with some of your comments in the past to me and to other people.

@Markiyan Hirnyk Nonlinearity has nothing to do with. The jacobian is the result of linearizing about the equilibrium point. There would be no point in linearizing if the system was linear.

I resent the tone of your comment, as I have done with some of your comments in the past to me and to other people.

@Markiyan Hirnyk Stability questions are sometimes difficult to settle, yes. But very often it is a very simple matter as is the case with the two points found by fsolve (still assuming the order given in my answer):

subs(sol1,[v1,v2,u1]);
J(op(%));
LinearAlgebra:-Eigenvalues(%);
#So that point is asymptotically stable.
subs(sol2,[v1,v2,u1]);
J(op(%));
LinearAlgebra:-Eigenvalues(%);
#So that point is unstable.

@Markiyan Hirnyk Stability questions are sometimes difficult to settle, yes. But very often it is a very simple matter as is the case with the two points found by fsolve (still assuming the order given in my answer):

subs(sol1,[v1,v2,u1]);
J(op(%));
LinearAlgebra:-Eigenvalues(%);
#So that point is asymptotically stable.
subs(sol2,[v1,v2,u1]);
J(op(%));
LinearAlgebra:-Eigenvalues(%);
#So that point is unstable.