В последнее время приходится решать такую задачу. Дана функция многих переменных, например

, которая содержит параметр

. Необходимо найти максимум этой функции при натуральных значениях переменных, или как можно точнее оценить максимальное ее значение сверху, например

при некоторых условиях типа равенство и неравенство на переменные, ограничениях области их значений (зависящих от

). При отсутствии такого параметра часто легко воспользоваться функцией Лагранжа. А вот с ним получается совсем трудно. Хотелось бы узнать, есть ли автоматизированные средства, позволяющие решать подобные задачи? Можно ли как-то еще решать?