Store equilibrium points of a 3D system in a (x,y,z) array

[marellasunny]

I have a complicated 3d nonlinear differential system which I solve algebraically using maple's solve command. I aim to find the equilibrium points of the system using MAPLE and then plot these equilibrium points as a function of the 3 parameters alpha,beta,gamma(a.k.a.bifurcation diagram). The problem is that I have many equilibrium points as the solution and they are very lengthy in nature. In MATLAB,I would have accessed these solutions using STRUCTURES. How do I do so in MAPLE? Below,you will find part of the code with an form of the equilibrium points.

L1:=equation1=0
L2:=equation2=0
L3:=equation3=0

sys := {L1, L2, L3}; var := {x, y, z};
f3 := solve([L1, L2, L3], {x, y, z});
The output looks like this:
x1=fn(alpha,beta,gamma);y1=fn(alpha,beta,gamma);z1=fn(alpha,beta,gamma)...

x29=fn(alpha,beta,gamma);y29=fn(alpha,beta,gamma);z29=fn(alpha,beta,gamma);

Thats a lot of solutions and I would like to access them in (x,y,z) terms.How do I do that in MAPLE?

 

[Future Bruno]

You can access the solutions of your equations from the f3 output in MAPLE by using the subs command. For example: f3_x1 := subs(f3[1], x);This will give you the value of x1 as a function of alpha,beta and gamma. You can then use this value to substitute into other equations. For example, if you want to find the value of y1 as a function of alpha,beta and gamma, you can use the following command:f3_y1 := subs(f3[1], y);In order to plot the equilibrium points as a function of the 3 parameters alpha,beta and gamma, you can use the fsolve command. This command will allow you to solve the system of equations for a certain set of values given for alpha, beta and gamma. For example: fsolve({L1, L2, L3}, {x, y, z},{alpha=2,beta=3,gamma=4});This will solve the system of equations for the given values of alpha, beta and gamma and return the values of x, y and z as a list. You can then use these values to plot the equilibrium points on a bifurcation diagram.