функции simplify существенно зависит от областей определения переменных. В следующем примере упрощение выражения не произошло, поскольку результат этой операции неоднозначен:
> restart;
> simplify(sqrt(х^4*у^2));
Однако, определив переменные как реальные или положительные, можно легко добиться желаемого упрощения:
> simplify(sqrt(х^4*у^2),assume=positive);
x² у
> simplify(sqrt(х^4*у^2),assume=real);
x²|y|
С помощью равенств можно задать свои правила преобразования, например:
> eq:=x^2+2*x*y+y^2;
eq:=х² +2ху + y²
> simplify(eq,{х=1));
y² + 2y + 1
> simplify(eq,{х^2=х*у, у^2=1});
3хy + 1
> simplify(eq,{х,у});
0
Обратите внимание на то, что указание в списке равенств только левой части равенства означает, что правая часть принимается равной нулю. Если функция simplify не способна выполнить упрощение выражения expr, то она просто его повторяет. Это сигнал к применению опций, уточняющих преобразования.
Сложность упрощаемых выражений зависит от объема ОЗУ и вида интерфейса. Очень большие выражения надо разбивать на подвыражения и работать с ними раздельно.
3.7.2. Расширение выражений — expand
Даже в жизни мы говорим: «не все так просто». Порою упрощенное выражение скрывает его особенности, знание которых является желательным. В этом случае можно говорить о полезности расширения или раскрытия выражения. Функция expand «расширяет» выражение expr и записывается в виде
expand(expr, expr1, expr2, ..., exprn)
где expr — расширяемое выражение, expr1, expr2, …, exprn — необязательные подвыражения — опции. Имеется также инертная форма данной функции — Ехpand(expr). Кроме того, возможно применение операторной конструкции frontend(expans,[expr]).
Функция expand
раскладывает рациональные выражения на простые дроби, полиномы на полиномиальные разложения, она способна раскрыть многие математические функции, такие как sin, cos, tan, sinh, cosh, tanh, det, erf, exp, factorial, GAMMA, ln, max, min, Psi, binomial, sum, product, int, limit, bernoulli, euler, abs, signum, pochhammer, polylog, BesselJ, BesselY, BesselI, BesselK, AngerJ, Beta, Hankel, Kelvin, Struve, WeberE и функция piecewise. С помощью дополнительных аргументов expr1, expr2, …, exprn можно задать расширение отдельных фрагментов в expr.
Примеры применения функции expand приведены ниже (файл expand):
> expand((х+2)*(х+3)*(х+4));
x³ + 9х² + 26х + 24
> expand(sin(2*х));
2sin(x)cos(x)
> expand(sin(х+у));
sin(x)cos(y) +cos(x)sin(y)
> expand([(a+b)*(a-b),tan(2*x)]);
> expand((a+d)*(b+d)*(c+d));
abc + abd + adc + ad² + dbc + d²b + d²с = d³
> expand((х+1)*(y+1));
xy + х + у + 1
> expand((у+1),(х+1));
y + 1
> expand( (х+1) *(у+z));
ху + xz + y +z
> expand((х+1)*(y+z), х+1);
(х + 1)y +(х + 1)z
> frontend(expand,[(a+b)^3]);
а³ + 3a²b + 3аb²+b³
3.7.3. Разложение целых и рациональных чисел — ifactor
Для разложения целых или рациональных чисел на множители в виде простых чисел служит функция
ifactor(n)
или
ifactor(n,method)
где n — число, method — параметр, задающий метод разложения. Другая библиотечная функция, ifactors(n), возвращает результат разложения в форме вложенных списков (файл factor):