Artifact 2f1bc82a9ad3ea21f052ad87f3b93042bafb42a287de413209667b9ba6d737bf:
- Executable file
r38/packages/solve/ineq.rlg
— part of check-in
[f2fda60abd]
at
2011-09-02 18:13:33
on branch master
— Some historical releases purely for archival purposes
git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 1758) [annotate] [blame] [check-ins using] [more...]
Tue Feb 10 12:26:27 2004 run on Linux % polynomial Inequality (Example where another system returned {1 <= x}) ineq_solve( (2*x^2+x-1)/(x-1) >= (x+1/2)^2 ,x); {x=( - 0.894358 .. 0.326583),x=(1 .. 2.56777)} ineq_solve({(2*x^2+x-1)/(x-1) >= (x+1/2)^2, x>0}); {x=(0 .. 0.326583),x=(1 .. 2.56777)} ineq_solve({(2*x^2+x-1)/(x-1) >= (x+1/2)^2, x<-1}); {} % Systems for determining indices of Jacobi polynomials (Winfried Neun). reg := {2*a - 3>=0, 3>=0, 3>=0, 1>=0, 1>=0, 5>=0, 4>=0, 2*a - 4>=0, 2>=0, 2>=0, 0>=0, 2*a - 2>=0, k + 1>=0, - 2*a + k - 3>=0, - 2*a + k - 2>=0, - 2*a + k>=0, k - 7>=0, 2*a - k + 4>=0, 2*a - k + 5>=0, 2*a - k + 3>=0}$ ineq_solve(reg,{k,a}); {a=(2 .. infinity),k=2*a + 3} reg:= {a + b - c>=0, a - b + c>=0, - a + b + c>=0, 0>=0, 2>=0, 2*c - 2>=0, a - b + c>=0, a + b - c>=0, - a + b + c - 2>=0, 2>=0, 0>=0, 2*b - 2>=0, k + 1>=0, - a - b - c + k>=0, - a - b - c + k + 2>=0, - 2*b + k>=0, - 2*c + k>=0, a + b + c - k>=0, 2*b + 2*c - k - 2>=0, a + b + c - k>=0}$ ineq_solve (reg,{k,a,b,c}); {c=(1 .. infinity), b=(1 .. infinity), a=(max( - b + c,b - c) .. b + c - 2), k=a + b + c} clear reg; % Example from Richard Liska. lvars:={a,b,d}$ lfcond := {d>=0, b + d>=0, 2 a - b + d + 2>=0, - a + 2 d + 1>=0, b>=0, 2 a - b>=0, - a + 2 d>=0, b - d>=0, 2 a - b - d - 2>=0, - a + 2 d - 1>=0}$ ineq_solve(lfcond,lvars); {d=(2 .. infinity), b=(d .. 3*d - 4), b + d + 2 a=(----------- .. 2*d - 1)} 2 clear lfcond,lvars; end; Time for test: 30 ms, plus GC time: 10 ms