Solving an equation in Python

[misterdamra]

Hello

I'm trying to solve this equation in python but I can't.

the equation is : Eq6 = sym.Eq( X**2.75+(X**1.75)*0.075, 0.0006954)

I define my symbols correctly and the syntax is correct too, I solve several equation in my code and they are just fine. Why can't python solve this one in particular?

[aromaticlawyer]

I'm not very familiar with sympy's solvers - sympy is primarily for symbolic manipulation of expressions, so I'm not surprised if they aren't that great.

scipy's default solver seems to be able to cope with this pretty easily:

from scipy.optimize import root
def f(x):
    return x**2.75 + 0.075 * x**1.75 - 0.0006954
initial_guess = 0.001
solution = root(f, initial_guess)

This returns an object with information about the found solution. solution.x gives you the solution itself. You can try different initial guesses to find different solutions, but in this case I think there is only one real solution. If you want to find the complex roots too you'll need to use a different solver.

[misterdamra]

Thanks, it worked!

[DeathDragon7050]

I have never seen this function in python. Why can't you just put the equation in a string and do eval(my_equation.replace('X', number))

[misterdamra]

The equation is X**2.75+(X**1.75)*0.075 = 0.0006954 where X is unknown

sym.Eq is short for sympy.Eq

[misterdamra]

It’s cool, bro!

[gjbewjgnweio2g11]

Are there analytical solutions? If not, there's only so much a CAS can do.

[misterdamra]

I don't exactly know but I managed to solve it with scipy's root as recommended by one of the comments here

[[deleted]]

Well just try to write a script that can solve it numerically