Как сделать линейную апроксимацию функции на python?

from matplotlib import pyplot as plt
import numpy as np
#функция
x = np.arange(0.1, 10, 1*(10**-2))
y = np.sin(x) + np.log(x)
#полином 1 степени по функции
p = np.polyfit(x,y, 1)
#подставляем значения x к полученному полиному
ya = np.polyval(p, x)

plt.plot(x, y)
plt.plot(x, ya)

plt.show()