Question: Bifurcation Diagram from my system of equations

Please help with the bifurcation diagram for the system and parameter values below

NULL

with(VectorCalculus)

[`&x`, `*`, `+`, `-`, `.`, `<,>`, `<|>`, About, AddCoordinates, ArcLength, BasisFormat, Binormal, ConvertVector, CrossProduct, Curl, Curvature, D, Del, DirectionalDiff, Divergence, DotProduct, Flux, GetCoordinateParameters, GetCoordinates, GetNames, GetPVDescription, GetRootPoint, GetSpace, Gradient, Hessian, IsPositionVector, IsRootedVector, IsVectorField, Jacobian, Laplacian, LineInt, MapToBasis, Nabla, Norm, Normalize, PathInt, PlotPositionVector, PlotVector, PositionVector, PrincipalNormal, RadiusOfCurvature, RootedVector, ScalarPotential, SetCoordinateParameters, SetCoordinates, SpaceCurve, SurfaceInt, TNBFrame, TangentLine, TangentPlane, TangentVector, Torsion, Vector, VectorField, VectorPotential, VectorSpace, Wronskian, diff, eval, evalVF, int, limit, series]

(1)

interface(imaginaryunit = F)

I

(2)

M := Pi*theta-S*c__1-S*lambda+S__v*v__2

Pi*theta-S*c__1-S*lambda+S__v*v__2

(3)

Y := -S__v*c__2*lambda+Pi*b__1+S*v__1-S__v*c__3

-S__v*c__2*lambda+Pi*b__1+S*v__1-S__v*c__3

(4)

P := S__v*alpha+`&rho;__A`*A+c__4*`&rho;__Q`*Q+I*`&rho;__I`-µ*V

Q*c__4*rho__Q+A*rho__A+I*rho__I+S__v*alpha-V*µ

(5)

R := S__v*c__2*lambda-E*c__5+S*lambda

S__v*c__2*lambda-E*c__5+S*lambda

(6)

U := E*a*delta+Q*k*`&rho;__Q`-A*c__6

E*a*delta+Q*k*rho__Q-A*c__6

(7)

L := c__7*E-I*c__8

E*c__7-I*c__8

(8)

X := q__E*E+I*q__I-c__9*Q

E*q__E+I*q__I-Q*c__9

(9)

solve({L = 0, M = 0, P = 0, R = 0, U = 0, X = 0, Y = 0}, {I, A, E, Q, S, S__v, V})

{A = (a*c__8*c__9*delta+c__7*k*q__I*rho__Q+c__8*k*q__E*rho__Q)*lambda*Pi*(b__1*c__1*c__2+b__1*c__2*lambda+c__2*lambda*theta+c__2*theta*v__1+b__1*v__2+c__3*theta)/(c__6*c__9*c__5*c__8*(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2)), E = lambda*Pi*(b__1*c__1*c__2+b__1*c__2*lambda+c__2*lambda*theta+c__2*theta*v__1+b__1*v__2+c__3*theta)/(c__5*(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2)), I = c__7*lambda*Pi*(b__1*c__1*c__2+b__1*c__2*lambda+c__2*lambda*theta+c__2*theta*v__1+b__1*v__2+c__3*theta)/(c__5*c__8*(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2)), Q = (c__7*q__I+c__8*q__E)*lambda*Pi*(b__1*c__1*c__2+b__1*c__2*lambda+c__2*lambda*theta+c__2*theta*v__1+b__1*v__2+c__3*theta)/(c__9*c__5*c__8*(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2)), S = Pi*(c__2*lambda*theta+b__1*v__2+c__3*theta)/(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2), S__v = Pi*(b__1*c__1+b__1*lambda+theta*v__1)/(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2), V = Pi*(a*b__1*c__1*c__2*c__8*c__9*delta*lambda*rho__A+a*b__1*c__2*c__8*c__9*delta*lambda^2*rho__A+a*c__2*c__8*c__9*delta*lambda^2*rho__A*theta+a*c__2*c__8*c__9*delta*lambda*rho__A*theta*v__1+b__1*c__1*c__2*c__4*c__6*c__7*lambda*q__I*rho__Q+b__1*c__1*c__2*c__4*c__6*c__8*lambda*q__E*rho__Q+b__1*c__1*c__2*c__7*k*lambda*q__I*rho__A*rho__Q+b__1*c__1*c__2*c__8*k*lambda*q__E*rho__A*rho__Q+b__1*c__2*c__4*c__6*c__7*lambda^2*q__I*rho__Q+b__1*c__2*c__4*c__6*c__8*lambda^2*q__E*rho__Q+b__1*c__2*c__7*k*lambda^2*q__I*rho__A*rho__Q+b__1*c__2*c__8*k*lambda^2*q__E*rho__A*rho__Q+c__2*c__4*c__6*c__7*lambda^2*q__I*rho__Q*theta+c__2*c__4*c__6*c__7*lambda*q__I*rho__Q*theta*v__1+c__2*c__4*c__6*c__8*lambda^2*q__E*rho__Q*theta+c__2*c__4*c__6*c__8*lambda*q__E*rho__Q*theta*v__1+c__2*c__7*k*lambda^2*q__I*rho__A*rho__Q*theta+c__2*c__7*k*lambda*q__I*rho__A*rho__Q*theta*v__1+c__2*c__8*k*lambda^2*q__E*rho__A*rho__Q*theta+c__2*c__8*k*lambda*q__E*rho__A*rho__Q*theta*v__1+a*b__1*c__8*c__9*delta*lambda*rho__A*v__2+a*c__3*c__8*c__9*delta*lambda*rho__A*theta+b__1*c__1*c__2*c__6*c__7*c__9*lambda*rho__I+b__1*c__2*c__6*c__7*c__9*lambda^2*rho__I+b__1*c__4*c__6*c__7*lambda*q__I*rho__Q*v__2+b__1*c__4*c__6*c__8*lambda*q__E*rho__Q*v__2+b__1*c__7*k*lambda*q__I*rho__A*rho__Q*v__2+b__1*c__8*k*lambda*q__E*rho__A*rho__Q*v__2+c__2*c__6*c__7*c__9*lambda^2*rho__I*theta+c__2*c__6*c__7*c__9*lambda*rho__I*theta*v__1+c__3*c__4*c__6*c__7*lambda*q__I*rho__Q*theta+c__3*c__4*c__6*c__8*lambda*q__E*rho__Q*theta+c__3*c__7*k*lambda*q__I*rho__A*rho__Q*theta+c__3*c__8*k*lambda*q__E*rho__A*rho__Q*theta+alpha*b__1*c__1*c__5*c__6*c__8*c__9+alpha*b__1*c__5*c__6*c__8*c__9*lambda+alpha*c__5*c__6*c__8*c__9*theta*v__1+b__1*c__6*c__7*c__9*lambda*rho__I*v__2+c__3*c__6*c__7*c__9*lambda*rho__I*theta)/(c__5*c__6*c__8*c__9*µ*(c__1*c__2*lambda+c__2*lambda^2+c__1*c__3+c__3*lambda-v__1*v__2))}

(10)

``

lambda := beta*(I+`&eta;__A`*A+`&eta;__Q`*Q)/N

beta*(I+eta__A*A+eta__Q*Q)/N

(11)

``

NULL

k := .15

.15

(12)

delta := .125

.125

(13)

mu := 0.464360344e-4

0.464360344e-4

(14)

pi := .464360344

.464360344

(15)

delta__Q := 0.6847e-3

0.6847e-3

(16)

beta := .1086

.1086

(17)

q__E := 0.18113e-3

0.18113e-3

(18)

rho__Q := 0.815e-1

0.815e-1

(19)

a := .16255

.16255

(20)

v__1 := 0.5e-1

0.5e-1

(21)

v__2 := 0.5e-1

0.5e-1

(22)

alpha := 0.57e-1

0.57e-1

(23)

lambda := 0.765e-2

0.765e-2

(24)

rho__A := 0.915e-1

0.915e-1

(25)

rho__I := 0.515e-1

0.515e-1

(26)

a := .16255

.16255

(27)

q__I := 0.1923e-2

0.1923e-2

(28)

q__A := 0.4013e-7

0.4013e-7

(29)

eta__A := .1213

.1213

(30)

eta__Q := 0.3808e-2

0.3808e-2

(31)

w := .5925

.5925

(32)

Download Bifurcation_Equation.mw

Please Wait...