最近在搞一个项目,利用到了最优化算法,在Matlab里边有相应的函数fmincon实现,现在需要移植到C++当中,但是一直无法实现最优化的这个算法,求论坛的各位大神麻烦帮忙解决啊。要求解的方程为

其中ai0、ai1、ai2为与i有关系的常数,求在y取得最大值时x1与x2的取值,x1与x2的取值范围均为-50到50,不能用for循环反复迭代,因为步长不好确定,精度不够,步长太小了,循环次数太多,所以希望大家能够帮忙给出利用最优化求解的C++方法啊,或者是开源库也可以啊,如果是开源库的话,希望能跟大致说说如何实现啊,谢谢了,各位大神帮忙啊!