optimize非线性求解

"""
计算非线性方程组：
5x1+3 = 0
4x0^2-2sin(x1x2)=0
x1x2-1.5=0
"""

import numpy as np
import math
from scipy import optimize
def fun(x):
    x0, x1, x2 = x.tolist()
    return[5 * x1 + 3, 4 * x0 * x0 - 2 * math.sin(x1 * x2), x1 * x2 - 1.5]

result = optimize.fsolve(fun,[1,1,1])
print(result)