Question: Expression not correct? Solving a system of equations

So I have this system of equations with which I am not sure if the result is the same or not using "series" and "limit" or what is going on here.

I hope it is clear what I mean.


 

restart; with(MathematicalFunctions); Assume(k__2H2O > 0, `k__HA+OH` > 0, `k__A+H2O` > 0, `k__H3O+OH` > 0, `k__HA+H2O` > 0, `k__H3O+A` > 0, HA__0 > 0, H2O > 0); sys := k__2H2O*H2O^2+`k__A+H2O`*H2O*(HA__0-HA)-(H3O*`k__H3O+OH`+HA*`k__HA+OH`)*OH = 0, k__2H2O*H2O^2+`k__HA+H2O`*H2O*HA-(`k__H3O+A`*(HA__0-HA)+`k__H3O+OH`*OH)*H3O = 0, (H2O*`k__HA+H2O`+OH*`k__HA+OH`)*HA-(H2O*`k__A+H2O`+H3O*`k__H3O+A`)*(HA__0-HA) = 0; sys := `~`[simplify]([eval(eval(sys, HA = HA__0+OH-H3O), HA__0 = x__HA0*H2O)]); sol := solve(sys, [OH, H3O]); sol := sol[1]; OH__sol := simplify(rhs(sol[1])); H3O__sol := simplify(rhs(sol[2])); simplify(OH__sol*H3O__sol); OHH3O := simplify(limit(%, `k__HA+OH` = 0)); series(OHH3O, x__HA0 = 0, 2); collect(convert(%, polynom), x__HA0, simplify, factor); r1 := limit(%, x__HA0 = 0); r2 := radnormal(limit(OHH3O, x__HA0 = 0)); simplify(r1-r2)

[`&Intersect`, `&Minus`, `&Union`, Assume, Coulditbe, Evalf, Get, Is, SearchFunction, Sequences, Series]

 

{H2O::(RealRange(Open(0), infinity))}, {HA__0::(RealRange(Open(0), infinity))}, {k__2H2O::(RealRange(Open(0), infinity))}, {`k__A+H2O`::(RealRange(Open(0), infinity))}, {`k__H3O+A`::(RealRange(Open(0), infinity))}, {`k__H3O+OH`::(RealRange(Open(0), infinity))}, {`k__HA+H2O`::(RealRange(Open(0), infinity))}, {`k__HA+OH`::(RealRange(Open(0), infinity))}

 

k__2H2O*H2O^2+`k__A+H2O`*H2O*(HA__0-HA)-(H3O*`k__H3O+OH`+HA*`k__HA+OH`)*OH = 0, k__2H2O*H2O^2+`k__HA+H2O`*H2O*HA-(`k__H3O+A`*(HA__0-HA)+`k__H3O+OH`*OH)*H3O = 0, (H2O*`k__HA+H2O`+OH*`k__HA+OH`)*HA-(H2O*`k__A+H2O`+H3O*`k__H3O+A`)*(HA__0-HA) = 0

 

[-OH^2*`k__HA+OH`+((-x__HA0*`k__HA+OH`-`k__A+H2O`)*H2O+H3O*(`k__HA+OH`-`k__H3O+OH`))*OH+k__2H2O*H2O^2+`k__A+H2O`*H2O*H3O = 0, (x__HA0*`k__HA+H2O`+k__2H2O)*H2O^2+`k__HA+H2O`*(OH-H3O)*H2O+(-`k__H3O+A`*H3O+OH*(`k__H3O+A`-`k__H3O+OH`))*H3O = 0, H2O^2*x__HA0*`k__HA+H2O`+((x__HA0*`k__HA+OH`+`k__A+H2O`+`k__HA+H2O`)*OH-H3O*(`k__A+H2O`+`k__HA+H2O`))*H2O+(OH-H3O)*(H3O*`k__H3O+A`+OH*`k__HA+OH`) = 0]

 

-RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)*H2O^2*(-`k__A+H2O`*RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)+k__2H2O)/(-`k__H3O+OH`*RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)+`k__A+H2O`)

 

-(k__2H2O*`k__H3O+A`^2-2*`k__A+H2O`^2*`k__H3O+A`+`k__A+H2O`^2*`k__H3O+OH`-2*`k__A+H2O`*`k__H3O+A`*`k__HA+H2O`+2*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`+`k__H3O+OH`*`k__HA+H2O`^2)*`k__A+H2O`*`k__HA+H2O`*H2O^2*x__HA0/((`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`)*(k__2H2O*`k__H3O+A`^2-`k__A+H2O`^2*`k__H3O+OH`-2*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-`k__H3O+OH`*`k__HA+H2O`^2))-(k__2H2O*`k__H3O+A`-`k__A+H2O`^2-`k__A+H2O`*`k__HA+H2O`)*H2O^2*(`k__A+H2O`+`k__HA+H2O`)/(`k__H3O+A`*(`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`))

 

-(k__2H2O*`k__H3O+A`-`k__A+H2O`^2-`k__A+H2O`*`k__HA+H2O`)*H2O^2*(`k__A+H2O`+`k__HA+H2O`)/(`k__H3O+A`*(`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`))

 

k__2H2O*H2O^2/`k__H3O+OH`

 

-`k__A+H2O`*(-(`k__A+H2O`+`k__HA+H2O`)^2*`k__H3O+OH`+k__2H2O*`k__H3O+A`^2)*H2O^2/(`k__H3O+OH`*`k__H3O+A`*((-`k__A+H2O`-`k__HA+H2O`)*`k__H3O+OH`+`k__A+H2O`*`k__H3O+A`))

(1)

``


 

Download Mapleprimes_-_Ionproduct.mw

Please Wait...