MuhammadSinan

Mr. Muhammad Sinan

5 Reputation

One Badge

2 years, 285 days
University of Swat
University
Mingora, Pakistan

Social Networks and Content at Maplesoft.com

I belong to Pakistan and a student of Mathematics at University of Swat.

MaplePrimes Activity


These are questions asked by MuhammadSinan

Please check this code and fix the error. And help me with the below two questions.

1- In the above image, I'm getting the error when I'm to plot my solutions.

2- How to get information from this system?

      The values of x_1 and corresponding y_1 and so on and errors etc

This is a downloadable link to my worksheet [Practice.mw]

Thanks!


>restart;
> ODEtools[declare]((S, L, B)(t), prime = t);
         ODEtools[declare](S(t), L(t), B(t), prime = t)
> DE1 := diff(S(t), t) = -L*S*beta-`μS`+mu;
                    d                               
            DE1 := --- S(t) = -L S beta - μS + mu
                    dt                              
> DE2 := diff(L(t), t) = L*S*beta-(gamma+mu)*S;
                  d                                  
          DE2 := --- L(t) = beta S L - (gamma + mu) S
                  dt                                 
> DE3 := diff(B(t), t) = L*gamma-gamma*mu;
                      d                           
              DE3 := --- B(t) = L gamma - gamma mu
                      dt                          
> init_conds := S(0) >= 0, L(0) >= 0, B(0) >= 0;
         init_conds := 0 <= S(0), 0 <= L(0), 0 <= B(0)
> sys := {init_conds, diff(B(t), t) = L*gamma-gamma*mu, diff(L(t), t) = L*S*beta-(gamma+mu)*S, diff(S(t), t) = -L*S*beta-`&mu;S`+mu};

> sol := dsolve(sys, numeric, parameters = [mu, beta, gamma, S(t), L(t), B(t)], method = rkf45);

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations


> sol(parameters = [mu = .234, beta = 2.345, gamma = 5.678, S(t) = 6.678, L(t) = 6.789, B(t) = 7.123]);
           sol(parameters = [mu = 0.234, beta = 2.345, gamma = 5.678, S(t) = 6.678, L(t) = 6.789, B(t) = 7.123])

Page 1 of 1