python计算两个非线性函数的交点坐标

以y1,y2两个函数为例

y1 = 0.6367212882840221 + (-0.009728567579803794 * x) + (-0.0001003343355083918      * x * x) + (1.465336233127298e-006 * x * x * x)

y2 = 0.1819125137064415 + 0.8106280309507739 / x

float数据无法直接进行==计算，所以求y1，y2做差后的绝对值小于一个精确度即可

import numpy as np

x = 1.0
percision = 0.001
x1 = np.linspace(1, 40, 1000)
for x in x1:
    # print(x)
    y1 = 0.6367212882840221 + (-0.009728567579803794 * x) + (-0.0001003343355083918 * x * x) + (
            1.465336233127298e-006 * x * x * x)
    y2 = 0.1819125137064415 + 0.8106280309507739 / x
    if abs(y1 - y2) <= percision:
        print(x, y1)