Здравствуйте, я пытаюсь найти глобальный максимум следующей функции, но не могу понять, как это сделать. Сначала попробовал с PuLP, но я думаю, что это не сработало, потому что он не линейный, затем попробовал библиотеку sympy с надеждой на найти 1-ю производную и затем равную 0, но безуспешно.
Может быть, вы, ребята, знаете лучший способ сделать это, дайте мне знать
y=(np.exp(intercept_B+coef_B*dsc1*x))*((dsc1)*x*NF-CR)\
+ (np.exp(intercept_B+coef_B*dsc2*x))*((dsc2)*x*NF-CR)\
+ (np.exp(intercept_B+coef_B*dsc3*x))*((dsc3)*x*NF-CR)\
+ (np.exp(intercept_B+coef_B*dsc4*x))*((dsc4)*x*NF-CR)\
+ (np.exp(intercept_B+coef_B*dsc5*x))*((dsc5)*x*NF-CR)\
+ (np.exp(intercept_B+coef_B*dsc6*x))*((dsc6)*x*NF-CR)`
единственная переменная здесь - это х и остальные - это константы, которые я получаю из входных данных заранее, моя цель - найти максимальную точку y.