Здравствуйте, есть функция и нужно найти ее экстремумы. Пытался решить эту задачу через Maple, после выражения x у меня появляются комплексные числа, можете пожалуйста объяснить что я делаю не так или подсказать альтернативные способы?
Код:
pa1 := x -> liq1*Q96*pb1/(liq1*Q96 + x*pb1);
pa2 := x -> pb2 + x*Q96/liq2;
fn2 := x -> liq1*(pb1 - pa1(x))/Q96;
fn1 := x -> liq2*Q96*(pb2 - pa2(x))/(pa2(x)*pb2);
сама функция, экстремумы которой нужно найти:
Код:
fn := fn2(fn1(x)) - x
liq1, liq2, pb1, pb2 - все это натуральные числа.
q96 - 2^96
Вот мои действия:
Код:
dfn := diff(fn, x) // получил дифференциал функции
Получаю вот такое уравнение:
Код:
dfn := liq1^2*pb1*(-Q96^2*pb1/((pb2 + x*Q96/liq2)*pb2) + Q96^3*x*pb1/((pb2 + x*Q96/liq2)^2*pb2*liq2))/(liq1*Q96 - Q96^2*x*pb1/((pb2 + x*Q96/liq2)*pb2))^2 - 1
Код:
solve(dfn, {x}) // выражаю x
Получил такие формулы:
Код:
{x = (-pb2 + pb1*I)*pb2*liq2*liq1/((liq1*pb2 - liq2*pb1)*Q96)},
{x = (-pb2 - pb1*I)*pb2*liq2*liq1/((liq1*pb2 - liq2*pb1)*Q96)}
Мне нужно получить уравнение, результат которого будет вещественное число, а почему возникло вдруг комплексное - не понимаю :)