0
$\begingroup$
V = x^4 - m*x^3 + x^2
sol = Solve[D[V, x] == 0, x] // Simplify
Plot[{x/.sol[[1]],x/.sol[[2]],x/.sol[[3]]},{m,0.1,4}, PlotStyle->{Blue,Blue,Blue}]

as far as I'm concerned sol[[1]]=0 at first is stable, when m>mbif it becomes unstable

$\endgroup$
1
  • $\begingroup$ From your question it is unclear what you are asking for. $\endgroup$
    – yarchik
    Jan 23 at 12:32

1 Answer 1

2
$\begingroup$

I do not remember by heart the bifurcation value, MB, of the parameter m in this case, and too lazy to calculate it now. Let us assume that mb=2. Later you can substitute the correct value. With mb=2 try this:

ClearAll["Global`*"] 
V = x^4 - m*x^3 + x^2; 
sol = Solve[D[V, x] == 0, x] // Simplify;
mb = 2;
Show[{
  Plot[x /. sol[[1]], {m, 0.1, mb}, PlotStyle -> Blue, 
   PlotRange -> {{0.1, 4}, {-0.1, 2.5}}],
  Plot[x /. sol[[1]], {m, mb, 4}, PlotStyle -> Red],
  Plot[x /. sol[[3]], {m, 0.1, mb}, PlotStyle -> Red],
  Plot[x /. sol[[3]], {m, mb, 4}, PlotStyle -> Blue],
  Plot[x /. sol[[2]], {m, 0.1, 4}, PlotStyle -> Orange]
  }]

yielding the following

enter image description here

Have fun!

$\endgroup$

Not the answer you're looking for? Browse other questions tagged or ask your own question.