Advanced Search

Question: Making a plot with the result of the greatest eigenvalue and variation of parameters of MAtrix AB inverse

Question:Making a plot with the result of the greatest eigenvalue and variation of parameters of MAtrix AB inverse

ogunmiloro 120 Maple 18
parameters linear-algebra plotting
April 03 2020
1


 

> 

with(LinearAlgebra); S[p] := 34.722406639004; alpha[1] := 0.2e-3; mu := 0.2041e-1; tau := .33; beta := .5; eta[1] := 0.96e-1; alpha[2] := .2; sigma := .9; e[o] := .33; delta := .2115; eta[2] := 0.2485e-2; A := Matrix(6, 6, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = alpha[1]*S[p], (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = beta, (4, 4) = 0, (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0}); B := Matrix(6, 6, {(1, 1) = mu, (1, 2) = tau, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = e[o], (2, 1) = 0, (2, 2) = tau+mu, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = beta+eta[1]+sigma+mu, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = (1-delta)*alpha[2]-mu, (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = alpha[2]*delta, (5, 5) = mu+eta[2], (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = mu+e[o]}); 1/B; VectorMatrixMultiply(A, 1/B)

S[p] := 34.722406639004

  

alpha[1] := 0.2e-3

  

mu := 0.2041e-1

  

tau := .33

  

beta := .5

  

eta[1] := 0.96e-1

  

alpha[2] := .2

  

sigma := .9

  

e[o] := .33

  

delta := .2115

  

eta[2] := 0.2485e-2

  

A := Matrix(6, 6, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0.6944481328e-2, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = .5, (4, 4) = 0, (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 0, (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = 0})

  

B := Matrix(6, 6, {(1, 1) = 0.2041e-1, (1, 2) = .33, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = .33, (2, 1) = 0, (2, 2) = .35041, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1.51641, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = .13729, (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0.4230e-1, (5, 5) = 0.22895e-1, (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = .35041})

  

Matrix(6, 6, {(1, 1) = 48.9955903968643, (1, 2) = -46.1417905623847, (1, 3) = 0., (1, 4) = 0., (1, 5) = 0., (1, 6) = -46.1417905623847, (2, 1) = 0., (2, 2) = 2.85379983447961, (2, 3) = 0., (2, 4) = 0., (2, 5) = 0., (2, 6) = -0., (3, 1) = 0., (3, 2) = 0., (3, 3) = .659452258953713, (3, 4) = 0., (3, 5) = 0., (3, 6) = -0., (4, 1) = 0., (4, 2) = 0., (4, 3) = 0., (4, 4) = 7.28385170077937, (4, 5) = 0., (4, 6) = -0., (5, 1) = 0., (5, 2) = 0., (5, 3) = 0., (5, 4) = -13.4573892528049, (5, 5) = 43.6776588774842, (5, 6) = -0., (6, 1) = 0., (6, 2) = 0., (6, 3) = 0., (6, 4) = 0., (6, 5) = 0., (6, 6) = 2.85379983447961})

  

Matrix(6, 6, {(1, 1) = 0., (1, 2) = 0., (1, 3) = 0., (1, 4) = 0., (1, 5) = 0., (1, 6) = -0., (2, 1) = 0., (2, 2) = 0., (2, 3) = 0., (2, 4) = 0., (2, 5) = 0., (2, 6) = -0., (3, 1) = 0., (3, 2) = 0., (3, 3) = 0.457955389901148e-2, (3, 4) = 0., (3, 5) = 0., (3, 6) = -0., (4, 1) = 0., (4, 2) = 0., (4, 3) = .329726129476857, (4, 4) = 0., (4, 5) = 0., (4, 6) = -0., (5, 1) = 0., (5, 2) = 0., (5, 3) = 0., (5, 4) = 0., (5, 5) = 0., (5, 6) = -0., (6, 1) = 0., (6, 2) = 0., (6, 3) = 0., (6, 4) = 0., (6, 5) = 0., (6, 6) = -0.})

 (1)
> 

LinearAlgebra:-Eigenvalues( (1) );

Vector(6, {(1) = 0.+0.*I, (2) = 0.+0.*I, (3) = 0.+0.*I, (4) = 0.457955389901148e-2+0.*I, (5) = 0.+0.*I, (6) = -0.+0.*I})

 (2)
> 


 

Download AB.mw
Please Wait...