{{{#Method 1 var('x') f(x) = x**2 - 3*x + 1 solve([f(x) == 0], x)

#Method 2 from scipy.optimize import fsolve import numpy as np

def f(x):

fsolve(f, np.random.rand(2))}}}