Norma Muhtar

10 Reputation

4 Badges

0 years, 294 days

MaplePrimes Activity


These are replies submitted by

@tomleslie  iam sory, i have done it in this sheet2_fit.mw
 

@tomleslie i am sory. how do I draw a graph with the values of K and r already found for t = 1...12.
2_fit.mw

@tomleslie i am sory there is a mistake for my program. i mistyped the exponential in my program. and I've fixed it. is this suitable?
can you help me to make the figure for t  t =10? thank yu very much
2_fit.mw

@tomleslie thank you for your program. maybe you cen help me if my equation is like this:
N(t) = K/(1+(((K/80)-1)*e^(-r*t))) . e is eksponensial. the program become error when use Exp. Thank for your help.
1_fit.mw1_fit.mw1_fit.mw1_fit.mw

@mmcdara . thank you for the program. i had saw the program and i confuse about rho and r. in this program you found rho = 0.0230521 , when rho=-r implies r = -0.0230521. based on the teory r should be positive because r is the population growth rate. 
let equation : 
N(t)= k/(1+((k/80)-1)e^(-r*t)))
in teory if t = 0 then N(t) = 80.
if t uproach to infinity then N(t) uproach to K or N(t) => K. but in this program obtained K = 10, far below the initial value N(t) = 80.
please advise
Thankyou Very much

@mmcdara thank you very much, you are very helpful.

@mmcdara is it the same if the kurve is

y=k/(1 + (k/80 - 1)*e^-(r*t)).

e is exp

@tomleslie 

thank you very much, i will learning it

@mmcdara 

thank you very much, i will learning it.

@Rouben Rostamian  thank you so much. i will edit my worksheet. Thank you very much for your help. 

 

NULL

restart

with(linalg); with(LinearAlgebra); with(plots); with(Student[VectorCalculus]); with(DEtools); with(plots, implicitplot); with(Student[Calculus1])

`αAP` := 2; `μA` := 1; `αlS` := 2; `μl` := .2; `βlA` := .5; `αS` := 100; `βSl` := 1; `βSA` = .25

2

 

1

 

2

 

.2

 

.5

 

100

 

1

 

`βSA` = .25

(1)

Sistem := [diff(A(t), t) = `αAP`*A(t)-`μA`*A(t), diff(l(t), t) = `αlS`*l(t)*S(t)-`μl`*l(t)-`βlA`*l(t)*A(t), diff(S(t), t) = `αS`*S(t)-`βSl`*l(t)*S(t)-.25*S(t)*A(t)]

[diff(A(t), t) = A(t), diff(l(t), t) = 2*l(t)*S(t)-.2*l(t)-.5*l(t)*A(t), diff(S(t), t) = 100*S(t)-l(t)*S(t)-.25*S(t)*A(t)]

(2)

IC := [[A(0) = 0, l(0) = 0, S(0) = 0], [A(0) = 5, l(0) = 10, S(0) = 0], [A(0) = 3, l(0) = 0, S(0) = 8], [A(0) = 0, l(0) = 5, S(0) = 10], [A(0) = 40, l(0) = 50, S(0) = 0], [A(0) = 5, l(0) = 5, S(0) = 50]]

NULL

RHS := eval(map(rhs, Sistem), [A(t) = A, l(t) = l, S(t) = S])

[A, 2*l*S-.2*l-.5*l*A, 100*S-l*S-.25*S*A]

(3)

with(plots)

NULL

p1 := fieldplot3d(RHS, A = 0 .. 50, l = 0 .. 50, S = 0 .. 50, arrows = THICK, color = red, grid = [5, 5, 5], axes = boxed)

 

p1 := fieldplot3d(RHS, A = 0 .. 50, l = 0 .. 50, S = 0 .. 50, arrows = THICK, color = red, grid = [5, 5, 5], axes = boxed)

p6 := DEplot3d(Sistem, [A(t), l(t), S(t)], t = 0 .. 50, IC, stepsize = 0.2e-2, linecolor = blue)

Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further right of 8.7375101, maxfun limit exceeded (see ?dsolve,maxfun for details)
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further right of 6.8278866, maxfun limit exceeded (see ?dsolve,maxfun for details)
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further right of 19.483533, maxfun limit exceeded (see ?dsolve,maxfun for details)
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further right of 6.6568647, maxfun limit exceeded (see ?dsolve,maxfun for details)
Warning, plot may be incomplete, the following errors(s) were issued:
   cannot evaluate the solution further right of 8.1564439, maxfun limit exceeded (see ?dsolve,maxfun for details)

 

display(p1, p6, view = [0 .. 50, 0 .. 50, 0 .. 50])

 

NULL

NULL

NULL

Download faseportrait_3dAIS.mw

@Rouben Rostamian  i have upload my sheet. Thankyou very much

Page 1 of 1