まず，次をご覧ください． (%i1) f:((x-1/x)^2+(y^2-4*y+5)^3+(z+1)^2/z^4<=1)%and(x<0); 2 (z + 1) 2 3 1 2 (%o1) (x < 0) %and (-------- + (y - 4 y + 5) + (x - -) <= 1) 4 x z 人間が見るとバレバレですが，ストレートな qepmax では (%i2) qe([],f); 2…

solve は非常に強力な関数で (%i1) solve(x=y,[x,y]); (%o1) [[x = %r1, y = %r1]] と解いてくれます（笑）．これを抑えるには (%i2) solve(x=y,[x,y]),linsolve_params:false; (%o2) [[x = y]] とすればよいのですが，次数が上がると，聞く耳を持たないご様…

選言標準形，連言標準形は代数で言えば，それぞれ展開，因数分解に当たる基本的なものであり，前者は所謂（排他的とは限らない）場合分けなので，方程式系の解の構成には欠かせません．その実装は簡単で，例えば，A ∧ (B ∨ C) を (A ∧ B) ∨ (A ∧ C) に可能な…