How can I solve each element of an array symbolically?

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]

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