How can I solve each element of an array symbolically?

Question

Demetri Pananos 2015-10-9

回答： Walter Roberson 2015-10-9

I'm trying to find bifurcation points of a dynamical system.

I calculate the Jacobian of my system (which I call JJ), substitute in fixed points and find the eigenvalues (which I call L).

I would like to then solve each element of L, and find the smallest solution. This will give me the bifurcation value.

The solution may return 0, 1, or many solutions. How can I do this efficiently?

Answer 1

Walter Roberson 2015-10-9

For loop over length(L), index L, solve, sort through the solutions

minsol = inf;
for K = 1 : length(L)
  s = solve(L(K));
  mins = min(s);
  if mins < minsol; minsol = mins; end
end