Pemudahijrah01

25 Reputation

4 Badges

3 years, 268 days

MaplePrimes Activity


These are questions asked by Pemudahijrah01

this is my model. Please give me how to find the DFE and basic reproduction number from maple.

please help me to fixed and find the analytic solution...

restart

with(linalg)

f1 := mu[1]*N-delta*r-beta*N*s(e+i)

mu[1]*N-delta*r-beta*N*s(e+i)

(1)

f2 := beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

(2)

``

f3 := pi*e-(mu[1]+mu[2]+sigma)*i

pi*e-(mu[1]+mu[2]+sigma)*i

(3)

f4 := gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

(4)

f5 := theta*q+xi*v-(mu[1]+delta)*r

theta*q+xi*v-(mu[1]+delta)*r

(5)

f6 := eta*N*s-(mu[1]+omega+xi)*v

eta*N*s-(mu[1]+omega+xi)*v

(6)

T := solve({f1, f2, f3, f4, f5, f6}, [s, e, i, q, r, v])

Error, (in solve) cannot solve for an unknown function with other operations in its arguments

 

``

Download cobacoba.mw

How I can solve this, because I can't find the solutions. The display show "Length of output exceeds limit of 1000000"

Please help me.. 

This is my model on picture

please, help me..this is my model

How I can change this solution or remove the RootOf, because i have to show variables and params in this solutions.
 

restart

with(linalg)

f1 := alpha*x*(1-x/N)-beta[1]*sqrt(x)*y1/(1+h[1]*beta[1]*sqrt(x))-beta[2]*sqrt(x)*y2/(1+h[2]*beta[2]*sqrt(x))-d*E*x

f2 := -omega[1]*y1+Mu[1]*beta[1]*sqrt(x)*y1*(1-y1/(beta[1]*sqrt(x)))/(1+h[1]*beta[1]*sqrt(x))

f3 := -omega[2]*y2+Mu[2]*beta[2]*sqrt(x)*y2*(1-y2/(beta[2]*sqrt(x)))/(1+h[2]*beta[2]*sqrt(x))

T := solve({f1, f2, f3}, [x, y1, y2])

[[x = -N*(E*d-alpha)/alpha, y1 = 0, y2 = 0], [x = RootOf(alpha*beta[2]*Mu[2]*h[2]*_Z^4+alpha*Mu[2]*_Z^3+(E*N*d*Mu[2]*beta[2]*h[2]-N*alpha*Mu[2]*beta[2]*h[2])*_Z^2+(-N*beta[2]^2*h[2]*omega[2]+E*N*d*Mu[2]+N*Mu[2]*beta[2]^2-N*alpha*Mu[2])*_Z-N*omega[2]*beta[2])^2, y1 = 0, y2 = -(RootOf(alpha*beta[2]*Mu[2]*h[2]*_Z^4+alpha*Mu[2]*_Z^3+(E*N*d*Mu[2]*beta[2]*h[2]-N*alpha*Mu[2]*beta[2]*h[2])*_Z^2+(-N*beta[2]^2*h[2]*omega[2]+E*N*d*Mu[2]+N*Mu[2]*beta[2]^2-N*alpha*Mu[2])*_Z-N*omega[2]*beta[2])*omega[2]*beta[2]*h[2]-RootOf(alpha*beta[2]*Mu[2]*h[2]*_Z^4+alpha*Mu[2]*_Z^3+(E*N*d*Mu[2]*beta[2]*h[2]-N*alpha*Mu[2]*beta[2]*h[2])*_Z^2+(-N*beta[2]^2*h[2]*omega[2]+E*N*d*Mu[2]+N*Mu[2]*beta[2]^2-N*alpha*Mu[2])*_Z-N*omega[2]*beta[2])*beta[2]*Mu[2]+omega[2])/Mu[2]], [x = RootOf(alpha*beta[1]*Mu[1]*h[1]*_Z^4+alpha*Mu[1]*_Z^3+(E*N*d*Mu[1]*beta[1]*h[1]-N*alpha*Mu[1]*beta[1]*h[1])*_Z^2+(-N*beta[1]^2*h[1]*omega[1]+E*N*d*Mu[1]+N*Mu[1]*beta[1]^2-N*alpha*Mu[1])*_Z-N*omega[1]*beta[1])^2, y1 = -(RootOf(alpha*beta[1]*Mu[1]*h[1]*_Z^4+alpha*Mu[1]*_Z^3+(E*N*d*Mu[1]*beta[1]*h[1]-N*alpha*Mu[1]*beta[1]*h[1])*_Z^2+(-N*beta[1]^2*h[1]*omega[1]+E*N*d*Mu[1]+N*Mu[1]*beta[1]^2-N*alpha*Mu[1])*_Z-N*omega[1]*beta[1])*omega[1]*beta[1]*h[1]-RootOf(alpha*beta[1]*Mu[1]*h[1]*_Z^4+alpha*Mu[1]*_Z^3+(E*N*d*Mu[1]*beta[1]*h[1]-N*alpha*Mu[1]*beta[1]*h[1])*_Z^2+(-N*beta[1]^2*h[1]*omega[1]+E*N*d*Mu[1]+N*Mu[1]*beta[1]^2-N*alpha*Mu[1])*_Z-N*omega[1]*beta[1])*beta[1]*Mu[1]+omega[1])/Mu[1], y2 = 0], [x = RootOf(alpha*beta[2]*beta[1]*Mu[2]*Mu[1]*h[2]*h[1]*_Z^5+(alpha*Mu[1]*Mu[2]*beta[1]*h[1]+alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^4+(E*N*d*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]-N*alpha*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]+alpha*Mu[1]*Mu[2])*_Z^3+(-N*Mu[1]*beta[1]*beta[2]^2*h[1]*h[2]*omega[2]-N*Mu[2]*beta[1]^2*beta[2]*h[1]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]*beta[1]*h[1]+E*N*d*Mu[1]*Mu[2]*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]^2*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]*beta[2]^2*h[1]-N*alpha*Mu[1]*Mu[2]*beta[1]*h[1]-N*alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^2+(-N*Mu[1]*beta[1]*beta[2]*h[1]*omega[2]-N*Mu[1]*beta[2]^2*h[2]*omega[2]-N*Mu[2]*beta[1]^2*h[1]*omega[1]-N*Mu[2]*beta[1]*beta[2]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]+N*Mu[1]*Mu[2]*beta[1]^2+N*Mu[1]*Mu[2]*beta[2]^2-N*alpha*Mu[1]*Mu[2])*_Z-N*omega[2]*beta[2]*Mu[1]-N*omega[1]*beta[1]*Mu[2])^2, y1 = (-RootOf(alpha*beta[2]*beta[1]*Mu[2]*Mu[1]*h[2]*h[1]*_Z^5+(alpha*Mu[1]*Mu[2]*beta[1]*h[1]+alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^4+(E*N*d*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]-N*alpha*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]+alpha*Mu[1]*Mu[2])*_Z^3+(-N*Mu[1]*beta[1]*beta[2]^2*h[1]*h[2]*omega[2]-N*Mu[2]*beta[1]^2*beta[2]*h[1]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]*beta[1]*h[1]+E*N*d*Mu[1]*Mu[2]*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]^2*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]*beta[2]^2*h[1]-N*alpha*Mu[1]*Mu[2]*beta[1]*h[1]-N*alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^2+(-N*Mu[1]*beta[1]*beta[2]*h[1]*omega[2]-N*Mu[1]*beta[2]^2*h[2]*omega[2]-N*Mu[2]*beta[1]^2*h[1]*omega[1]-N*Mu[2]*beta[1]*beta[2]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]+N*Mu[1]*Mu[2]*beta[1]^2+N*Mu[1]*Mu[2]*beta[2]^2-N*alpha*Mu[1]*Mu[2])*_Z-N*omega[2]*beta[2]*Mu[1]-N*omega[1]*beta[1]*Mu[2])*omega[1]*beta[1]*h[1]+RootOf(alpha*beta[2]*beta[1]*Mu[2]*Mu[1]*h[2]*h[1]*_Z^5+(alpha*Mu[1]*Mu[2]*beta[1]*h[1]+alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^4+(E*N*d*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]-N*alpha*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]+alpha*Mu[1]*Mu[2])*_Z^3+(-N*Mu[1]*beta[1]*beta[2]^2*h[1]*h[2]*omega[2]-N*Mu[2]*beta[1]^2*beta[2]*h[1]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]*beta[1]*h[1]+E*N*d*Mu[1]*Mu[2]*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]^2*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]*beta[2]^2*h[1]-N*alpha*Mu[1]*Mu[2]*beta[1]*h[1]-N*alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^2+(-N*Mu[1]*beta[1]*beta[2]*h[1]*omega[2]-N*Mu[1]*beta[2]^2*h[2]*omega[2]-N*Mu[2]*beta[1]^2*h[1]*omega[1]-N*Mu[2]*beta[1]*beta[2]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]+N*Mu[1]*Mu[2]*beta[1]^2+N*Mu[1]*Mu[2]*beta[2]^2-N*alpha*Mu[1]*Mu[2])*_Z-N*omega[2]*beta[2]*Mu[1]-N*omega[1]*beta[1]*Mu[2])*beta[1]*Mu[1]-omega[1])/Mu[1], y2 = (-RootOf(alpha*beta[2]*beta[1]*Mu[2]*Mu[1]*h[2]*h[1]*_Z^5+(alpha*Mu[1]*Mu[2]*beta[1]*h[1]+alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^4+(E*N*d*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]-N*alpha*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]+alpha*Mu[1]*Mu[2])*_Z^3+(-N*Mu[1]*beta[1]*beta[2]^2*h[1]*h[2]*omega[2]-N*Mu[2]*beta[1]^2*beta[2]*h[1]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]*beta[1]*h[1]+E*N*d*Mu[1]*Mu[2]*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]^2*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]*beta[2]^2*h[1]-N*alpha*Mu[1]*Mu[2]*beta[1]*h[1]-N*alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^2+(-N*Mu[1]*beta[1]*beta[2]*h[1]*omega[2]-N*Mu[1]*beta[2]^2*h[2]*omega[2]-N*Mu[2]*beta[1]^2*h[1]*omega[1]-N*Mu[2]*beta[1]*beta[2]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]+N*Mu[1]*Mu[2]*beta[1]^2+N*Mu[1]*Mu[2]*beta[2]^2-N*alpha*Mu[1]*Mu[2])*_Z-N*omega[2]*beta[2]*Mu[1]-N*omega[1]*beta[1]*Mu[2])*omega[2]*beta[2]*h[2]+RootOf(alpha*beta[2]*beta[1]*Mu[2]*Mu[1]*h[2]*h[1]*_Z^5+(alpha*Mu[1]*Mu[2]*beta[1]*h[1]+alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^4+(E*N*d*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]-N*alpha*Mu[1]*Mu[2]*beta[1]*beta[2]*h[1]*h[2]+alpha*Mu[1]*Mu[2])*_Z^3+(-N*Mu[1]*beta[1]*beta[2]^2*h[1]*h[2]*omega[2]-N*Mu[2]*beta[1]^2*beta[2]*h[1]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]*beta[1]*h[1]+E*N*d*Mu[1]*Mu[2]*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]^2*beta[2]*h[2]+N*Mu[1]*Mu[2]*beta[1]*beta[2]^2*h[1]-N*alpha*Mu[1]*Mu[2]*beta[1]*h[1]-N*alpha*Mu[1]*Mu[2]*beta[2]*h[2])*_Z^2+(-N*Mu[1]*beta[1]*beta[2]*h[1]*omega[2]-N*Mu[1]*beta[2]^2*h[2]*omega[2]-N*Mu[2]*beta[1]^2*h[1]*omega[1]-N*Mu[2]*beta[1]*beta[2]*h[2]*omega[1]+E*N*d*Mu[1]*Mu[2]+N*Mu[1]*Mu[2]*beta[1]^2+N*Mu[1]*Mu[2]*beta[2]^2-N*alpha*Mu[1]*Mu[2])*_Z-N*omega[2]*beta[2]*Mu[1]-N*omega[1]*beta[1]*Mu[2])*beta[2]*Mu[2]-omega[2])/Mu[2]]]

(1)

NULL

NULL


 

Download titikkesetimbangan.mw

1 2 Page 1 of 2