Question: Computation of left and right eigenvectors

Someone should please help me compute the left and right eigenvectors of the system below. The purpose is to compute values for 'a' and 'b' in the bifurcation formula.

Thank you

``

with(VectorCalculus)

 

(1)

interface(imaginaryunit = I)

I

(2)

I

I

(3)
 

diff(S(t), t) := `Λ__p`-(`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p]+µ__C)*S+`ω__B`*I__B

Lambda__p-(`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p]+µ__C)*S+omega__B*I__B

(4)

diff(I__B(t), t) := `#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S/N[p]-`ω__B`*I__B-(`σ__B`+µ__C)*I__B

`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S/N[p]-omega__B*I__B-(sigma__B+µ__C)*I__B

(5)

NULL

``

(6)

diff(S__A(t), t) := `Λ__A`-(µ__A+`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p])*S__A+`δ__A`*I__A

Lambda__A-(µ__A+`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p])*S__A+delta__A*I__A

(7)

diff(I__A(t), t) := `#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S__A/N[p]-(µ__A+`δ__A`)*I__A

`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S__A/N[p]-(µ__A+delta__A)*I__A

(8)

NULL

``

(9)

Download CBD2.mw

Please Wait...