Question: How do i run and plot this code involving system of nonlinear differential equations. Also, how do i debug this code with maple codes

 

restart

with(numapprox):

with*plots:

I

 

Warning, The imaginary unit, I, has been renamed _I

 

M__h := .50; 1; beta[o] := 0.34e-1; 1; beta[1] := 0.25e-1; 1; mu[r] := 0.4e-3; 1; sigma := .7902; 1; alpha := .11; 1; psi := 0.136e-3; 1; xi := 0.5e-1; 1; gamma := .7; 1; M__c := .636; 1; mu[b] := 0.5e-2; 1; `ϖ` := .134

.50

 

0.34e-1

 

0.25e-1

 

0.4e-3

 

.7902

 

.11

 

0.136e-3

 

0.5e-1

 

.7

 

.636

 

0.5e-2

 

.134

(1)

B(0) := .50;

.50

 

.30

 

.21

 

.14

 

.70

 

.45

 

.14

(2)

ODEs := {diff(J[1](T), T) = M__h-beta[1]*psi*(J[1](T)+J[2](T))*J[7](T)-sigma*psi*beta[1]*J[4](T)*J[7](T)-mu[r]*J[1](T), diff(J[2](T), T) = beta[1]*psi*(J[1](T)+J[2](T))*J[7](T)-(alpha+xi+mu[r])*J[2](T), diff(J[3](T), T) = alpha*J[2](T)-(`ϖ`+mu[r])*J[3](T), diff(J[4](T), T) = `ϖ`*J[3](T)-(gamma+mu[r])*J[4](T), diff(J[5](T), T) = gamma*J[4](T)+sigma*psi*beta[1]*J[5](T)*J[7](T)-mu[r]*J[5](T), diff(J[6](T), T) = M__c-psi*beta[o]*J[6](T)*J[3](T)-mu[b]*J[6](T), diff(J[7](T), T) = psi*beta[o]*J[6](T)*J[3](T)-mu[b]*J[7](T)}

{diff(J[1](T), T) = .50-0.3400e-5*(J[1](T)+J[2](T))*J[7](T)-0.26866800e-5*J[4](T)*J[7](T)-0.4e-3*J[1](T), diff(J[2](T), T) = 0.3400e-5*(J[1](T)+J[2](T))*J[7](T)-.1604*J[2](T), diff(J[3](T), T) = .11*J[2](T)-.1344*J[3](T), diff(J[4](T), T) = .134*J[3](T)-.7004*J[4](T), diff(J[5](T), T) = .7*J[4](T)+0.26866800e-5*J[5](T)*J[7](T)-0.4e-3*J[5](T), diff(J[6](T), T) = .636-0.4624e-5*J[6](T)*J[3](T)-0.5e-2*J[6](T), diff(J[7](T), T) = 0.4624e-5*J[6](T)*J[3](T)-0.5e-2*J[7](T)}

(3)

ic1 := {J[1](0) = B(0), J[2](0) = C(0), J[3](0) = DD(0), J[4](0) = E(0), J[5](0) = F(0), J[6](0) = G(0), J[7](0) = H(0)};

{J[1](0) = .50, J[2](0) = .30, J[3](0) = .21, J[4](0) = .14, J[5](0) = .70, J[6](0) = .45, J[7](0) = .14}

(4)

sol1 := dsolve(`union`(ODEs, ic1), {J[1](T), J[2](T), J[3](T), J[4](T), J[5](T), J[6](T), J[7](T)}, type = numeric, output = listprocedure)