Advanced Search

Question: How to assign parameters to a matrix?

Question:How to assign parameters to a matrix?

Peak King 5 Maple
assign
May 06 2011
1

I want to calculate  an eigenvector of a jacobian matrix J with fixed parameters from the solution of some equations, but the following code does not work. When I use the command in the last 3rd line, it works well.
But I want Maple to calculate J automatically, what shall I do?

restart;
with(linalg);

F[1] := A-(B+1)*X[1]+X[1]^2*X[2];
F[2] := B*X[1]-X[1]^2*X[2];

J := jacobian([F[1], F[2]], [X[1], X[2]]);

sol := solve({F[1] = 0, F[2] = 0, trace(J) = 0}, {B, X[1], X[2]});

assign(sol);
assume(A > 0);
# J := matrix(2, 2, [A^2, A^2, -1-A^2, -A^2])
ev := eigenvectors(J, 'radical');
q := ev[1][3][1];
Please Wait...