cineqs4 使用例
作動確認バージョン.
Maxima 5.39.0 using Lisp CMU Common Lisp 21b (21B Unicode)
高次不等式.
(%i1) ep(1,0)$cineqs4('( (x-1)*(x^2-2)*(x^3-3)*(x^2+4)<=0 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 856 bytes. Evaluation took 0.0600 seconds (0.0600 elapsed) using 1.337 MB. (%o2) [[-sqrt(2) <= x,x <= 1],[sqrt(2) <= x,x <= 3^(1/3)]]
代数的な表示が出来ない場合は,近似値(内部では小数第6位以下を切り捨てています)で答えます.
(%i3) ep(1,0)$cineqs4('( (x-1)*(x^2-2)*(x^3-3)*(x^2+4)<=1/10 )); Evaluation took 0.0000 seconds (0.0100 elapsed) using 720 bytes. Evaluation took 0.0100 seconds (0.0100 elapsed) using 1.091 MB. (%o4) [[-1.414631932468004 <= x,x <= 1.010331702011963], [1.374072765807135 <= x,x <= 1.473559962228517]]
分数関数で表された不等式.
(%i5) ep(1,0)$cineqs4('( (x-1)/x<2 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0000 seconds (0.0100 elapsed) using 520.242 KB. (%o6) [[x < -1],[0 < x]]
連言.
(%i7) ep(1,0)$cineqs4('( (x-1)*(x^2-2)*(x^3-3)<=0 and (x-1)/x<2 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0400 seconds (0.0500 elapsed) using 2.571 MB. (%o8) [[-sqrt(2) <= x,x < -1],[0 < x,x <= 1],[sqrt(2) <= x,x <= 3^(1/3)]]
選言.
(%i9) ep(1,0)$cineqs4('( (x-1)*(x^2-2)*(x^3-3)<=0 or (x-1)/x<2 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0300 seconds (0.0300 elapsed) using 2.804 MB. (%o10) [[true]]
こんな場合は...
(%i11) ep(1,0)$cineqs4('( x^4-x-1<=0 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.2500 seconds (0.2500 elapsed) using 24.248 MB. (%o12) [[-(sqrt(sqrt(6)*sqrt((sqrt(849)+9)^(1/3)*(3*sqrt(849)+155)^(1/3) *(2^(8/3)*18^(2/3)*sqrt(849) -9*2^(8/3)*18^(2/3)) +(sqrt(849)+9)^(1/3)*(3*sqrt(849)+155)^(2/3) *(2^(1/3)*18^(2/3)*sqrt(849) -9*2^(1/3)*18^(2/3)) +(sqrt(849)+9)^(1/3) *(32*18^(2/3)*sqrt(849)-16*18^(5/3))) +3*2^(11/3)*(3^(3/2)*sqrt(283)-27)^(1/3) -8*18^(2/3)*(sqrt(3)*sqrt(283)+9)^(1/3)) -sqrt(6)*sqrt(2^(1/3)*(3*sqrt(849)-27)^(1/3)*(3*sqrt(849)+155)^(1/3) -2^(8/3)*(3*sqrt(849)-27)^(1/3))) /24 <= x, x <= (sqrt(sqrt(6)*sqrt((sqrt(849)+9)^(1/3)*(3*sqrt(849)+155)^(1/3) *(2^(8/3)*18^(2/3) *sqrt(849) -9*2^(8/3)*18^(2/3)) +(sqrt(849)+9)^(1/3)*(3*sqrt(849)+155)^(2/3) *(2^(1/3)*18^(2/3) *sqrt(849) -9*2^(1/3)*18^(2/3)) +(sqrt(849)+9)^(1/3) *(32*18^(2/3)*sqrt(849)-16*18^(5/3))) +3*2^(11/3)*(3^(3/2)*sqrt(283)-27)^(1/3) -8*18^(2/3)*(sqrt(3)*sqrt(283)+9)^(1/3)) +sqrt(6)*sqrt(2^(1/3)*(3*sqrt(849)-27)^(1/3) *(3*sqrt(849)+155)^(1/3) -2^(8/3)*(3*sqrt(849)-27)^(1/3))) /24]]
コントローラー ep の第一引数を 0 にすると,近似値で答えます.
(%i13) ep(0,0)$cineqs4('( x^4-x-1<=0 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 736 bytes. Evaluation took 0.0000 seconds (0.0100 elapsed) using 175.711 KB. (%o14) [[-0.7244919786096257 <= x,x <= 1.220744081172491]]
無理関数で表された不等式.
(%i15) ep(1,0)$cineqs4('( sqrt(x)<x-1 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0100 seconds (0.0100 elapsed) using 674.156 KB. (%o16) [[(sqrt(5)+3)/2 < x]]
含意.
(%i17) ep(1,0)$cineqs4('( 1<sqrt(x) implies 1<x )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0100 seconds (0.0100 elapsed) using 428.797 KB. (%o18) [[true]]
根号が複数あっても.
(%i19) ep(1,0)$cineqs4('( sqrt(x)+sqrt(3-x)<2 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0300 seconds (0.0200 elapsed) using 1.524 MB. (%o20) [[0 <= x,x < -(2^(3/2)-3)/2],[(2^(3/2)+3)/2 < x,x <= 3]] (%i21) ep(0,0)$cineqs4('( sqrt(x)+sqrt(3-x)+sqrt(3+x)<22/5 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 736 bytes. Evaluation took 0.0200 seconds (0.0200 elapsed) using 2.500 MB. (%o22) [[0 <= x,x < 0.9647335423197492],[2.937906564163217 < x,x <= 3]]
ネストされても.
(%i23) ep(1,0)$cineqs4('( 1/sqrt(sqrt(1+x)-x)<2 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0700 seconds (0.0800 elapsed) using 4.835 MB. (%o24) [[-1 <= x,x < 5/4]]
有理数冪でも.
(%i25) ep(1,0)$cineqs4('( x^4-x^(1/3)-1>0 )); Evaluation took 0.0000 seconds (0.0000 elapsed) using 720 bytes. Evaluation took 0.0300 seconds (0.0400 elapsed) using 1.039 MB. (%o26) [[x < -0.6196702671972711],[1.198342827550492 < x]]