Advanced Search

Question: Query: Same Equation Giving Two Different Solutions

Question:Query: Same Equation Giving Two Different Solutions

Andiguys 85 Maple , MaplePrimes
syntax error
7 hours ago
0

I am getting two different solutions for the same set of equations. Notice that G1 is identical to G3 and G2 is identical to G4. Why is this happening? Is there a syntax issue causing this, and how can I fix it?


 

> 

restart
> 

with(Optimization); with(plots); with(LinearAlgebra)
> 

_local(Pi)

Pi

 (1)
> 

DATA := [delta = .6, sigma = .4, a = 0.5e-2, d = 0.2e-1, g = 0.3e-1, Ce = 0.5e-2, c = 0.2e-1, Cn = .25, Crm = .1, Cr = 0.7e-1, rho0 = .5, i2 = 0.5e-1, Pr = .6, t = 0.5e-1, alpha = .9]

[delta = .6, sigma = .4, a = 0.5e-2, d = 0.2e-1, g = 0.3e-1, Ce = 0.5e-2, c = 0.2e-1, Cn = .25, Crm = .1, Cr = 0.7e-1, rho0 = .5, i2 = 0.5e-1, Pr = .6, t = 0.5e-1, alpha = .9]

 (2)
> 

``
> 

``
> 

NULL
> 

subs(DATA, (2*Cr*Pr*alpha*d*rho0^2-2*Cr*alpha*d*delta*rho0^2+Ce*alpha*d*delta*rho0+2*Cr*alpha*d*rho0^2-Crm*alpha*d*delta*rho0+Pr*alpha*d*delta*rho0-alpha*c*d*delta*rho0-alpha*delta*g*i2*rho0-Ce*alpha*d*rho0-Ce*d*delta*rho0+Crm*alpha*d*rho0-Pr*alpha*d*rho0+4*Pr*alpha*rho0^2+2*a*alpha*delta*rho0+alpha*c*d*rho0-alpha*c*delta*rho0-4*alpha*delta*rho0^2+alpha*g*i2*rho0+Ce*d*rho0-Cn*d*delta-Pr*d*delta-2*a*alpha*rho0+alpha*c*rho0+4*alpha*rho0^2+d*delta^2+Cn*d+Pr*d-2*d*delta+d)/(2*(Cr*alpha*d*rho0^2+2*alpha*rho0^2-d*delta+d))); subs(DATA, -(Ce*Cr*alpha*d^2*rho0^2+Cr*Crm*alpha*d^2*rho0^2-Cr*Pr*alpha*d^2*rho0^2-Cr*alpha*c*d^2*rho0^2-Cr*alpha*d*g*i2*rho0^2+Ce*Cr*d^2*rho0^2+Cr*alpha*c*d*rho0^2+2*Ce*alpha*d*rho0^2+Cn*Cr*d^2*rho0-Cr*Pr*d^2*rho0+Cr*d^2*delta*rho0+2*Crm*alpha*d*rho0^2-2*Pr*alpha*d*rho0^2-2*alpha*c*d*rho0^2-2*alpha*g*i2*rho0^2-2*Ce*d^2*delta+2*Ce*d*rho0^2-Cr*d^2*rho0+2*alpha*c*rho0^2+2*c*d^2*delta+2*d*delta*g*i2+2*Ce*d^2+2*Cn*d*rho0-2*Pr*d*rho0-2*a*d*delta-2*c*d^2+2*d*delta*rho0-2*d*g*i2+2*a*d-2*d*rho0)/(2*d*(Cr*alpha*d*rho0^2+2*alpha*rho0^2-d*delta+d))); subs(DATA, (((g*i2-a)*Cr-c-(Ce*Cr*alpha*d^2*rho0^2+Cr*Crm*alpha*d^2*rho0^2-Cr*Pr*alpha*d^2*rho0^2-Cr*alpha*c*d^2*rho0^2-Cr*alpha*d*g*i2*rho0^2+Ce*Cr*d^2*rho0^2+Cr*alpha*c*d*rho0^2+2*Ce*alpha*d*rho0^2+Cn*Cr*d^2*rho0-Cr*Pr*d^2*rho0+Cr*d^2*delta*rho0+2*Crm*alpha*d*rho0^2-2*Pr*alpha*d*rho0^2-2*alpha*c*d*rho0^2-2*alpha*g*i2*rho0^2-2*Ce*d^2*delta+2*Ce*d*rho0^2-Cr*d^2*rho0+2*alpha*c*rho0^2+2*c*d^2*delta+2*d*delta*g*i2+2*Ce*d^2+2*Cn*d*rho0-2*Pr*d*rho0-2*a*d*delta-2*c*d^2+2*d*delta*rho0-2*d*g*i2+2*a*d-2*d*rho0)/(2*d*(Cr*alpha*d*rho0^2+2*alpha*rho0^2-d*delta+d))+Ce)*d+g*i2-a)/(d*(Cr*d+2)))

.9938295713

  

.6118434920

  

.2106517897

 (3)
> 

G08 := subs(Pn = (2*Cr*Pr*alpha*d*rho0^2-2*Cr*alpha*d*delta*rho0^2+Ce*alpha*d*delta*rho0+2*Cr*alpha*d*rho0^2-Crm*alpha*d*delta*rho0+Pr*alpha*d*delta*rho0-alpha*c*d*delta*rho0-alpha*delta*g*i2*rho0-Ce*alpha*d*rho0-Ce*d*delta*rho0+Crm*alpha*d*rho0-Pr*alpha*d*rho0+4*Pr*alpha*rho0^2+2*a*alpha*delta*rho0+alpha*c*d*rho0-alpha*c*delta*rho0-4*alpha*delta*rho0^2+alpha*g*i2*rho0+Ce*d*rho0-Cn*d*delta-Pr*d*delta-2*a*alpha*rho0+alpha*c*rho0+4*alpha*rho0^2+d*delta^2+Cn*d+Pr*d-2*d*delta+d)/(2*(Cr*alpha*d*rho0^2+2*alpha*rho0^2-d*delta+d)), 1-(Pn-Pr)/(1-delta)); G09 := subs([i1 = -(Ce*Cr*alpha*d^2*rho0^2+Cr*Crm*alpha*d^2*rho0^2-Cr*Pr*alpha*d^2*rho0^2-Cr*alpha*c*d^2*rho0^2-Cr*alpha*d*g*i2*rho0^2+Ce*Cr*d^2*rho0^2+Cr*alpha*c*d*rho0^2+2*Ce*alpha*d*rho0^2+Cn*Cr*d^2*rho0-Cr*Pr*d^2*rho0+Cr*d^2*delta*rho0+2*Crm*alpha*d*rho0^2-2*Pr*alpha*d*rho0^2-2*alpha*c*d*rho0^2-2*alpha*g*i2*rho0^2-2*Ce*d^2*delta+2*Ce*d*rho0^2-Cr*d^2*rho0+2*alpha*c*rho0^2+2*c*d^2*delta+2*d*delta*g*i2+2*Ce*d^2+2*Cn*d*rho0-2*Pr*d*rho0-2*a*d*delta-2*c*d^2+2*d*delta*rho0-2*d*g*i2+2*a*d-2*d*rho0)/(2*d*(Cr*alpha*d*rho0^2+2*alpha*rho0^2-d*delta+d))], d*i1-g*i2+a)
> 

G1 := subs(DATA, G08); G2 := subs(DATA, G09)

0.154260718e-1

  

0.1573686984e-1

 (4)
> 

-G1*got*we+Actual*value
> 

``
> 

NULL

Error, illegal use of an object as a name

"a+d*i1-g*i2:=0.005+ 0.02*0.2106517897-0.03*0.05;"

  
> 

``
> 

``


 

Download q_subs.mw
Please Wait...