Advanced Search

ogunmiloro

120 Reputation

6 Badges

6 years, 180 days
Contact ogunmiloro

MaplePrimes Activity


These are replies submitted by ogunmiloro

Help on the parameter estimation of a co...

ogunmiloro 120
April 23 2020
0 0

@tomleslie @Carl Love @Kitonum please these are problems I want to adress as earlier posed by @mmcdara 10007. The problems deal with parameter estimation. I intend to fit the data to the model and obtain the unknown values and also the plot of the error. Help, kindly see to this

Update on data...

ogunmiloro 120
April 22 2020
0 0

@tomleslie Thank you for the job done. However, i experienced an error when trying to open it.~

Error, (in ExcelTools:-Import) Could not open the file.
                      someData[4 .. (), 8]
                     someData[4 .. (), 22]


Also, Please i want you to take just an example of a year and perform a time series and fitting to a SIR model. That will be of great guide to me in subsequent data analysis. Thanks in Advance

> 

  restart;
> 

#
# Read the data - NB will have ti change the path
# name and file name to something appropriate for
# his/her installation
#
  someData:=ExcelTools:-Import("C:/Users/TomLeslie/Desktop/testDAta.xls"):
#
# There are many ways to "organise" this data into
# ways which are convenient to access - the following
# is just one possibility
#
  data:= table
         ( [ seq
             ( 2013+j=table
                      ( [ cases=someData[4..,6+4*j],
                          labCfd=someData[4..,7+4*j],
                          deaths=someData[4..,8+4*j],
                          CFR=someData[4..,9+4*j]
                        ]
                      ),
               j=0..4
             )
           ]
         ):
#
# Access data by year and category, where the former is
# one of 2013, 2014,....2017 and the latter is one of
# cases, labCfd, deaths, CFR
#
  data[2013][deaths];
  data[2017][cases];

_rtable[18446744074355802590]

  

_rtable[18446744074355804270]

 (1)
> 

 

> 

 


 

Download dataOrg.mw

A proper data imported from excel...

ogunmiloro 120
April 20 2020
0 0

@Carl Love thanks very much for your effort. I have imported a proper data from excel.
This data involves  diseases cases from 2012 - 2017. In the data the following terms are properly defined
(1) Susp---Suspected cases
(2)Labcfd---Laboratory confirmed of confirmation
(3) Deaths
(4)CFR-----Cases fatality rate.
Thanks i look forward to the properly time series analysisdisease_data-_2012-2017_in_regions..xlsdisease_data-_2012-2017_in_regions..xlsdisease_data-_2012-2017_in_regions..xls

awaiting your help...

ogunmiloro 120
April 18 2020
0 0

@mmcdara

please i stll expect your help towards the actualizations of the computations

System solutions repeating itself...

ogunmiloro 120
October 22 2019
0 0

#
 

> 

restart
> 

N__1 := 9;

9

  

3

  

1

  

4

  

0.4e-1

  

.4

  

.3

  

0.1e-1

  

.1

  

.2

  

.6

  

0.3e-1

  

.8

  

proc (t) options operator, arrow; N__1 end proc

  

proc (t) options operator, arrow; N__2 end proc

  

proc (t) options operator, arrow; N__3 end proc

  

proc (t) options operator, arrow; N__4 end proc

  

3

 (1)
> 

for n from 0 to m do `x__n+1` := N__1+lambda*t-d*(int(add(x__n__, n = 0 .. m), t = 0 .. t))-beta*(int(add(x__n, n = 0 .. m)*x__n*v__n, t = 0 .. t)); `w__n+1` := N__2+(1-q)*beta*(int(add(x__n*v__n, n = 0 .. m), t = 0 .. t))-e*(int(add(w__n, n = 0 .. m), t = 0 .. t))-delta*(int(add(w__n, n = 0 .. m), t = 0 .. t)); `y__n+1` := N__3+q*beta*(int(add(x__n*v__n, n = 0 .. m), t = 0 .. t))-a*(int(add(y[n], n = 0 .. m), t = 0 .. t))+delta*(int(add(w__n, n = 0 .. m), t = 0 .. t)); `v__n+1` := N__4+k*(int(add(y[n], n = 0 .. m), t = 0 .. t))-mu*(int(add(v[n], n = 0 .. m), t = 0 .. t)) end do;

9+.4*t-0.4e-1*x__n__*t-.16*x__n^2*v__n*t

  

3+0.32e-1*x__n*v__n*t-1.6*w__n*t

  

1+.128*x__n*v__n*t-.2*y[0]*t-.2*y[1]*t-.2*y[2]*t-.2*y[3]*t+1.2*w__n*t

  

4+.6*y[0]*t+.6*y[1]*t+.6*y[2]*t+.6*y[3]*t-mu*(t*v[0]+t*v[1]+t*v[2]+t*v[3])

  

9+.4*t-0.4e-1*x__n__*t-.16*x__n^2*v__n*t

  

3+0.32e-1*x__n*v__n*t-1.6*w__n*t

  

1+.128*x__n*v__n*t-.2*y[0]*t-.2*y[1]*t-.2*y[2]*t-.2*y[3]*t+1.2*w__n*t

  

4+.6*y[0]*t+.6*y[1]*t+.6*y[2]*t+.6*y[3]*t-mu*(t*v[0]+t*v[1]+t*v[2]+t*v[3])

  

9+.4*t-0.4e-1*x__n__*t-.16*x__n^2*v__n*t

  

3+0.32e-1*x__n*v__n*t-1.6*w__n*t

  

1+.128*x__n*v__n*t-.2*y[0]*t-.2*y[1]*t-.2*y[2]*t-.2*y[3]*t+1.2*w__n*t

  

4+.6*y[0]*t+.6*y[1]*t+.6*y[2]*t+.6*y[3]*t-mu*(t*v[0]+t*v[1]+t*v[2]+t*v[3])

  

9+.4*t-0.4e-1*x__n__*t-.16*x__n^2*v__n*t

  

3+0.32e-1*x__n*v__n*t-1.6*w__n*t

  

1+.128*x__n*v__n*t-.2*y[0]*t-.2*y[1]*t-.2*y[2]*t-.2*y[3]*t+1.2*w__n*t

  

4+.6*y[0]*t+.6*y[1]*t+.6*y[2]*t+.6*y[3]*t-mu*(t*v[0]+t*v[1]+t*v[2]+t*v[3])

 (2)
> 

NULL
> 

NULL


 

Download MPGH_sim.mw

@Carl Love thanks very much. I corrected that, but i still have the system solutions repeating itself. This #is what i got.

This is what i observe and still need so...

ogunmiloro 120
October 16 2019
0 0

@tomleslie thank you very much for the job well done. This is the problem i encounter
(1) My maple version is not running the graphs as it ought to be, but it ran through in yours

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/someODEs2.mw .
 

Download someODEs2.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/someODEs2.mw .
 

Download someODEs2.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/someODEs2.mw .
 

Download someODEs2.mw

Maple Worksheet - Error


Maple Worksheet - Error

.
 

Download someODEs2.mw

 

 

Download someODEs2.mw

 

 

 

 

. you will see it in the code i uploaded. i mean it was giving me something like
PLOT( )
PLOT( ) and so on.

(2) i want to obtain the numerical values  from the series code you wrote as

(3)Also, if the answers are thesame i think the combined graphs should fit together
i mean the code you wrote as
display([v1, a1]); display([v2, a2]); display([v4, a4])

@Christian Wolinski  I initialized ...

ogunmiloro 120
October 15 2019
0 0

@Christian Wolinski 
I initialized at B(0) := .50; C(0) := .30; DD(0) := .21; E(0) := .14; F(0) := .70; G(0) := .45; H(0) := .14 but it wasnt running smoothly
 

The first two points raised by you are v...

ogunmiloro 120
October 06 2019
0 0

@mmcdara thanks for your help. The first two points raised by you are very correct. B(T) is the priori model of the empirical observations iV. So I want iV to be considered as a observation for C(T) or DD(T). The model is resumed by 16 empirical observations by 7 matrix dimensions. Also, I want to know if transmissibility parameters like beta1 and betao values can be fitted to the experimental observations iV. You can please go forward with the coding/computations. Thanks.

There is no relationship between i(T) a...

ogunmiloro 120
October 05 2019
0 0

@mmcdara
There is no relationshipt between i(t) and B(T), C(T),....H(T)....it is wrong. i intend to write either B(T) or....H(T) where i wrongly wrote i(T).
 pt := odeplot(N, [t, i(t)], 1 .. len, color = blue)
My intention is to fit the data values with each state variables and see the behavior on the two on the same graph

Fitting real life assumed data and compa...

ogunmiloro 120
October 05 2019
0 0

@mmcdara 777 thanks a lot for your expalantions and the code correction.
These are my intentions;
(1) The system of equations describes an epidemic transmission which i got an assumed data values for infected individuals per year in different localities earlier listed in the begining of the code. i intend to fit this data to the model equations which still retains its parameters and variable values.
I want to fit this data values with each variables to see the relationship/behavior between real life data and parameters/variable values involved in the model system equations.
Or if you have another fitting method you can help with that and i will gladly apprreciate.
Please i need help on that.
Thanks
 

power series code not working for my ver...

ogunmiloro 120
October 04 2019
0 0

@tomleslie 
The code to obtain the power series solution to this system is not working, Please i need help on a new code that will support my maple version


 

> 

restart; with(plots); interface(version); _local(gamma)

`Standard Worksheet Interface, Maple 18.00, Windows 8, February 10 2014 Build ID 922027`

 (1)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*(B(T)+C(T))*H(T)-sigma*psi*beta__1*E(T)*H(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*(B(T)+C(T))*H(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(`ϖ`+xi+mu__r)*DD(T), diff(E(T), T) = `ϖ`*DD(T)-(gamma+mu__r)*E(T), diff(F(T), T) = gamma*E(T)+sigma*psi*beta__1*F(T)*H(T)-mu__r*F(T), diff(G(T), T) = M__c-psi*beta__o*G(T)*DD(T)-mu__b*G(T), diff(H(T), T) = psi*beta__o*G(T)*DD(T)-mu__b*H(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 

 
> 

interface(rtablesize = 22); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

M := Matrix(22, 8, {(1, 1) = T, (1, 2) = B(T), (1, 3) = C(T), (1, 4) = DD(T), (1, 5) = E(T), (1, 6) = F(T), (1, 7) = G(T), (1, 8) = H(T), (2, 1) = 0., (2, 2) = 50., (2, 3) = 30., (2, 4) = 21., (2, 5) = 14., (2, 6) = 70., (2, 7) = 45., (2, 8) = 14., (3, 1) = 1., (3, 2) = 50.4757829913694, (3, 3) = 25.5575079474562, (3, 4) = 20.2413350094118, (3, 5) = 8.93476021674514, (3, 6) = 77.8009301456023, (3, 7) = 45.4056689838066, (3, 8) = 13.9344799352524, (4, 1) = 2., (4, 2) = 50.9517042861457, (4, 3) = 21.7732057953036, (4, 4) = 19.1991185729657, (4, 5) = 6.33111424301372, (4, 6) = 83.0173704221051, (4, 7) = 45.8094678095317, (4, 8) = 13.8691335287751, (5, 1) = 3., (5, 2) = 51.4276669563117, (5, 3) = 18.5495796746541, (5, 4) = 17.9820385233544, (5, 5) = 4.92813671736454, (5, 6) = 86.8798586692581, (5, 7) = 46.2114574379362, (5, 8) = 13.8039082811349, (6, 1) = 4., (6, 2) = 51.9036130320714, (6, 3) = 15.8035626969457, (6, 4) = 16.6714516583213, (6, 5) = 4.10887104570361, (6, 6) = 89.9876659437646, (6, 7) = 46.6116795393824, (6, 8) = 13.7387709410969, (7, 1) = 5., (7, 2) = 52.3795051400209, (7, 3) = 13.4643946140523, (7, 4) = 15.3273237666627, (7, 5) = 3.57371565743505, (7, 6) = 92.6321391296507, (7, 7) = 47.0101611197031, (7, 8) = 13.6737028799647, (8, 1) = 6., (8, 2) = 52.8553172609914, (8, 3) = 11.4717985735174, (8, 4) = 13.9929712055072, (8, 5) = 3.17880170203220, (8, 6) = 94.9557749671457, (8, 7) = 47.4069181944613, (8, 8) = 13.6086964175305, (9, 1) = 7., (9, 2) = 53.3310300135796, (9, 3) = 9.77442795393849, (9, 4) = 12.6988352778328, (9, 5) = 2.85616909435746, (9, 6) = 97.0300691802980, (9, 7) = 47.8019586919789, (9, 8) = 13.5437519192551, (10, 1) = 8., (10, 2) = 53.8066281916241, (10, 3) = 8.32854334230931, (10, 4) = 11.4654771278926, (10, 5) = 2.57446425033359, (10, 6) = 98.8932984213762, (10, 7) = 48.1952847329760, (10, 8) = 13.4788755168356, (11, 1) = 9., (11, 2) = 54.2820994312045, (11, 3) = 7.09688550193379, (11, 4) = 10.3059458447231, (11, 5) = 2.31968445942390, (11, 6) = 100.568614896467, (11, 7) = 48.5868944080828, (11, 8) = 13.4140773308998, (12, 1) = 10., (12, 2) = 54.7574334499874, (12, 3) = 6.04771539478717, (12, 4) = 9.22764331200851, (12, 5) = 2.08581976619023, (12, 6) = 102.072526790591, (12, 7) = 48.9767831521078, (12, 8) = 13.3493700969435, (13, 1) = 11., (13, 2) = 55.2326215818715, (13, 3) = 5.15399607440913, (13, 4) = 8.23378670861787, (13, 5) = 1.87036441689728, (13, 6) = 103.418731306890, (13, 7) = 49.3649447962597, (13, 8) = 13.2847681133127, (14, 1) = 12., (14, 2) = 55.7076564691196, (14, 3) = 4.39269644819759, (14, 4) = 7.32454792482105, (14, 5) = 1.67221333782958, (14, 6) = 104.619737022156, (14, 7) = 49.7513723627469, (14, 8) = 13.2202864468074, (15, 1) = 13., (15, 2) = 56.1825318417621, (15, 3) = 3.74419753068560, (15, 4) = 6.49793765842352, (15, 5) = 1.49071216478514, (15, 6) = 105.687467337294, (15, 7) = 50.1360586566244, (15, 8) = 13.1559403410368, (16, 1) = 14., (16, 2) = 56.6572423491652, (16, 3) = 3.19178693040782, (16, 4) = 5.75048524495814, (16, 5) = 1.32526003144654, (16, 6) = 106.633415749484, (16, 7) = 50.5189966968304, (16, 8) = 13.0917447855840, (17, 1) = 15., (17, 2) = 57.1317834251097, (17, 3) = 2.72122859254014, (17, 4) = 5.07775680378739, (17, 5) = 1.17516896160884, (17, 6) = 107.468624579589, (17, 7) = 50.9001800213339, (17, 8) = 13.0277142110589, (18, 1) = 16., (18, 2) = 57.6061511758188, (18, 3) = 2.32039527095494, (18, 4) = 4.47474747095627, (18, 5) = 1.03963680345053, (18, 6) = 108.203610436062, (18, 7) = 51.2796028953429, (18, 8) = 12.9638622810886, (19, 1) = 17., (19, 2) = 58.0803422857694, (19, 3) = 1.97895651531399, (19, 4) = 3.93617229075042, (19, 5) = .917766153838763, (19, 6) = 108.848289685285, (19, 7) = 51.6572604433541, (19, 8) = 12.9002017604666, (20, 1) = 18., (20, 2) = 58.5543539371233, (20, 3) = 1.68811238325572, (20, 4) = 3.45667982032798, (20, 5) = .808597352505364, (20, 6) = 109.411923234287, (20, 7) = 52.0331487242844, (20, 8) = 12.8367444402178, (21, 1) = 19., (21, 2) = 59.0281837405625, (21, 3) = 1.44036675537223, (21, 4) = 3.03100499830172, (21, 5) = .711142050392071, (21, 6) = 109.903085135445, (21, 7) = 52.4072647634383, (21, 8) = 12.7735011058283, (22, 1) = 20., (22, 2) = 59.5018296756537, (22, 3) = 1.22933413066306, (22, 4) = 2.65407535007048, (22, 5) = .624412016612247, (22, 6) = 110.329653167156, (22, 7) = 52.7796065526334, (22, 8) = 12.7104815373142})

 (2)
> 

Order := 12; ans2 := `~`[convert](`~`[evalf](dsolve*({ODEs, bcs}, indets([ODEs], function(name)), 'series')), polynom); for j in indets([ODEs], function(name)) do display*([odeplot*(ans, [T, j], T = 0 .. 10, thickness = 2, color = red), plot*(eval(j, ans2), T = 0 .. 10, thickness = 2, color = blue)], title = cat(convert(j, string), " (Numerical and PowerSeries compared)"), titlefont = [tims, bold, 20]) end do

dsolve*({B(0) = 50., C(0) = 30., DD(0) = 21., E(0) = 14., F(0) = 70., G(0) = 45., H(0) = 14., diff(B(T), T) = .50-0.3400e-5*(B(T)+C(T))*H(T)-0.26866800e-5*E(T)*H(T)-0.4e-3*B(T), diff(C(T), T) = 0.3400e-5*(B(T)+C(T))*H(T)-.1604*C(T), diff(DD(T), T) = .11*C(T)-.1844*DD(T), diff(E(T), T) = .134*DD(T)-.7004*E(T), diff(F(T), T) = .7*E(T)+0.26866800e-5*F(T)*H(T)-0.4e-3*F(T), diff(G(T), T) = .636-0.4624e-5*G(T)*DD(T)-0.5e-2*G(T), diff(H(T), T) = 0.4624e-5*G(T)*DD(T)-0.5e-2*H(T)}, {B(T), C(T), DD(T), E(T), F(T), G(T), H(T)}, series)

  

Error, invalid input: eval received dsolve*({B(0) = 50., C(0) = 30., DD(0) = 21., E(0) = 14., F(0) = 70., G(0) = 45., H(0) = 14., diff(B(T), T) = .50-0.3400e-5*(B(T)+C(T))*H(T)-0.26866800e-5*E(T)*H(T)-0.4e-3*B(T), diff(C(T), T) = 0.3400e-5*(B(T)+C(T))*H(T)-.1604*C(T), diff(DD(T), T) = .11*C(T)-.1844*DD(T), diff(E(T), T) = .134*DD(T)-.7004*E(T), diff(F(T), T) = .7*E(T)+0.26866800e-5*F(T)*H(T)-0.4e-3*F(T), diff(G(T), T) = .636-0.4624e-5*G(T)*DD(T)-0.5e-2*G(T), diff(H(T), T) = 0.4624e-5*G(T)*DD(T)-0.5e-2*H(T)}, {B(T), C(T), DD(T), E(T), F(T), G(T), H(T)}, series), which is not valid for its 2nd argument, eqns

  
> 

``


 

Download Phd11.mw

complex or singularity error...

ogunmiloro 120
October 03 2019
0 0

@tomleslie good day, thanks for your attention. What can cause complex or singularity error in this case as shown below
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*(B(T)+C(T))*V(T)-sigma*psi*beta__1*E(T)*V(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*(B(T)+C(T))*V(T)-(alpha+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(`ϖ`+xi+mu__r)*DD(T), diff(E(T), T) = `ϖ`*DD(T)-(gamma+mu__r)*E(T), diff(F(T), T) = gamma*E(T)-sigma*psi*beta__1*E(T)*G(T)-mu__r*F(T), diff(G(T), T) = M[c]-psi*beta__o*F(T)*C(T)-mu__b*G(T), diff(V(T), T) = psi*beta__o*F(T)*C(T)-mu__b*V(T):
> 

ans := dsolve([ODEs, bcs], numeric):

Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)

  

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

  

 
> 

``


 

Download MPL1.mw

Variation of parameters...

ogunmiloro 120
September 27 2019
0 0

@tomleslie Thanks once again for your good work and special attention to this program.
My question is
(1) You have varied the parameter "xi" at three different values. I checked the code for alteration so that i can also vary other parameters like "Mh, psi, mu" etc. but i couldnt get it.
How do i vary other parameters.

Salient points raised on the maple codes...

ogunmiloro 120
September 26 2019
0 0

@tomleslie thank you very much for taking your time to attend to my problems, i'm very grateful. Here are some points i need to raise;
(1)A particular code you wrote as
interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])
is not running on my maple sheet. I dont know if my version cant support that or better still, another code that can give the matrix is appreciated.
(2) The code i wrote for the series is the differential transformation method(DTM) series applied to model equations. My intentions were to see wether the series solution for DTM will agree with the numerical values/solutions using RK4
(3) Also, i need codes for parametric plots, e.g., for Mh, beta1 etc.
Once again, I appreciate @tomleslie
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw
 

> 

restart; with(plots); _local(gamma)
> 

M__h := .50:
> 

ODEs := diff(B(T), T) = M__h-beta__1*psi*B(T)*G(T)-mu__r*B(T), diff(C(T), T) = beta__1*psi*B(T)*G(T)-sigma*psi*beta__1*E(T)*DD(T)-(alpha+xi+mu__r)*C(T), diff(DD(T), T) = alpha*C(T)-(gamma+mu__r)*DD(T), diff(E(T), T) = gamma*DD(T)+sigma*psi*beta__1*E(T)*G(T)-mu__r*E(T), diff(F(T), T) = M__c-psi*beta__o*F(T)*C(T)-mu__b*F(T), diff(G(T), T) = psi*beta__o*F(T)*C(T)-mu__b*G(T):
> 

ans := dsolve([ODEs, bcs], numeric):

 

 

 

 

 

 
> 

interface(rtablesize = [22, 10]); M := Matrix([`~`[lhs](ans(0)), seq(`~`[rhs](ans(j)), j = 0 .. 20)])

Error, (in interface) rtablesize must be a positive integer, or infinity

  

M := Vector(4, {(1) = ` 22 x 7 `*Matrix, (2) = `Data Type: `*anything, (3) = `Storage: `*rectangular, (4) = `Order: `*Fortran_order})

 (1)
> 

#
# Generate power series solutions of order 12
# (Maple default is 6) for all ODEs
#
  Order:=12:
  ans2:= convert~
         ( evalf~
           ( dsolve
             ( {ODEs, bcs},
               indets([ODEs], function(name)),
               'series'
             )
           ),
           polynom
         );
#
# Plot the "numeric" solution obtained earlier
# and the equivalent power series solution
# obtaine above on th same graph.
#
# Note that the range of the plots is reduced,
# because even with 12-th order polynomials,
# the power series soluton and the numerical
# solution start to deviate for T>~6
#
  for j in indets([ODEs], function(name)) do
      display
      ( [ odeplot
          ( ans,
            [T, j],
            T = 0 .. 10,
            thickness = 4, style = point,
            color = red
          ),
          plot
          ( eval(j, ans2),
            T = 0 .. 10,
            thickness = 2, style = point,
            color=blue
          )
        ],
        title = cat(convert(j,string), " (Numerical and PowerSeries compared)"),
        titlefont = [tims, bold, 20]
      )
   end do

{B(T) = .5000000000+.4997992350*T-0.1003402817e-3*T^2+0.1470042520e-7*T^3-0.4186689729e-11*T^4+0.1797722185e-13*T^5-0.7652371692e-15*T^6+0.2923183302e-16*T^7-0.8844182475e-18*T^8+0.2215403975e-19*T^9-0.4589870880e-21*T^10+0.5389509965e-23*T^11, C(T) = .3000000000-0.4811931399e-1*T+0.3859529419e-2*T^2-0.2063396439e-3*T^3+0.8267055682e-5*T^4-0.2631388485e-6*T^5+0.6554216926e-8*T^6-0.5551714044e-10*T^7-0.1520371564e-10*T^8+0.2777703195e-11*T^9-0.3923000816e-12*T^10+0.4970214040e-13*T^11, DD(T) = .2100000000-.1140840000*T+0.3730565453e-1*T^2-0.8568110732e-2*T^3+0.1494601849e-2*T^4-0.2091819518e-3*T^5+0.2441368229e-4*T^6-0.2442660302e-5*T^7+0.2138541461e-6*T^8-0.1664279070e-7*T^9+0.1165691616e-8*T^10-0.7422668732e-10*T^11, E(T) = .1400000000+.1469441693*T-0.3995870043e-1*T^2+0.8709964152e-2*T^3-0.1500287681e-2*T^4+0.2093639084e-3*T^5-0.2441847491e-4*T^6+0.2442759169e-5*T^7-0.2138545263e-6*T^8+0.1664257448e-7*T^9-0.1165658906e-8*T^10+0.7422262550e-10*T^11, F(T) = .7000000000+.6324990290*T-0.1581608397e-2*T^2+0.2679492283e-5*T^3-0.6093270788e-8*T^4+0.1272029637e-9*T^5-0.4253521584e-11*T^6+0.1189725836e-12*T^7-0.2702521878e-14*T^8+0.3074973933e-16*T^9+0.3481499101e-17*T^10-0.6346661500e-18*T^11, G(T) = .4500000000-0.2249028960e-2*T+0.5983397429e-5*T^2-0.5345061622e-7*T^3+0.2810718704e-8*T^4-0.1239204116e-9*T^5+0.4250786124e-11*T^6-0.1189706297e-12*T^7+0.2702520657e-14*T^8-0.3074973865e-16*T^9-0.3481499102e-17*T^10+0.6346661500e-18*T^11}

  

Warning, inserted missing semicolon at end of statement

  

 

 

 

 

 

 
> 

``


 

Download MAPQ12.mw

 

 

 

 

 

thank you very much @tomleslie, the nume...

ogunmiloro 120
September 26 2019
0 0

thank you very much @tomleslie, the numerical values of the dependent variable values you sent, how can it agree with the values of this series solution below. its only the values of initial conditions that are the same in the matrix values you sent and the one i send below. Thanks

 

 

 

 

> 

restart:

I

  

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

  

100+.4447000000*T[i]+.2498768924*T[i]^2+.1666205504*T[i]^3+.1249769369*T[i]^4+0.9998616086e-1*T[i]^5+0.8332410682e-1*T[i]^6+0.7142198089e-1*T[i]^7+0.6249505700e-1*T[i]^8

  

90-14.42529422*T[i]+1.155030935*T[i]^2-0.6143642353e-1*T[i]^3+0.2421975133e-2*T[i]^4-0.6730874690e-4*T[i]^5+0.3600667425e-5*T[i]^6+0.1822546276e-5*T[i]^7+0.1328718119e-5*T[i]^8

  

45-21.6180*T[i]+6.777232420*T[i]^2-1.539906728*T[i]^3+.2679481662*T[i]^4-0.3748089566e-1*T[i]^5+0.4374035893e-2*T[i]^6-0.4375969523e-3*T[i]^7+0.3833667319e-4*T[i]^8

  

38+31.48939422*T[i]-7.570686050*T[i]^2+1.582057225*T[i]^3-.2695938568*T[i]^4+0.3752776352e-1*T[i]^5-0.4374513525e-2*T[i]^6+0.4375776774e-3*T[i]^7-0.3830497045e-4*T[i]^8

  

70+.256868800*T[i]+.3173129456*T[i]^2+.2114336177*T[i]^3+.1587169689*T[i]^4+.1270300292*T[i]^5+.1058866358*T[i]^6+0.9077614676e-1*T[i]^7+0.7943924211e-1*T[i]^8

  

45+.195868800*T[i]+0.4447896434e-3*T[i]^2-0.3678604877e-4*T[i]^3-0.1878503860e-4*T[i]^4-0.1127262504e-4*T[i]^5-0.7515283005e-5*T[i]^6-0.5368149890e-5*T[i]^7-0.4026158289e-5*T[i]^8

  

 
> 

> 

``

 

 

Download MAPQQ.mw

 

1 2 3 4 5 6 Page 5 of 6