Advanced Search

Question: Help ! trying to fit assumed data values to a model system but ran into coding problems.

Question:Help ! trying to fit assumed data values to a model system but ran into coding problems.

ogunmiloro 120 Maple 18
ode syntax homework dsolve numeric differential-equations duplicate-question
October 05 2019
0


#I'm trying to fit data to a system equations with variables and parameters but couldnt run through.

``

> 

restart; with(plots); with(DEtools); with(LinearAlgebra); with(Statistics); with(CurveFitting); with(Optimization)
> 

iV := Vector([300870, 83142, 155094, 146394, 177041, 94876, 150116, 158735, 228170, 144377, 167341, 44166, 136422, 125252, 161935, 138512], datatype = float):
> 

len := Dimension(iV):
> 

tV := Vector([seq(j, j = 1 .. len)]):
> 

pp := plot(tV, iV, style = point):

````   

> 

eq1 := diff(B(T), T) = M__h-beta__1*psi*(B(T)+C(T))*H(T)-sigma*psi*beta__1*E(T)*DD(T)-mu__r*B(T):
> 

eq2 := diff(C(T), T) = beta__1*psi*(B(T)+C(T))*H(T)-(alpha+xi+mu__r)*C(T):
> 

eq3 := diff(DD(T), T) = alpha*C(T)-(varpi+xi+mu__r)*DD(T):
> 

eq4 := diff(E(T), T) = varpi*DD(T)-(Ggamma+mu__r)*E(T):
> 

eq5 := diff(F(T), T) = Ggamma*E(T)+sigma*psi*beta__1*F(T)*H(T)-mu__r*F(T):
> 

eq6 := diff(G(T), T) = M__c-psi*beta__o*G(T)*DD(T)-mu__b*G(T):
> 

eq7 := diff(H(T), T) = psi*beta__o*G(T)*DD(T)-mu__b*H(T):
> 

B0 := .50

.50

 (1)
> 

C0 := .30

.30

 (2)
> 

DD0 := .21

.21

 (3)
> 

E0 := .14

.14

 (4)
> 

F0 := .70

.70

 (5)
> 

G0 := .45

.45

 (6)
> 

H0 := .14

.14

 (7)

 

> 

sys := {eq1, eq2, eq3, eq4, eq5, eq6, eq7, B(0) = B0, C(0) = C0, DD(0) = DD0, E(0) = E0, F(0) = F0, G(0) = G0, H(0) = H0}:

 

> 

N := dsolve(sys, numeric, parameters = [M__h = .352, beta__o = 0.79e-3, beta__1 = 0.198e-2, mu__r = 0.96e-2, sigma = .7902, alpha = .65, psi = 0.2e-2, xi = 0.5e-1, Ggamma = .7, M__c = .144, mu__b = 0.68e-1, varpi = .134], output = listprocedure)

Error, (in dsolve/numeric) 'parameters' must be specified as a list of unique unassigned names

  
> 

ip := subs(N, i(t))

Error, invalid input: subs received N, which is not valid for its 1st argument

  
> 

"SSQ:=proc(`M__h`=.352, `beta__o`=0.00079, `beta__1`=0.00198, `mu__r`=0.0096, sigma=.7902, alpha=.65, psi=0.002, xi=0.5e-1, Ggamma=.7, `M__c`=.144, `mu__b`=0.068,varpi=.134 ) local j; N(parameters=[`beta__1`,alpha ]); add((ip(j)-iV[j])^(2),j=1..len) end proc:"

Error, invalid procedure parameters

"SSQ:=proc(M__h=.352, beta__o=0.00079, beta__1=0.00198, mu__r=0.0096, sigma=.7902, alpha=.65, psi=0.002, xi=0.5e-1, Ggamma=.7, M__c=.144, mu__b=0.068,varpi=.134 ) local j; N(parameters=[`beta__1`,alpha ]); add((ip(j)-iV[j])^2,j=1..len) end proc:"

  

NULL

> 

plot3d(SSQ, 0 .. 0.2e-3, 0 .. .1, axes = boxed, labels = [beta, alpha, SSQ])

Error, (in plot3d) procedure expected, as ranges contain no plotting variables

  
> 

contourplot(SSQ, 0 .. 0.2e-3, 0 .. .1, labels = [beta, alpha], contours = 40)

Error, (in plot/iplot2d/levelcurve) could not evaluate expression

  
> 

res := Minimize(SSQ, 0 .. 0.2e-3, 0 .. .1)

Error, (in Optimization:-LPSolve) invalid input: too many and/or wrong type of arguments passed to ActiveSet; first unused argument is 0 .. 0.2e-3

  
> 

res[2]

res[2]

 (8)
> 

N(parameters = convert(res[2], list)):
> 

pt := odeplot(N, [t, i(t)], 1 .. len, color = blue)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

  
> 

display(pp, pt)

Error, (in plots:-display) expecting plot structure but received: pt

  
> 

``


 

Download Fit2.mw
Please Wait...