Fix bug #1848: taytorat leaks internal gensyms from multivar expansions

[maxima.git] / doc / info / ja / to_poly_solve.texi

@menu
* Functions and Variables for to_poly_solve::
@end menu

@node Functions and Variables for to_poly_solve,  , to_poly_solve, to_poly_solve
@section Functions and Variables for to_poly_solve

パッケージ @code{to_poly}と @code{to_poly_solve}は実験的なものです;
これらのパッケージの関数の仕様は変更の可能性があり、
これらのパッケージの関数のいくつかは他のMaxima関数にマージされるかもしれません。

Barton Willis (University of Nebraska at Kearneyの数学科の教授)が
@code{to_poly}と @code{to_poly_solve}パッケージと
これらのパッケージのための英語のユーザードキュメンテーションを書きました。

@deffn {Operator} %and
@ifinfo
@fnindex Logical conjunction
@end ifinfo

演算子 @code{%and}は整理化(simplifying)非短絡論理積です。
Maximaは
@code{%and}式をtrueかfalseか論理的に同値ですが整理された式に整理します。
演算子 @code{%and}は結合的で可換でべき等的です。
そして、
@code{%and}が名詞形を返す時、
@code{%and}の引数は非冗長なソートされたリストを形成します;
例えば、

@example
(%i1) a %and (a %and b);
(%o1)                       a %and b
@end example

もし積の引数の1つが @i{明示的に}もう一つの引数の否定なら、
@code{%and}はfalseを返します:

@example
(%i2) a %and (not a);
(%o2)                         false
@end example

積のいずれかのメンバーがfalseなら、
たとえ他のメンバーが明白に非ブーリアンでも積はfalseに整理されます;
例えば、

@example
(%i3) 42 %and false;
(%o3)                         false
@end example

inequation(すなわち、不等もしくは等式)の@code{%and}式の任意の引数は
Fourier消去パッケージを使って整理されます。
Fourier消去整理器は
すべてでないですがいくつかの非線形inequationを線形inequationに変換する
プリプロセッサを持ちます;
例えば、Fourier消去コードは
@code{abs(x) + 1 > 0}を trueに整理します。だから、

@example
(%i4) (x < 1) %and (abs(x) + 1 > 0);
(%o4)                         x < 1
@end example

@b{注釈}  
@itemize @bullet
@item
オプション変数 @code{prederror}は
@code{%and}式の整理を変え@i{ません}。

@item
演算子の優先順位エラーを避けるため、
演算子 @code{%and, %or}と @code{not}を含む複合式は完全に括弧でくくらなくてはいけません。

@item
Maxima演算子 @code{and}と @code{or}はどちらも短絡的です。
なので @code{and}は結合的でなかったり可換ではなかったりします。

@end itemize

@b{制限} 積 @code{%and}はinequationを
@i{大域的にではなく局所的に}整理します。
これは以下のような積はfalseに整理@i{されない}ことを意味します。

@example
(%i5) (x < 1) %and (x > 1);
(%o5)                 (x > 1) %and (x < 1)
@end example

また、 Fourier消去コードは事実データベースを@i{無視します}。

@example
(%i6) assume(x > 5);
(%o6)                        [x > 5]
(%i7) (x > 1) %and (x > 2);
(%o7)                 (x > 1) %and (x > 2)
@end example

最終的に、同値の線形inequationに容易に変換されない非線形inequationは整理されません。

@code{%or}上に @code{%and}を分配するためのサポートはありません;
@code{%and}上に論理否定を分配するサポートもありません。

@b{利用のためには} @file{load(to_poly_solve)}

@b{関連関数} @code{%or, %if, and, or, not}

@b{状況} 演算子 @code{%and}は実験的です;
この関数の仕様は変更の可能性があり、
その機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {演算子} %if (@var{bool}, @var{a}, @var{b})
@ifinfo
@fnindex conditional evaluation
@end ifinfo

演算子 @code{%if}は整理化条件文です。
@i{条件} @var{bool}はブーリアン値でなければいけません。
条件がtrueの時、二番目の引数を返します;
条件がfalseの時、三番目の引数を返します;
その他の場合、名詞形を返します。

Maxima inequations (不等式か等式)はブーリアン値では@i{ありません};
例えば、 Maximaは @math{5 < 6}をtrueに整理 @i{しません}し、
@math{5 = 6}をfalseに整理しません;
しかしながら、 to false; however, in
the context of a conditional to an 
@code{%if}文への条件の文脈では、
Maximaは @i{自動的に}inequationの真値を決定しようとします。
例:

@example
(%i1) f : %if(x # 1, 2, 8);
(%o1)                 %if(x - 1 # 0, 2, 8)
(%i2) [subst(x = -1,f), subst(x=1,f)];
(%o2)                        [2, 8]
@end example

もし条件がinequationを含むなら、
MaximaはFourier消去パッケージを使ってそれを整理します。

@b{注釈} 

@itemize @bullet
@item
もし条件が明白に非ブーリアンなら、Maximaは名詞形を返します:
@end itemize

@example
(%i3) %if(42,1,2);
(%o3)                     %if(42, 1, 2)
@end example

@itemize @bullet
@item
Maxima演算子 @code{if}はn項演算子ですが、演算子 @code{%if}は n項演算子@i{ではないです}。
@end itemize

@b{制限} Fourier消去コードは
同値の線形inequationに容易に変換可能な非線形inequationだけを
整理します。

@b{利用には:} @file{load(to_poly_solve)}

@b{状況:} 演算子 @code{%if}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn
  
@deffn {演算子} %or
@ifinfo
@fnindex Logical disjunction
@end ifinfo

演算子 @code{%or}は整理化非短絡論理和です。
Maximaは
@code{%or}式をtrueかfalseか論理的に同値ですが整理された式に整理します。
演算子 @code{%or}は結合的で可換でべき等です。
そして
@code{%or}が名詞形を返す時、
@code{%or}の引数は非冗長なソートされたリストを形成します;
例えば、

@example
(%i1) a %or (a %or b);
(%o1)                        a %or b
@end example

もし和のメンバーの１つが@i{明示的に}他のメンバーの否定なら、
@code{%or}はtrueを返します:

@example
(%i2) a %or (not a);
(%o2)                         true
@end example

もし和の任意のメンバーがtrueなら、
たとえ和の他のメンバーが明白に非ブーリアンでも和はtrueに整理されます。
例えば

@example
(%i3) 42 %or true;
(%o3)                         true
@end example

inequation(不等式か等式)である @code{%or}式の任意の引数は
Fourier消去パッケージを使って整理されます。
Fourier消去コードは @code{abs(x) + 1 > 0}をtrueに整理します。
なので、

@example
(%i4) (x < 1) %or (abs(x) + 1 > 0);
(%o4)                         true
@end example

@b{注釈}  
@itemize @bullet
@item
オプション変数 @code{prederror}は
@code{%or}式の整理を変更@i{しません}。

@item
演算子 @code{%and, %or}, @code{not}を含む複合式は括弧でくくらなければいけません;
これらの演算子の結合力は期待と一致しないことがあります。

@item
Maxima演算子 @code{and}と @code{or}はどちらも短絡的です。
なので @code{or}は結合的でないか可換でないかです。

@end itemize

@b{制限} 和 @code{%or}はinequationを
@i{大域的にではなく、局所的に}整理します。
これは以下のような和はtrueに整理@i{されない}ことを意味します。

@c TODO: IN MAXIMA 5.24POST THIS SIMPLIFIES TO TRUE.

@example
(%i1) (x < 1) %or (x >= 1);
(%o1) (x > 1) %or (x >= 1)
@end example

更に Fourier消去コードは事実データベースを無視します:

@example
(%i2) assume(x > 5);
(%o2)                        [x > 5]
(%i3) (x > 1) %and (x > 2);
(%o3)                 (x > 1) %and (x > 2)
@end example

最終的に、同値の線形inequationに容易に変換されない非線形inequationは整理されません。

両方ともfalseである項を探すアルゴリズムは弱いものです;
また、@code{%and}上に @code{%or}を分配するためのサポートはありません;
@code{%or}上に論理否定を分配するサポートもありません。

@b{利用のためには} @file{load(to_poly_solve)}

@b{関連関数} @code{%or, %if, and, or, not}

@b{状況} 演算子 @code{%or}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} complex_number_p (@var{x})

もし引数が
 @code{a + %i * b}, @code{a}, @code{%i b}, @code{%i}のいずれかなら、
述語論理 @code{complex_number_p}はtrueを返します。
ここで @code{a}と @code{b}は有理数か(多倍長浮動小数点数を含む)浮動小数点数です;
他のすべての入力に対して、
@code{complex_number_p}はfalseを返します; 例えば

@example
(%i1) map('complex_number_p,[2/3, 2 + 1.5 * %i, %i]);
(%o1)                  [true, true, true]
(%i2) complex_number_p((2+%i)/(5-%i));
(%o2)                         false
(%i3) complex_number_p(cos(5 - 2 * %i));
(%o3)                         false
@end example

@b{関連関数} @code{isreal_p}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{complex_number_p}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} compose_functions (@var{l})

関数コール @code{compose_functions(l)}は
リスト @var{l}の中の関数の合成をラムダ形を返します。
関数は@i{右}から @i{左}へ適用されます; 例えば

@example
(%i1) compose_functions([cos, exp]);
                                        %g151
(%o1)             lambda([%g151], cos(%e     ))
(%i2) %(x);
                                  x
(%o2)                       cos(%e )
@end example

関数リストが空のとき、恒等関数を返します:

@example
(%i3) compose_functions([]);
(%o3)                lambda([%g152], %g152)
(%i4)  %(x);
(%o4)                           x
@end example

@b{注釈} 
@itemize @bullet
@item
Maximaがリストメンバーがシンボルやラムダ形でないか決める時、
(@code{compose_functions}@i{ではなく}) @code{funmake}がエラーをシグナルします:
@end itemize

@example
(%i5) compose_functions([a < b]);

funmake: first argument must be a symbol, subscripted symbol,
string, or lambda expression; found: a < b
#0: compose_functions(l=[a < b])(to_poly_solve.mac line 40)
 -- an error. To debug this try: debugmode(true);
@end example

@itemize @bullet
@item
名前の衝突を避けるために、独立変数は
関数 @code{new_variable}が決めます。

@end itemize

@example
(%i6) compose_functions([%g0]);
(%o6)              lambda([%g154], %g0(%g154))
(%i7) compose_functions([%g0]);
(%o7)              lambda([%g155], %g0(%g155))
@end example

@itemize
@item
 独立変数が違っても、
Maximaはこれらのラムダ形が意味的に等しいことを演繹できます:
@end itemize

@example
(%i8) is(equal(%o6,%o7));
(%o8)                         true
@end example

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{compose_functions}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。
@end deffn

@deffn {関数} dfloat (@var{x})

関数 @code{dfloat}は
@code{float}に似ていますが、
@code{float}が IEEE 倍精度浮動小数点数に評価するのを失敗した時、
関数 @code{dfloat}は
@code{rectform}を適用します; 例えば

@example
(%i1) float(4.5^(1 + %i));
                               %i + 1
(%o1)                       4.5
(%i2) dfloat(4.5^(1 + %i));
(%o2)        4.48998802962884 %i + .3000124893895671
@end example

@b{注釈} 

@itemize @bullet
@item
式の直交形は
数値評価には不完全で合わないかもしれません--例えば、
直交形は不必要に浮動小数点数の差(桁落ち)を含むかもしれません。

@item
識別子 @code{float}は(デフォルト値がfalseの)オプション変数であり
かつ関数名です。

@end itemize

@b{関連関数} @code{float, bfloat}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{dfloat}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} elim (@var{l}, @var{x})

関数 @code{elim}は
集合もしくはリスト @code{l}内の等式から
集合もしくはリスト@code{x}内の変数を消去します。
@code{x}のメンバーそれぞれはシンボルでなければいけません;
 must be a symbol; the members of 
@code{l}のメンバーは等式か零に等しいと仮定される式であり得ます。

関数 @code{elim}は2つのリストのリストを返します;
一番目は変数が消去された式のリストです;
二番目はピボットのリストです;
二番目のリストは
@code{elim}が変数を消去するのに使った式のリストです。

以下は線形方程式を消去する例です:

@example
(%i1) elim(set(x + y + z = 1, x - y  - z = 8, x - z = 1), 
           set(x,y));
(%o1)            [[2 z - 7], [y + 7, z - x + 1]]
@end example

@code{x}と @code{y}の消去は一つの等式 @code{2 z - 7 = 0}をもたらします;
等式 @code{y + 7 = 0}と @code{z - z + 1 = 1}がピボットとして使われました。
これらの等式から3つの変数すべての消去は、線形系を三角化します:

@example
(%i2) elim(set(x + y + z = 1, x - y  - z = 8, x - z = 1),
           set(x,y,z));
(%o2)           [[], [2 z - 7, y + 7, z - x + 1]]
@end example

もちろん、等式は線形である必要はありません:

@example
(%i3) elim(set(x^2 - 2 * y^3 = 1,  x - y = 5), [x,y]);
                     3    2
(%o3)       [[], [2 y  - y  - 10 y - 24, y - x + 5]]
@end example

ユーザーは変数が消去される順序を制御しません。
確かに、アルゴリズムは
最良のピボットと最良の消去順序を選ぼうとする発見的方法を使います。

@b{注釈} 

@itemize @bullet
@item
関連関数 @code{eliminate}と違い、
等式の数が変数の数と等しい時、関数 @code{elim}は
@code{solve}を呼び出し@i{ません}。

@item
関数 @code{elim}は終結式を適用することで機能します;
オプション変数 @code{resultant}は
Maximaがどのアルゴリズムを使うか決めます。
@code{sqfr}を使って、 Maximaは
それぞれの終結式を因数分解し、
多重零点を抑制します。

@item
@code{elim}は多項式方程式の非線形一式を三角化します;
三角化された集合の解集合は
三角化されていない集合の解集合よりも大きい@i{かもしれません}。
例えば、、
三角化された方程式は@i{見せかけの}解を持ち得ます。
@end itemize

@b{関連関数} @i{elim_allbut, eliminate_using, eliminate}

@b{オプション変数} @i{resultant}

@b{利用するには} @file{load(to_poly)}

@b{状況} 関数 @code{elim}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn
 
@deffn {関数} elim_allbut (@var{l}, @var{x})

この関数は
 @code{elim}に似ています。
等式のリスト @code{l}内の
リスト @code{x}内の変数を @i{除いて}
変数すべてを消去することを除いて

@example
(%i1) elim_allbut([x+y = 1, x - 5*y = 1],[]);
(%o1)                 [[], [y, y + x - 1]]
(%i2) elim_allbut([x+y = 1, x - 5*y = 1],[x]);
(%o2)                [[x - 1], [y + x - 1]]
@end example

@b{利用するには} @file{load(to_poly)}

@b{オプション変数} @i{resultant}

@b{関連関数} @i{elim, eliminate_using, eliminate}

@b{状況} 関数 @code{elim_allbut}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} eliminate_using (@var{l}, @var{e}, @var{x})

Using @code{e} as the pivot, eliminate the symbol @code{x} from the
list or set of equations in @code{l}.  The function @code{eliminate_using}
returns a set.

@example
(%i1) eq : [x^2 - y^2 - z^3 , x*y - z^2 - 5, x - y + z];
               3    2    2     2
(%o1)      [- z  - y  + x , - z  + x y - 5, z - y + x]
(%i2) eliminate_using(eq,first(eq),z);
        3              2      2      3    2
(%o2) @{y  + (1 - 3 x) y  + 3 x  y - x  - x , 
                        4    3  3       2  2             4
                       y  - x  y  + 13 x  y  - 75 x y + x  + 125@}
(%i3) eliminate_using(eq,second(eq),z);
        2            2       4    3  3       2  2             4
(%o3) @{y  - 3 x y + x  + 5, y  - x  y  + 13 x  y  - 75 x y + x
                                                           + 125@}
(%i4) eliminate_using(eq, third(eq),z);
        2            2       3              2      2      3    2
(%o4) @{y  - 3 x y + x  + 5, y  + (1 - 3 x) y  + 3 x  y - x  - x @}
@end example

@b{オプション変数} @i{resultant}

@b{関連関数} @i{elim, eliminate, elim_allbut}

@b{利用するには} @file{load(to_poly)}

@b{状況} 関数 @code{elimimate_using}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} fourier_elim ([@var{eq1}, @var{eq2}, @dots{}], [@var{var1}, @var{var}, @dots{}])

Fourier消去は
Gauss消去の線形inequation(等式か不等式)に関する類似物です。
関数コール @code{fourier_elim([eq1, eq2, ...], [var1, var2, ...]}は
線形inequationsのリスト @code{[eq1, eq2, ...]}上の
変数 @code{[var1, var2, ...]}に関するFourier消去です;
例えば

@example
(%i1) fourier_elim([y-x < 5, x - y < 7, 10 < y],[x,y]);
(%o1)            [y - 5 < x, x < y + 7, 10 < y]
(%i2) fourier_elim([y-x < 5, x - y < 7, 10 < y],[y,x]);
(%o2)        [max(10, x - 7) < y, y < x + 5, 5 < x]
@end example

Eliminating first with respect to
最初に @math{x}を次に @math{y}を消去することは
@math{y}に依存する、@math{x}に関する下限と上限、
そして@math{y}に関する下限と上限(それらは数です)ををもたらします。
他の順序での消去は
@math{x}依存の@math{y}の下限と上限、
そして@math{x}の数値的下限と上限を与えます。

必要な時、 @code{fourier_elim}はinequationのリストの@emph{論理和}を返します:

@example
(%i3) fourier_elim([x # 6],[x]);
(%o3)                  [x < 6] or [6 < x]
@end example

解集合が空の時、 @code{fourier_elim}は @code{emptyset}を返し、
解集合が実数すべての時、 @code{fourier_elim}は @code{universalset}を返します;
例えば

@example
(%i4) fourier_elim([x < 1, x > 1],[x]);
(%o4)                       emptyset
(%i5) fourier_elim([minf < x, x < inf],[x]);
(%o5)                     universalset
@end example

非線形inequationに対して、 @code{fourier_elim}は
(幾分) 整理されたinequationsのリストを返します:

@example
(%i6) fourier_elim([x^3 - 1 > 0],[x]);
@group
               2                             2
(%o6) [1 < x, x  + x + 1 > 0] or [x < 1, - (x  + x + 1) > 0]
@end group
(%i7) fourier_elim([cos(x) < 1/2],[x]);
(%o7)                  [1 - 2 cos(x) > 0]
@end example

inequationのリストの代わりに、
@code{fourier_elim}の一番目の引数は
論理和か論理積であるかもしれません:

@example
(%i8) fourier_elim((x + y < 5) and (x - y >8),[x,y]);
                                              3
(%o8)            [y + 8 < x, x < 5 - y, y < - -]
                                              2
(%i9) fourier_elim(((x + y < 5) and x < 1) or  (x - y >8),[x,y]);
(%o9)          [y + 8 < x] or [x < min(1, 5 - y)]
@end example

関数 @code{fourier_elim}は
inequation演算子 @code{<, <=, >, >=, #}, @code{=}をサポートします。

Fourier消去コードは
絶対値や最小、最大関数を含むいくつかの非線形inequationを
線形inequationに変換する
プリプロセッサを持ちます．
加えて、プリプロセッサは
線形項の積か商であるいくつかの式を扱います:

@example
(%i10) fourier_elim([max(x,y) > 6, x # 8, abs(y-1) > 12],[x,y]);
(%o10) [6 < x, x < 8, y < - 11] or [8 < x, y < - 11]
 or [x < 8, 13 < y] or [x = y, 13 < y] or [8 < x, x < y, 13 < y]
 or [y < x, 13 < y]
(%i11) fourier_elim([(x+6)/(x-9) <= 6],[x]);
(%o11)           [x = 12] or [12 < x] or [x < 9]
(%i12) fourier_elim([x^2 - 1 # 0],[x]);
(%o12)      [- 1 < x, x < 1] or [1 < x] or [x < - 1]
@end example

@b{利用するには} @file{load(fourier_elim)}

@end deffn

@deffn {関数} isreal_p (@var{e})

述語論理 @code{isreal_p}は、
Maximaが
@code{e}が実数線@i{全体}上で実数値であることを
決められる時
trueを返します;
Maximaが
実数線の空でないある部分集合上で@code{e}が実数値で @i{ない}ことを
決められる時
falseを返します;
他の場合、名詞形を返します。

@example
(%i1) map('isreal_p, [-1, 0, %i, %pi]);
(%o1)               [true, true, false, true]
@end example

Maxima変数は実と仮定されます; たとえば

@example
(%i2) isreal_p(x);
(%o2)                         true
@end example

関数 @code{isreal_p}は事実データベースを検査します:

@example
(%i3) declare(z,complex)$

(%i4) isreal_p(z);
(%o4)                      isreal_p(z)
@end example

@b{制限}
頻繁すぎるほど @code{isreal_p}は
falseを返すことができる時に名詞形を返します;
対数函数は実数線全体で実数値ではありません。
だから @code{isreal_p(log(x))}はfalseを返すべきです;
しかしながら

@example
(%i5) isreal_p(log(x));
(%o5)                   isreal_p(log(x))
@end example

@b{利用するには} @file{load(to_poly_solve)}

@b{関連関数} @i{complex_number_p}

@b{状況} 関数 @code{real_p}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。
@end deffn


@deffn {関数} new_variable (type)

 @code{%[z,n,r,c,g]k}という形のユニークなシンボルを返します。
ここで @code{k}は整数です
@math{type}に許される値は @i{integer, natural_number, real, natural_number,} @i{general}です。
(自然数(natural number)は @i{非負整数}を意味します;
なので零は自然数です
自然数の定義のすべてではありませんがいくつかは零を@i{除外します}。)

@math{type}が許される値の１つでない時、
@math{type}は @math{general}にデフォルト設定されます。
整数、自然数、複素数に対して、
Maximaは自動的にこの情報を事実データベースに追加します。

@example
(%i1) map('new_variable,
          ['integer, 'natural_number, 'real, 'complex, 'general]);
(%o1)          [%z144, %n145, %r146, %c147, %g148]
(%i2) nicedummies(%);
(%o2)               [%z0, %n0, %r0, %c0, %g0]
(%i3) featurep(%z0, 'integer);
(%o3)                         true
(%i4) featurep(%n0, 'integer);
(%o4)                         true
(%i5) is(%n0 >= 0);
(%o5)                         true
(%i6) featurep(%c0, 'complex);
(%o6)                         true
@end example

@b{注釈} 一般に、@code{new_variable}の引数はクォートすべきです。
クォートは以下に似たエラーを防ぎます。

@example
(%i7) integer : 12$

(%i8) new_variable(integer);
(%o8)                         %g149
(%i9) new_variable('integer);
(%o9)                         %z150
@end example

@b{関連関数} @i{nicedummies}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{new_variable}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} nicedummies

関数 @code{nicedummies}は
@code{new_variable}が導入した式の中の変数を零から再インデックスします;

@example
(%i1) new_variable('integer) + 52 * new_variable('integer);
(%o1)                   52 %z136 + %z135
(%i2) new_variable('integer) - new_variable('integer);
(%o2)                     %z137 - %z138
(%i3) nicedummies(%);
(%o3)                       %z0 - %z1
@end example

@b{関連関数} @i{new_variable}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{nicedummies}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} parg (@var{x})

関数 @code{parg}は
複素偏角関数 @code{carg}の整理化バージョンです;
例えば

@example
(%i1) map('parg,[1,1+%i,%i, -1 + %i, -1]);
                        %pi  %pi  3 %pi
(%o1)               [0, ---, ---, -----, %pi]
                         4    2     4
@end example

一般に非定数入力に対して @code{parg}は名詞形を返します;例えば

@example
(%i2) parg(x + %i * sqrt(x));
(%o2)                 parg(x + %i sqrt(x))
@end example

@code{sign}を使って入力が正の実数か負の実数か決めることができる時
@code{parg}は非定数入力に対して非名詞形を返します。
以下は2つの例です:

@c TODO: THE FIRST RESULT IS A NOUNFORM IN MAXIMA 5.24POST

@example
(%i3) parg(abs(x));
(%o3) 0
(%i4) parg(-x^2-1);
(%o4)                          %pi
@end example

@b{注釈} @code{sign}関数は
複素数と宣言された変数(@code{declare(x,complex)})を
ほとんど無視します;
複素数と宣言された変数に対して、
@code{parg}は間違った値を返す可能性があります;
例えば

@c TODO: IN MAXIMA 5.24POST THE RESULT IS A NOUNFORM.

@example
(%i1) declare(x,complex)$

(%i2) parg(x^2 + 1);
(%o2) 0
@end example

@b{Related function} @i{carg, isreal_p}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{parg}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} real_imagpart_to_conjugate (@var{e})

関数 @code{real_imagpart_to_conjugate}は
@code{realpart}と @code{imagpart}をすべて
@code{conjugate}を含む代数的に同値な式に置き換えます。

@example
(%i1) declare(x, complex)$

(%i2) real_imagpart_to_conjugate(realpart(x) +  imagpart(x) = 3);
          conjugate(x) + x   %i (x - conjugate(x))
(%o2)     ---------------- - --------------------- = 3
                 2                     2
@end example

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{real_imagpart_to_conjugate}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} rectform_log_if_constant (@var{e})

関数 @code{rectform_if_constant}は
@code{ log(c)}の形の項をすべて
@code{rectform(log(c))}に変換します。
ここで @code{c}は宣言された定数式か明示的に宣言された定数です。

@example
(%i1) rectform_log_if_constant(log(1-%i) - log(x - %i));
                                 log(2)   %i %pi
(%o1)            - log(x - %i) + ------ - ------
                                   2        4
(%i2) declare(a,constant, b,constant)$

(%i3) rectform_log_if_constant(log(a + %i*b));
                       2    2
                  log(b  + a )
(%o3)             ------------ + %i atan2(b, a)
                       2
@end example

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{rectform_log_if_constant}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} simp_inequality (@var{e})

関数 @code{simp_inequality}は
inequationの論理積と論理和にいくつかの整理を適用します。

@b{制限} 関数 @code{simp_inequality}は少なくとも2つの点で制限されます;
第一に整理が局所的です; 例えば

@c TODO: IN MAXIMA 5.24POST THE RESULT IS SIMPLIFIED.

@example
(%i1) simp_inequality((x > minf) %and (x < 0));
(%o1) (x>1) %and (x<1)
@end example

第二に @code{simp_inequality}は事実データベースを参照しません:

@example
(%i2) assume(x > 0)$

(%i3) simp_inequality(x > 0);
(%o3)                         x > 0
@end example

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{simp_inequality}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} standardize_inverse_trig (@var{e})

この関数は
恒等式 @code{cot(x) = atan(1/x), acsc(x) = asin(1/x),}と
 @code{asec, acoth, acsch} @code{asech}に関する類似
の恒等式を式に適応します。
Abramowitz and Stegun, 
Eqs. 4.4.6 から 4.4.8までと 4.6.4から 4.6.6までを参照してください。

@b{利用するには} @file{load(to_poly_solve)}

@b{状況} 関数 @code{standardize_inverse_trig}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。
@end deffn

@deffn {関数} subst_parallel (@var{l}, @var{e})

@code{l}が単一の等式か等式のリストの時、
それぞれの等式の右辺側を左辺に代入します。
代入は並列に行われます; 例えば

@example
(%i1) load(to_poly_solve)$

(%i2) subst_parallel([x=y,y=x], [x,y]);
(%o2)                        [y, x]
@end example

これを順にされる代入と比較します:

@example
(%i3) subst([x=y,y=x],[x,y]);
(%o3)                        [x, x]
@end example

関数 @code{subst_parallel}は
@code{subst_parallel}は非アトムの代入を許すことを除いて
@code{sublis}と似ています;
例えば

@example
(%i4) subst_parallel([x^2 = a, y = b], x^2 * y);
(%o4)                          a b
(%i5) sublis([x^2 = a, y = b], x^2 * y);

                                                             2
sublis: left-hand side of equation must be a symbol; found: x
 -- an error. To debug this try: debugmode(true);
@end example

@code{subst_parallel}が行う代入は文字通りであり意味論的ではありません;
なので
@code{subst_parallel}は
@math{x * y}が @math{x^2 * y}の部分式であることを認識@i{しません}

@example
(%i6) subst_parallel([x * y = a], x^2 * y);
                               2
(%o6)                         x  y
@end example

関数 @code{subst_parallel}は整理の@i{前に}すべての代入を完了します。
これは、もし整理が前に行われていたらエラーが起こるような
条件式への代入を許します:

@example
(%i7) subst_parallel([x = 0], %if(x < 1, 5, log(x)));
(%o7)                           5
(%i8) subst([x = 0], %if(x < 1, 5, log(x)));

log: encountered log(0).
 -- an error. To debug this try: debugmode(true);
@end example

@b{関連関数} @i{subst, sublis, ratsubst}

@b{利用するには} @file{load(to_poly_solve_extra.lisp)}

@b{状況} 関数 @code{subst_parallel}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} to_poly (@var{e}, @var{l})

関数 @code{to_poly}は
等式 @code{e}を
inequation制約を伴う多項式系に変換しようとします;
制約を満たす多項式系の解は
等式 @code{e}の解です。
非公式には @code{to_poly}は
等式 @var{e}を多項式化しようとします; 
例がわかりやすいかもしれません:

@example
(%i1) load(to_poly_solve)$

(%i2) to_poly(sqrt(x) = 3, [x]);
                            2
(%o2) [[%g130 - 3, x = %g130 ], 
                      %pi                               %pi
                   [- --- < parg(%g130), parg(%g130) <= ---], []]
                       2                                 2
@end example

条件 @code{-%pi/2<parg(%g6),parg(%g6)<=%pi/2}は
@code{%g6}が平方根関数の範囲内であることを告げます。
これがtrueの時、
@code{sqrt(x) = 3}に設定された解は
@code{%g6-3,x=%g6^2}に設定された解と同じです。

三角式を多項式化するには、
非代数的代入の導入が必要です;
これらの非代数的代入は
@code{to_poly}が返す三番目のリスト内で返されます;
例えば

@example
(%i3) to_poly(cos(x),[x]);
                2                                 %i x
(%o3)    [[%g131  + 1], [2 %g131 # 0], [%g131 = %e    ]]
@end example

数字1が変数リストのメンバーでない限り定数項は多項式化されません;
例えば

@example
(%i4) to_poly(x = sqrt(5),[x]);
(%o4)                [[x - sqrt(5)], [], []]
(%i5) to_poly(x = sqrt(5),[1,x]);
                            2
(%o5) [[x - %g132, 5 = %g132 ], 
                      %pi                               %pi
                   [- --- < parg(%g132), parg(%g132) <= ---], []]
                       2                                 2
@end example

根の1つとして@math{sqrt(5) + sqrt(7)}を持つ多項式を生成するには、
以下のコマンドを使います。

@example
(%i6) first(elim_allbut(first(to_poly(x = sqrt(5) + sqrt(7),
                                      [1,x])), [x]));
                          4       2
(%o6)                   [x  - 24 x  + 4]
@end example

@b{関連関数} @i{to_poly_solve}

@b{利用するには} @file{load(to_poly)}

@b{状況:} 関数 @code{to_poly}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。

@end deffn

@deffn {関数} to_poly_solve (@var{e}, @var{l}, [options])

関数 @code{to_poly_solve}は
等式 @math{e}を
変数 @math{l}について解こうとします。
等式 @math{e}は単一の式か式の集合もしくはリストでありえます;
同様に @math{l}は単一のシンボルかシンボルの集合のリストでありえます;
@math{e}のメンバーが明示的に等式でない時(例えば @math{x^2 -1})、
ソルバは式は0に等しいと仮定します。

@code{to_poly_solve}の基本戦略は
入力を多項式形に変換するために @code{to_poly}を使い、
多項式系上で @code{algsys}をコールします。
なので Thus user options that affect 
@code{algsys}に影響を与えるユーザーオプション、特に
@code{algexact}は @code{to_poly_solve}にも影響を与えます。
@code{algexact}のデフォルト値はfalseですが、
@code{to_poly_solve}のためには、、一般的に @code{algexact}はtrueであるべきです。
関数 @code{to_poly_solve}は
@code{algexact}を局所的にtrueに設定しません。
@code{algsys}が厳密は解を決定することができない時
trueに設定すると近似解を見つけることが不可能になるためです。

@code{to_poly_solve}が解集合を決定できる時、
解集合のそれぞれのメンバーは
@code{%union}オブジェクト内のリストです:

@example
(%i1) load(to_poly_solve)$

(%i2) to_poly_solve(x*(x-1) = 0, x);
(%o2)               %union([x = 0], [x = 1])
@end example

@code{to_poly_solve}が解集合を決定@i{できない}時
@code{%solve}名詞形が返されます。(この場合、警告が印字されます。)

@example
(%i3) to_poly_solve(x^k + 2* x + 1 = 0, x);

Nonalgebraic argument given to 'to_poly'
unable to solve
                          k
(%o3)            %solve([x  + 2 x + 1 = 0], [x])
@end example

@code{%solve}名詞形への代入は時々解に帰着します。

@example
(%i4) subst(k = 2, %);
(%o4)                   %union([x = - 1])
@end example

特に三角等式の場合、
ソルバは時々任意整数の導入を必要とします。
これらの任意整数は
@code{%zXXX}という形を持ちます。
ここで @code{XXX}は整数です;
例えば

@example
(%i5) to_poly_solve(sin(x) = 0, x);
(%o5)   %union([x = 2 %pi %z33 + %pi], [x = 2 %pi %z35])
@end example

これらの変数を零から再インデックスするには @code{nicedummies}を使います:

@example
(%i6) nicedummies(%);
(%o6)    %union([x = 2 %pi %z0 + %pi], [x = 2 %pi %z1])
@end example

時折、ソルバは
@code{%cXXX}という形の任意複素数か
@code{%rXXX}という形の任意実数を導入します。
関数 @code{nicedummies}はこれらの識別子を零から再インデックスします。

解集合は時々
論理積、論理和、含意に関してそれぞれ
@code{%and}か, @code{%or}, @code{%if}を含む様々な論理演算子の
整理化バージョンを含みます;
例えば

@example
(%i7) sol : to_poly_solve(abs(x) = a, x);
(%o7) %union(%if(isnonnegative_p(a), [x = - a], %union()), 
                      %if(isnonnegative_p(a), [x = a], %union()))
(%i8) subst(a = 42, sol);
(%o8)             %union([x = - 42], [x = 42])
(%i9) subst(a = -42, sol);
(%o9)                       %union()
@end example

空集合は @code{%union}によって表現されます。

関数 @code{to_poly_solve}は
有理数や非有理べき、絶対値、三角関数、最小最大を含む等式の(すべてではありませんが)
いくつかを解くことができます。
Lambert W 函数を使って解くことができるいくつかの等式も解くことができます;
例:

@example
(%i1) load(to_poly_solve)$

(%i2) to_poly_solve(set(max(x,y) = 5, x+y = 2), set(x,y));
(%o2)      %union([x = - 3, y = 5], [x = 5, y = - 3])
(%i3) to_poly_solve(abs(1-abs(1-x)) = 10,x);
(%o3)             %union([x = - 10], [x = 12])
(%i4) to_poly_solve(set(sqrt(x) + sqrt(y) = 5, x + y = 10),
                    set(x,y));
                     3/2               3/2
                    5    %i - 10      5    %i + 10
(%o4) %union([x = - ------------, y = ------------], 
                         2                 2
                                3/2                 3/2
                               5    %i + 10        5    %i - 10
                          [x = ------------, y = - ------------])
                                    2                   2
(%i5) to_poly_solve(cos(x) * sin(x) = 1/2,x,
                    'simpfuncs = ['expand, 'nicedummies]);
                                         %pi
(%o5)              %union([x = %pi %z0 + ---])
                                          4
(%i6) to_poly_solve(x^(2*a) + x^a + 1,x);
                                        2 %i %pi %z81
                                        -------------
                                  1/a         a
                  (sqrt(3) %i - 1)    %e
(%o6) %union([x = -----------------------------------], 
                                  1/a
                                 2
@group
                                                  2 %i %pi %z83
                                                  -------------
                                            1/a         a
                          (- sqrt(3) %i - 1)    %e
                     [x = -------------------------------------])
                                           1/a
                                          2
@end group
(%i7) to_poly_solve(x * exp(x) = a, x);
(%o7)              %union([x = lambert_w(a)])
@end example

@i{線形}不等式に対して、 @code{to_poly_solve}は自動的にFourier消去を行います:

@example
(%i8) to_poly_solve([x + y < 1, x - y >= 8], [x,y]);
                               7
(%o8) %union([x = y + 8, y < - -], 
                               2
                                                              7
                                 [y + 8 < x, x < 1 - y, y < - -])
                                                              2
@end example

@code{to_poly_solve}のオプション引数のそれぞれは等式でなければいけません;
一般にこれらのオプションの順序は問題になりません．

@itemize @bullet
@item
@code{simpfuncs = l}, 
ここで @code{l}は函数のリストです。
lのメンバーの合成をそれぞれの解に適用します。
@end itemize

@example
(%i1) to_poly_solve(x^2=%i,x);
                               1/4             1/4
(%o1)       %union([x = - (- 1)   ], [x = (- 1)   ])
(%i2) to_poly_solve(x^2= %i,x, 'simpfuncs = ['rectform]);
                      %i         1             %i         1
(%o2) %union([x = - ------- - -------], [x = ------- + -------])
                    sqrt(2)   sqrt(2)        sqrt(2)   sqrt(2)
@end example

@itemize
@item
時々追加の整理は整理を戻すことがあります;
例えば
@end itemize

@example
(%i3) to_poly_solve(x^2=1,x);
(%o3)              %union([x = - 1], [x = 1])
(%i4) to_poly_solve(x^2= 1,x, 'simpfuncs = [polarform]);
                                        %i %pi
(%o4)            %union([x = 1], [x = %e      ]
@end example

@itemize
@item
Maximaは
関数リスト @code{l}のそれぞれのメンバーが純粋に整理であるかチェックしようとしません;
例えば
@end itemize

@example
(%i5) to_poly_solve(x^2 = %i,x, 'simpfuncs = [lambda([s],s^2)]);
(%o5)                   %union([x = %i])
@end example

@itemize
@item
それぞれの解を倍精度浮動小数点数に変換するには
@code{simpfunc = ['dfloat]}を使います:
@end itemize

@example
(%i6) to_poly_solve(x^3 +x + 1 = 0,x, 
                    'simpfuncs = ['dfloat]), algexact : true;
(%o6) %union([x = - .6823278038280178], 
[x = .3411639019140089 - 1.161541399997251 %i], 
[x = 1.161541399997251 %i + .3411639019140089])
@end example

@itemize @bullet
@item
@code{use_grobner = true}
このオプションを使うと
解を企てる前に関数 @code{poly_reduced_grobner}が
等式に適用されます
主として、このオプションは
関数 @code{algsys}の弱みに対する応急措置を提供します。
以下はそんな応急措置の例です:
@end itemize

@example
(%i7) to_poly_solve([x^2+y^2=2^2,(x-1)^2+(y-1)^2=2^2],[x,y],
                    'use_grobner = true);
@group
                    sqrt(7) - 1      sqrt(7) + 1
(%o7) %union([x = - -----------, y = -----------], 
                         2                2
@end group
                                 sqrt(7) + 1        sqrt(7) - 1
                            [x = -----------, y = - -----------])
                                      2                  2
(%i8) to_poly_solve([x^2+y^2=2^2,(x-1)^2+(y-1)^2=2^2],[x,y]);
(%o8)                       %union()
@end example

@itemize @bullet
@item
@code{maxdepth = k},
ここで @code{k}は正の整数です。
この関数はソルバの最大再帰深さを制御します。
@code{maxdepth}のデフォルト値は5です。
再帰深さを越えた時、ソルバはエラーをシグナルします:
@end itemize

@example
(%i9) to_poly_solve(cos(x) = x,x, 'maxdepth = 2);

Unable to solve
Unable to solve
(%o9)        %solve([cos(x) = x], [x], maxdepth = 2)
@end example

@itemize @bullet
@item
@code{parameters = l},
ここで @code{l}はシンボルのリストです。
ソルバは
リスト @code{l}のメンバーすべてに有効な解を返そうとします;
例えば:
@end itemize

@example
(%i10) to_poly_solve(a * x = x, x);
(%o10)                   %union([x = 0])
(%i11) to_poly_solve(a * x = x, x, 'parameters = [a]);
(%o11) %union(%if(a - 1 = 0, [x = %c111], %union()), 
                               %if(a - 1 # 0, [x = 0], %union()))
@end example

@itemize
@item
@code{(%o2)}で、ソルバはダミー変数を導入しました;
これらのダミー変数を再インデックスするには、
関数 @code{nicedummies}を使います:
@end itemize

@example
(%i12) nicedummies(%);
(%o12) %union(%if(a - 1 = 0, [x = %c0], %union()), 
                               %if(a - 1 # 0, [x = 0], %union()))
@end example

@math{f(a) = f(b)}という形の等式を解くために
@code{to_poly_solve}は
ハッシュ配列 @code{one_to_one_reduce}に記憶されたデータを使います。
割り当て @code{one_to_one_reduce['f,'f] : lambda([a,b], a=b)}は
@code{to_poly_solve}に
@math{f(a) = f(b)}の解集合が@math{a=b}の解集合に等しいことを告げます;
例えば

@example
(%i13) one_to_one_reduce['f,'f] : lambda([a,b], a=b)$

(%i14) to_poly_solve(f(x^2-1) = f(0),x);
(%o14)             %union([x = - 1], [x = 1])
@end example

さらに一般的に、
割り当て @code{one_to_one_reduce['f,'g] : lambda([a,b], w(a, b) = 0}は
@code{to_poly_solve}に
@math{f(a) = f(b)}の解集合が @math{w(a,b) = 0}の解集合に等しいことを告げます;
例えば

@example
(%i15) one_to_one_reduce['f,'g] : lambda([a,b], a = 1 + b/2)$

(%i16) to_poly_solve(f(x) - g(x),x);
(%o16)                   %union([x = 2])
@end example

加えて、@math{f(a) = b}という形の等式を解くために
関数 @code{to_poly_solve}は
ハッシュ配列 @code{function_inverse}に記憶されたデータを使います。
割り当て @code{function_inverse['f] : lambda([s], g(s))}は
@code{to_poly_solve}に
@code{f(x) = b}に設定された解が
@code{x = g(b)}に設定された解に等しいことを知らせます;
2つの例:

@example
(%i17) function_inverse['Q] : lambda([s], P(s))$

(%i18) to_poly_solve(Q(x-1) = 2009,x);
(%o18)              %union([x = P(2009) + 1])
(%i19) function_inverse['G] : lambda([s], s+new_variable(integer));
(%o19)       lambda([s], s + new_variable(integer))
(%i20) to_poly_solve(G(x - a) = b,x);
(%o20)             %union([x = b + a + %z125])
@end example


@b{注釈}

@itemize
@item
解く変数はシンボルである必要はありません;
@code{fullratsubst}が適切に代入できる時、
解く変数は非シンボルであり得ます:
@end itemize

@example
(%i1) to_poly_solve([x^2 + y^2 + x * y = 5, x * y = 8],
                    [x^2 + y^2, x * y]);
                                  2    2
(%o1)           %union([x y = 8, y  + x  = - 3])
@end example

@itemize
@item
複素共役を含む等式に対して、
ソルバは自動的に共役等式を追加します;
例えば
@end itemize

@example
(%i1) declare(x,complex)$

(%i2) to_poly_solve(x + (5 + %i) * conjugate(x) = 1, x);
                                   %i + 21
(%o2)              %union([x = - -----------])
                                 25 %i - 125
(%i3) declare(y,complex)$

(%i4) to_poly_solve(set(conjugate(x) - y = 42 + %i,
                        x + conjugate(y) = 0), set(x,y));
                           %i - 42        %i + 42
(%o4)        %union([x = - -------, y = - -------])
                              2              2
@end example

@itemize
@item
絶対値関数を含む等式に対して、
@code{to_poly_solve}は
絶対値の引数が複素数値か決めるため
事実データベースに問い合わせます。
@end itemize

@example
(%i1) to_poly_solve(abs(x) = 6, x);
(%o1)              %union([x = - 6], [x = 6])
(%i2) declare(z,complex)$

(%i3) to_poly_solve(abs(z) = 6, z);
(%o3) %union(%if((%c11 # 0) %and (%c11 conjugate(%c11) - 36 = 
                                       0), [z = %c11], %union()))
@end example

@itemize
@item
@i{これはソルバが事実データベースを参照する唯一の状況です。例えば、もし解く変数が整数と宣言されていたら、@code{to_poly_solve}はこの宣言を無視します}。
@end itemize

@b{Relevant option variables} @i{algexact, resultant, algebraic}

@b{関連関数} @i{to_poly}

@b{利用するには} @file{load(to_poly_solve)}

@b{状況:} 関数 @code{to_poly_solve}は実験的です;
仕様は変更の可能性があり、
機能は他のMaxima関数にマージされるかもしれません。
@end deffn