Advanced Search

Pemudahijrah01

25 Reputation

4 Badges

5 years, 20 days
Contact Pemudahijrah01

MaplePrimes Activity


These are replies submitted by Pemudahijrah01

Okay.. i will try it.. thanks a lot ...

Pemudahijrah01 25
June 08 2020
0 0

Okay.. i will try it.. thanks a lot @Carl Love 

Sorry, this is my worksheet..thank you ...

Pemudahijrah01 25
June 08 2020
0 0

Sorry, this is my worksheet..thank you

restart; with(DEtools); with(plots); with(linalg); alpha := .75; omega[1] := .15; omega[2] := .20; N := .8; beta[1] := .65; beta[2] := .70; Mu[1] := .75; Mu[2] := .74; h[1] := 0.5e-2; h[2] := 0.4e-2; d := 0.5e-1; E := .35

sys := [(D(x))(t) = alpha*x(t)*(1-x(t)/N)-beta[1]*sqrt(x(t))*y1(t)/(1+h[1]*beta[1]*sqrt(x(t)))-beta[2]*sqrt(x(t))*y2(t)/(1+h[2]*beta[2]*sqrt(x(t)))-d*E*x(t), (D(y1))(t) = -omega[1]*y1(t)+Mu[1]*beta[1]*sqrt(x(t))*y1(t)*(1-y1(t)/(beta[1]*sqrt(x(t))))/(1+h[1]*beta[1]*sqrt(x(t))), (D(y2))(t) = -omega[2]*y2(t)+Mu[2]*beta[2]*sqrt(x(t))*y2(t)*(1-y2(t)/(beta[2]*sqrt(x(t))))/(1+h[2]*beta[2]*sqrt(x(t)))]; RHS := eval(map(rhs, sys), [x(t) = x, y1(t) = y1, y2(t) = y2])

[.75*x*(1-1.250000000*x)-.65*x^(1/2)*y1/(1+0.325e-2*x^(1/2))-.70*x^(1/2)*y2/(1+0.280e-2*x^(1/2))-0.175e-1*x, -.15*y1+.4875*x^(1/2)*y1*(1-1.538461538*y1/x^(1/2))/(1+0.325e-2*x^(1/2)), -.20*y2+.1050*x^(1/2)*y2*(1-1.428571429*y2/x^(1/2))/(1+0.280e-2*x^(1/2))]

 (1)

plots:-fieldplot3d(sys, x = 0 .. 1, y1 = () .. 1, 0*y2 = 0 .. 1, grid = [5, 5, 5], axes = boxed)

``


Download phaseportrait.mw

1 2 Page 2 of 2