Advanced Search

Question: why my solutions lost?

Question:why my solutions lost?

Yee Voon 20 Maple
solve
December 26 2017
1


I am trying to find for the equilibrium but why my solutions lost?

> 

restart
> 

interface(imaginaryunit = j);

I

 (1)
> 

unprotect(Pi)
> 

lambda := k*tau*(C*Upsilon+I)/N;

k*tau*(C*Upsilon+I)/N

 (2)
> 

eqn1 := (1-p)*Pi+phi*V+delta*R-(mu+lambda+vartheta)*S;

(1-p)*Pi+phi*V+delta*R-(mu+k*tau*(C*Upsilon+I)/N+vartheta)*S

 (3)
> 

eqn2 := p*Pi+vartheta*S-(epsilon*lambda+mu+phi)*V;

p*Pi+vartheta*S-(epsilon*k*tau*(C*Upsilon+I)/N+mu+phi)*V

 (4)
> 

eqn3 := rho*lambda*S+rho*epsilon*lambda*V+(1-q)*eta*I-(mu+beta+chi)*C;

rho*k*tau*(C*Upsilon+I)*S/N+rho*epsilon*k*tau*(C*Upsilon+I)*V/N+(1-q)*eta*I-(mu+beta+chi)*C

 (5)
> 

eqn4 := (1-rho)*lambda*S+(1-rho)*epsilon*lambda*V+chi*C-(mu+alpha+eta)*I;

(1-rho)*k*tau*(C*Upsilon+I)*S/N+(1-rho)*epsilon*k*tau*(C*Upsilon+I)*V/N+chi*C-(mu+alpha+eta)*I

 (6)
> 

eqn5 := beta*C+q*eta*I-(mu+delta)*R;

beta*C+q*eta*I-(mu+delta)*R

 (7)
> 

mu := 0.1e-1;

0.1e-1

  

116.1

  

0.8e-2

  

0.25e-2

  

0.2e-2

  

0.5e-1

  

0.115e-1

  

0.598e-2

  

.5

  

.2

  

.1

  

0.57e-2

  

.2

  

11610

 (8)
> 

Equilibria := solve({eqn1 = 0, eqn2 = 0, eqn3 = 0, eqn4 = 0, eqn5 = 0}, {C, I, R, S, V});

Warning, solutions may have been lost

  

{C = 0., I = 0., R = 0., S = 5946.585366, V = 5663.414634}

 (9)
> 

``


 

Download Equilibria.mw
 

> 

restart
> 

interface(imaginaryunit = j);

I

 (1)
> 

unprotect(Pi)
> 

lambda := k*tau*(C*Upsilon+I)/N;

k*tau*(C*Upsilon+I)/N

 (2)
> 

eqn1 := (1-p)*Pi+phi*V+delta*R-(mu+lambda+vartheta)*S;

(1-p)*Pi+phi*V+delta*R-(mu+k*tau*(C*Upsilon+I)/N+vartheta)*S

 (3)
> 

eqn2 := p*Pi+vartheta*S-(epsilon*lambda+mu+phi)*V;

p*Pi+vartheta*S-(epsilon*k*tau*(C*Upsilon+I)/N+mu+phi)*V

 (4)
> 

eqn3 := rho*lambda*S+rho*epsilon*lambda*V+(1-q)*eta*I-(mu+beta+chi)*C;

rho*k*tau*(C*Upsilon+I)*S/N+rho*epsilon*k*tau*(C*Upsilon+I)*V/N+(1-q)*eta*I-(mu+beta+chi)*C

 (5)
> 

eqn4 := (1-rho)*lambda*S+(1-rho)*epsilon*lambda*V+chi*C-(mu+alpha+eta)*I;

(1-rho)*k*tau*(C*Upsilon+I)*S/N+(1-rho)*epsilon*k*tau*(C*Upsilon+I)*V/N+chi*C-(mu+alpha+eta)*I

 (6)
> 

eqn5 := beta*C+q*eta*I-(mu+delta)*R;

beta*C+q*eta*I-(mu+delta)*R

 (7)
> 

mu := 0.1e-1;

0.1e-1

  

116.1

  

0.8e-2

  

0.25e-2

  

0.2e-2

  

0.5e-1

  

0.115e-1

  

0.598e-2

  

.5

  

.2

  

.1

  

0.57e-2

  

.2

  

11610

 (8)
> 

Equilibria := solve({eqn1 = 0, eqn2 = 0, eqn3 = 0, eqn4 = 0, eqn5 = 0}, {C, I, R, S, V});

Warning, solutions may have been lost

  

{C = 0., I = 0., R = 0., S = 5946.585366, V = 5663.414634}

 (9)
> 

``


 

Download Equilibria.mw

 
Please Wait...