Mathematica MathematicaSimultaneous Inequalities?

  • Thread starter Thread starter Saladsamurai
  • Start date Start date
  • Tags Tags
    Inequalities
AI Thread Summary
The discussion revolves around finding the range of values for K such that three inequalities involving variables A, B, and C, all greater than zero, are satisfied. Users suggest using Mathematica's built-in functions like Reduce to determine the limits for K. The inequalities are defined as A = 576 - 11819*K, B = 15848*K, and C = 397119*K - 33.6*B/A. The results indicate that K must be in the range 0 < K < 0.0486216 for all three variables to remain positive. The conversation highlights the effectiveness of Mathematica in solving simultaneous inequalities, despite some complexity in the output.
Saladsamurai
Messages
3,009
Reaction score
7
Hey guys,

I have some variables A,B, and C such that all must be greater then zero. A,B and C are all functions of K. I would like to find the range of values of K such that the 3 inequalities are satisfied.

Is there a function built into Mathematica that will do this? I could write a code, but I am kind of sick of coding right now :-p

If not..then code it is!
 
Physics news on Phys.org
Well, as a quick response I know Reduce can sort of find the limits if you know what you're looking for.

Code:
a[k_] = (-1) k^3 + k^2 + 6;
b[k_] = Sin[k];
c[k_] = Cos[k];
Plot[{a[k], b[k], c[k]}, {k, 0, 3}]
Reduce[a[k] > 0 && b[k] > 0 && c[k] > 0, k] // FullSimplify

So the result :
c1 elem of Z OR c1 <= 0 OR 0 <k - 2 pi c1 < pi/2

Means its true if (looking at the second two) if c1(some constant) is negative, so that

k - 2 pi n < pi/2
basically the first solution being
0< k < pi/2 (which is the region that's positive where all three are positive).

It's not perfect, but its correct.
 
Hey there Hepth,

this is what I have:

A = 576 - 11819*K
B = 15848*K
C = C = 397119*K - 33.6*B/A

and I need to find the range of K such that A,B, and C all remain greater than 0.
 
A = 576 - 11819*k;
B = 15848*k;
CC = 397119*k - 33.6*B/A;
Plot[{A, B, CC}, {k, 0, 0.04862}, PlotRange -> {0, 1000}]
Reduce[A > 0 && B > 0 && CC > 0, k] // FullSimplify
Reduce[A > y && B > y && CC > y && y > 0, y] // FullSimplify

Gives for k : 0<k<0.0486216 (WHERE they're all positive)
And for Y (bounded by this region):
<br /> y&gt;0\land \left((k&gt;0\land y&lt;15848. k\land k\leq 0.020819)\lor (k&gt;0.020819\land 11819. k+y&lt;576.\land k\leq 0.0486216)\lor \left(k&gt;0.0486216\land
<br /> \landk&lt;0.0486216\land y&lt;\frac{k (2.34677\times 10^{10} k-1.14104\times 10^9)}{59095. k-2880.}\right)\right)<br />

Which, while ugly, describes that range and domain bounded by the functions. I think you're only asking for the first one though.
 
Last edited:
That is great Hepth! I think that will do the trick :smile: Thanks for your help.
 

Similar threads

Replies
5
Views
3K
Replies
1
Views
1K
Replies
2
Views
2K
Replies
5
Views
2K
Replies
7
Views
2K
Replies
13
Views
2K
Back
Top