제약조건에서 주어진 비선형 수식이 최대값을 갖도록 하는 변수 C++로 찾고싶어요

주어진 제약조건에서 아래 수식을 최대값이 나오게 하는것을 c++로 계산하는 건데 라이브러리를 써도 좋거든요

수식은 ln(1+a*x+b*y+c+z) : 자연로그

(a,b,c는 사용자가 입력하는 숫자이고
x,y,z가 도출하고자 하는 변수)

(변수 제약조건은 x,y,z가 양수이고 x+y+z<=1)

- dc official App