Python Forum
Sympy nonlinsolve does not return explicit solution
Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Sympy nonlinsolve does not return explicit solution
#2
Can you explain the problem, what the variables symbolize and the equations in the equation system? A print of the solution is very confusing, and it is obvious that the unknowns you hope to eliminate are not eliminated.
SbSol  = (-E*Re**2*Rl + E*Re*Rl**2 + 2*Re**2*c1 + Re*Rl**2*c1 - 2*Re*Rl*c1 - Re*Rl*c2b*teta + Re*Rl*c2b + 2*Re*c1*teta - 2*Re*c1 + Rl**2*c1*teta - Rl**2*c1 + Rl**2*c2b*teta - Rl**2*c2b - 2*Rl*c1*teta + 2*Rl*c1)/(Re*(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))
cSol   = Complement(FiniteSet(Re*c1*(Re + teta - 1)/(E*Re**2 - E*Re*Rl - Re*Rl*c1 + 2*Re*c1 + Re*c2b*teta - Re*c2b - Rl*c1*teta + Rl*c1 - Rl*c2b*teta + Rl*c2b + 2*c1*teta - 2*c1)), FiniteSet(0))
DSol   = (E*Re**2 - E*Re*Rl - Re*Rl*c1 + 2*Re*c1 + Re*c2b*teta - Re*c2b - Rl*c1*teta + Rl*c1 - Rl*c2b*teta + Rl*c2b + 2*c1*teta - 2*c1)/(Re*(Re + teta - 1))
KSol   = (E*Re*Rl + E*Re*teta - E*Re + Re*Rl*c1 - 2*Re*c1 - Re*c2b*teta + Re*c2b + Rl*c1*teta - Rl*c1 + Rl*c2b*teta - Rl*c2b - 2*c1*teta + 2*c1)/(Re*(Re + teta - 1))
LbSol  = (E*Re**2 - E*Re*Rl - Re*Rl*c1 + Re*c2b*teta - Re*c2b - Rl*c1*teta + Rl*c1 - Rl*c2b*teta + Rl*c2b)/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl)
LfSol  = Rl*c1/(Re - Rl)
SfSol  = (E*Re**2*(Re + teta - 1) - E*Re*Rl*(Re + teta - 1) - Re*Rl*c1*(Re + teta - 1) - (Re - Rl)*(E*Re**2 - E*Re*Rl - Re*Rl*c1 + 2*Re*c1 + Re*c2b*teta - Re*c2b - Rl*c1*teta + Rl*c1 - Rl*c2b*teta + Rl*c2b + 2*c1*teta - 2*c1))/(Re*(Re - Rl)*(Re + teta - 1))
c1Sol  = c1
c2Sol  = exp((teta*log(c1) + teta*log(c2b) - teta*log((E*Re + c2b*teta - c2b)/(Re + teta - 1)) - log(c2b))/(teta - 1))
c1bSol = (E*Re + c2b*teta - c2b)/(Re + teta - 1)
c2bSol = Complement(FiniteSet(c2b), ConditionSet(c2b, Eq(c2b**(2*Re**2*Rl**3*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(2*Re**3*Rl**4*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(2*Re**4*Rl**2*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(2*Re**2*Rl**3*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(4*Re**2*Rl**4*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(4*Re**3*Rl**2*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(2*Re**3*Rl**3*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*c2b**(2*Re**4*Rl**3*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(2*Re**2*Rl**3*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(2*Re**2*Rl**4*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(2*Re**3*Rl**2*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(2*Re**3*Rl**4*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(2*Re**4*Rl**2*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**2*Rl**4*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**3*Rl**2*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**4*Rl**3*c1**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(4*Re**2*Rl**3*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**2*Rl**4*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**3*Rl**2*c1**2*teta**2/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**3*Rl**3*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**3*Rl**4*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl))*(c1**(teta/(teta - 1))*c2b**(teta/(teta - 1))*c2b**(-1/(teta - 1))*((E*Re + c2b*teta - c2b)/(Re + teta - 1))**(-teta/(teta - 1)))**(2*Re**4*Rl**2*c1**2*teta/(Re**2 - Re*Rl + Re*teta - Re - Rl*teta + Rl)), 0), Complexes))
Reply


Messages In This Thread
RE: Sympy nonlinsolve does not return explicit solution - by Serafim - Jan-20-2021, 11:38 AM

Possibly Related Threads…
Thread Author Replies Views Last Post
  Python sympy problem for symbolic equation definition joellapointe_engineering 0 287 Mar-24-2024, 11:09 PM
Last Post: joellapointe_engineering
  plot differential equation sympy dan_adi 0 3,035 Oct-13-2020, 10:44 AM
Last Post: dan_adi
  Attribute error - Sympy VictorG8 1 5,049 Apr-09-2020, 06:03 PM
Last Post: snippsat
  How to use list of symbols for sympy calculation in python? tanmaybhise 1 2,834 Mar-01-2020, 10:36 PM
Last Post: Gribouillis
  Next zero of a function on sympy VicenteM 0 2,109 Aug-21-2019, 04:58 PM
Last Post: VicenteM
  sympy psosmol 9 5,131 Mar-11-2019, 11:09 PM
Last Post: scidam
  How to return values from For Loop via return in a function AykutRobotics 14 8,422 Jan-08-2019, 04:12 AM
Last Post: stullis

Forum Jump:

User Panel Messages

Announcements
Announcement #1 8/1/2020
Announcement #2 8/2/2020
Announcement #3 8/6/2020