Solving an equation with constraintsSelect one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation valuesHow to plot a function with two parameters determined by another parameter?How to “Control” the Type of Solutions Mathematica Produces Using the Solve FunctionEquation Solving with NSolveNeed help setting up a numerical simulation and getting output into a tableGet random solution for equation with $n$ parameters and infinite solutions?Solving equation with NSolveGood practice when working with multiple solutions of a system of equationsSolving equation with constraintsUsing constraints when solving a trigometric equationSelect one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values

What does the same-ish mean?

How to coordinate airplane tickets?

Machine learning testing data

Is it a bad idea to plug the other end of ESD strap to wall ground?

Why are UK visa biometrics appointments suspended at USCIS Application Support Centers?

How could indestructible materials be used in power generation?

If a warlock makes a Dancing Sword their pact weapon, is there a way to prevent it from disappearing if it's farther away for more than a minute?

Processor speed limited at 0.4 Ghz

What are the G forces leaving Earth orbit?

Using "tail" to follow a file without displaying the most recent lines

Does Dispel Magic work on Tiny Hut?

One verb to replace 'be a member of' a club

Is there a hemisphere-neutral way of specifying a season?

Getting extremely large arrows with tikzcd

In the UK, is it possible to get a referendum by a court decision?

How can a day be of 24 hours?

Was the old ablative pronoun "med" or "mēd"?

GFCI outlets - can they be repaired? Are they really needed at the end of a circuit?

What historical events would have to change in order to make 19th century "steampunk" technology possible?

What is the opposite of "eschatology"?

Why do I get negative height?

Does int main() need a declaration on C++?

Why is the sentence "Das ist eine Nase" correct?

Do Iron Man suits sport waste management systems?



Solving an equation with constraints


Select one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation valuesHow to plot a function with two parameters determined by another parameter?How to “Control” the Type of Solutions Mathematica Produces Using the Solve FunctionEquation Solving with NSolveNeed help setting up a numerical simulation and getting output into a tableGet random solution for equation with $n$ parameters and infinite solutions?Solving equation with NSolveGood practice when working with multiple solutions of a system of equationsSolving equation with constraintsUsing constraints when solving a trigometric equationSelect one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values













2












$begingroup$


This is a question that follows up from this post:



Select one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values



I was trying to apply the answer to a more complicated case I am working on, and I'm getting an error message.



To start, here is the code of my more complicated function:



eq = 1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))) == 0


Among multiple solutions of $r$, I would like to select the positive real $r$ and simulate it with the following simulation values, for example, $d=0.8$, $s=2$, $t=0.02$ and $0<c<1$ and $1<q<2$. I would like to 3DPlot $r$ against $c$ and $q$.



My Mathematica code is as follows:



With[d = 0.8, s = 2, t = 0.02,
sol = Solve[eq, 0 <= c <= 1, 1 <= q <= 2, r, r > 0]];
Plot3D[Evaluate[r /. sol], c, 0, 1, q, 1, 2]


And once run, I get a very long warning message including something like:




Warning: r>0 is not a valid domain specification. Assuming it is a variable to eliminate.




Can anyone help?










share|improve this question











$endgroup$







  • 1




    $begingroup$
    You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
    $endgroup$
    – MarcoB
    4 hours ago










  • $begingroup$
    d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
    $endgroup$
    – Bill
    4 hours ago











  • $begingroup$
    @Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    @MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    (expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
    $endgroup$
    – Bill
    2 hours ago















2












$begingroup$


This is a question that follows up from this post:



Select one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values



I was trying to apply the answer to a more complicated case I am working on, and I'm getting an error message.



To start, here is the code of my more complicated function:



eq = 1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))) == 0


Among multiple solutions of $r$, I would like to select the positive real $r$ and simulate it with the following simulation values, for example, $d=0.8$, $s=2$, $t=0.02$ and $0<c<1$ and $1<q<2$. I would like to 3DPlot $r$ against $c$ and $q$.



My Mathematica code is as follows:



With[d = 0.8, s = 2, t = 0.02,
sol = Solve[eq, 0 <= c <= 1, 1 <= q <= 2, r, r > 0]];
Plot3D[Evaluate[r /. sol], c, 0, 1, q, 1, 2]


And once run, I get a very long warning message including something like:




Warning: r>0 is not a valid domain specification. Assuming it is a variable to eliminate.




Can anyone help?










share|improve this question











$endgroup$







  • 1




    $begingroup$
    You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
    $endgroup$
    – MarcoB
    4 hours ago










  • $begingroup$
    d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
    $endgroup$
    – Bill
    4 hours ago











  • $begingroup$
    @Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    @MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    (expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
    $endgroup$
    – Bill
    2 hours ago













2












2








2





$begingroup$


This is a question that follows up from this post:



Select one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values



I was trying to apply the answer to a more complicated case I am working on, and I'm getting an error message.



To start, here is the code of my more complicated function:



eq = 1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))) == 0


Among multiple solutions of $r$, I would like to select the positive real $r$ and simulate it with the following simulation values, for example, $d=0.8$, $s=2$, $t=0.02$ and $0<c<1$ and $1<q<2$. I would like to 3DPlot $r$ against $c$ and $q$.



My Mathematica code is as follows:



With[d = 0.8, s = 2, t = 0.02,
sol = Solve[eq, 0 <= c <= 1, 1 <= q <= 2, r, r > 0]];
Plot3D[Evaluate[r /. sol], c, 0, 1, q, 1, 2]


And once run, I get a very long warning message including something like:




Warning: r>0 is not a valid domain specification. Assuming it is a variable to eliminate.




Can anyone help?










share|improve this question











$endgroup$




This is a question that follows up from this post:



Select one among multiple solutions that satisfies a certain condition and 3D Plot it with varying simulation values



I was trying to apply the answer to a more complicated case I am working on, and I'm getting an error message.



To start, here is the code of my more complicated function:



eq = 1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))) == 0


Among multiple solutions of $r$, I would like to select the positive real $r$ and simulate it with the following simulation values, for example, $d=0.8$, $s=2$, $t=0.02$ and $0<c<1$ and $1<q<2$. I would like to 3DPlot $r$ against $c$ and $q$.



My Mathematica code is as follows:



With[d = 0.8, s = 2, t = 0.02,
sol = Solve[eq, 0 <= c <= 1, 1 <= q <= 2, r, r > 0]];
Plot3D[Evaluate[r /. sol], c, 0, 1, q, 1, 2]


And once run, I get a very long warning message including something like:




Warning: r>0 is not a valid domain specification. Assuming it is a variable to eliminate.




Can anyone help?







plotting equation-solving






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited 3 hours ago









m_goldberg

88.1k872199




88.1k872199










asked 5 hours ago









pppppp

306111




306111







  • 1




    $begingroup$
    You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
    $endgroup$
    – MarcoB
    4 hours ago










  • $begingroup$
    d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
    $endgroup$
    – Bill
    4 hours ago











  • $begingroup$
    @Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    @MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    (expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
    $endgroup$
    – Bill
    2 hours ago












  • 1




    $begingroup$
    You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
    $endgroup$
    – MarcoB
    4 hours ago










  • $begingroup$
    d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
    $endgroup$
    – Bill
    4 hours ago











  • $begingroup$
    @Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    @MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
    $endgroup$
    – ppp
    3 hours ago










  • $begingroup$
    (expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
    $endgroup$
    – Bill
    2 hours ago







1




1




$begingroup$
You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
$endgroup$
– MarcoB
4 hours ago




$begingroup$
You should write: Solve[equations, 0 <= c <= 1, 1 <= q <= 2, r > 0, r]. As the error suggests, r > 0 does not belong where you have it in your code.
$endgroup$
– MarcoB
4 hours ago












$begingroup$
d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
$endgroup$
– Bill
4 hours ago





$begingroup$
d=0.8;s=2;t=0.02;NMinimize[(expr)^2,0<c<1&&1<q<2,r,c,q] rapidly finds a solution near r->0.6457, c->0.1130, q->1.5452 Hopefully that will help you, but it isn't clear to me how you intend to Plot3D that.
$endgroup$
– Bill
4 hours ago













$begingroup$
@Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
$endgroup$
– ppp
3 hours ago




$begingroup$
@Bill: I'm trying to do simulation to see how the solution for r changes for different values of c and q. Would this be possible? In your code, what does (expr)^2 mean?
$endgroup$
– ppp
3 hours ago












$begingroup$
@MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
$endgroup$
– ppp
3 hours ago




$begingroup$
@MarcoB: Would it be possible to simulate (3D Plot) to see how the solution for r changes for different values of c and q? By the way, I'm trying your suggestion, and Mathematica has been running for quite a while. Did you get the answer quickly?
$endgroup$
– ppp
3 hours ago












$begingroup$
(expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
$endgroup$
– Bill
2 hours ago




$begingroup$
(expr)^2 was the square of your really long expression but with the trailing ==0 removed. So you wanted to find where your expression equaled zero and this found where the square of your expression was almost exactly zero. I am not sure there are still solutions for expr==0 if you make tiny changes in c or q or both. There may be some solution in r somewhere else, but not close by. I don't know. I'm guessing that is what you are supposed to learn from your exercise.
$endgroup$
– Bill
2 hours ago










1 Answer
1






active

oldest

votes


















3












$begingroup$

Mathematica's kernel seemed to slow down to a crawl when constraints on c and q were given along with your equation. However, the constraints are enforced in the s plot domain specification, so I thought it would worth trying to solve the equation without them which succeeded. Rationalization of the parameters suppressed all warnings.



Making the plot was a little slow because there seem some regions with very steep gradients that Plot3D has trouble with.



eq = (1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))));

Module[eqn, soln,
eqn = With[d = 4/5, s = 2, t = 2/100, Evaluate @ eq];
soln = Solve[Evaluate @ eqn == 0, r];
Plot3D[Evaluate[r /. soln], c, 0, 1, q, 1, 2, PlotPoints -> 50]]


plot



Edit



This addresses a issue raised by the OP in a comment below.



Here is way to get numerical values of r.



soln =
Module[eqn,
eqn = With[d = 0.8, s = 2, t = 0.02, Evaluate@eq];
NSolve[Evaluate @ eqn == 0, r]];
Do[rVal[i][c_, q_] = soln[[i, 1, 2]], i, Length[soln]];


Then to get the value of r for the 1st solution with c = .5 and `q = 1.5:



rVal[1][.5, 1.5]



-0.485047




To get all the value of r for all seven solutions at c = .25 and `q = 1.



Through[Array[rVal, 7][.25, 1.]]



-0.524898, 0.519435, 0.792331, 0.986375, 1.19151, 
0.0176256 - 0.0442494 I, 0.0176256 + 0.0442494 I



Notice that in this last case there are five real solutions.



2nd Edit



Here is an exploration of the seven root objects. Note that two of plots are are empty because two the root objects are never real in the specified domain. Also, you should be able to see where the icicles come from.



Column[
Table[
Plot3D[Evaluate[rVal[i][c, q]], c, 0, 1, q, 1, 2, PlotPoints -> 50],
i, 7]]


plots






share|improve this answer











$endgroup$












  • $begingroup$
    Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
    $endgroup$
    – ppp
    2 hours ago











  • $begingroup$
    @ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
    $endgroup$
    – m_goldberg
    2 hours ago










  • $begingroup$
    @ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
    $endgroup$
    – m_goldberg
    1 hour ago










  • $begingroup$
    Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
    $endgroup$
    – ppp
    1 hour ago










  • $begingroup$
    One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
    $endgroup$
    – ppp
    1 hour ago











Your Answer





StackExchange.ifUsing("editor", function ()
return StackExchange.using("mathjaxEditing", function ()
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
);
);
, "mathjax-editing");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "387"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);













draft saved

draft discarded


















StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f194463%2fsolving-an-equation-with-constraints%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown

























1 Answer
1






active

oldest

votes








1 Answer
1






active

oldest

votes









active

oldest

votes






active

oldest

votes









3












$begingroup$

Mathematica's kernel seemed to slow down to a crawl when constraints on c and q were given along with your equation. However, the constraints are enforced in the s plot domain specification, so I thought it would worth trying to solve the equation without them which succeeded. Rationalization of the parameters suppressed all warnings.



Making the plot was a little slow because there seem some regions with very steep gradients that Plot3D has trouble with.



eq = (1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))));

Module[eqn, soln,
eqn = With[d = 4/5, s = 2, t = 2/100, Evaluate @ eq];
soln = Solve[Evaluate @ eqn == 0, r];
Plot3D[Evaluate[r /. soln], c, 0, 1, q, 1, 2, PlotPoints -> 50]]


plot



Edit



This addresses a issue raised by the OP in a comment below.



Here is way to get numerical values of r.



soln =
Module[eqn,
eqn = With[d = 0.8, s = 2, t = 0.02, Evaluate@eq];
NSolve[Evaluate @ eqn == 0, r]];
Do[rVal[i][c_, q_] = soln[[i, 1, 2]], i, Length[soln]];


Then to get the value of r for the 1st solution with c = .5 and `q = 1.5:



rVal[1][.5, 1.5]



-0.485047




To get all the value of r for all seven solutions at c = .25 and `q = 1.



Through[Array[rVal, 7][.25, 1.]]



-0.524898, 0.519435, 0.792331, 0.986375, 1.19151, 
0.0176256 - 0.0442494 I, 0.0176256 + 0.0442494 I



Notice that in this last case there are five real solutions.



2nd Edit



Here is an exploration of the seven root objects. Note that two of plots are are empty because two the root objects are never real in the specified domain. Also, you should be able to see where the icicles come from.



Column[
Table[
Plot3D[Evaluate[rVal[i][c, q]], c, 0, 1, q, 1, 2, PlotPoints -> 50],
i, 7]]


plots






share|improve this answer











$endgroup$












  • $begingroup$
    Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
    $endgroup$
    – ppp
    2 hours ago











  • $begingroup$
    @ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
    $endgroup$
    – m_goldberg
    2 hours ago










  • $begingroup$
    @ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
    $endgroup$
    – m_goldberg
    1 hour ago










  • $begingroup$
    Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
    $endgroup$
    – ppp
    1 hour ago










  • $begingroup$
    One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
    $endgroup$
    – ppp
    1 hour ago















3












$begingroup$

Mathematica's kernel seemed to slow down to a crawl when constraints on c and q were given along with your equation. However, the constraints are enforced in the s plot domain specification, so I thought it would worth trying to solve the equation without them which succeeded. Rationalization of the parameters suppressed all warnings.



Making the plot was a little slow because there seem some regions with very steep gradients that Plot3D has trouble with.



eq = (1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))));

Module[eqn, soln,
eqn = With[d = 4/5, s = 2, t = 2/100, Evaluate @ eq];
soln = Solve[Evaluate @ eqn == 0, r];
Plot3D[Evaluate[r /. soln], c, 0, 1, q, 1, 2, PlotPoints -> 50]]


plot



Edit



This addresses a issue raised by the OP in a comment below.



Here is way to get numerical values of r.



soln =
Module[eqn,
eqn = With[d = 0.8, s = 2, t = 0.02, Evaluate@eq];
NSolve[Evaluate @ eqn == 0, r]];
Do[rVal[i][c_, q_] = soln[[i, 1, 2]], i, Length[soln]];


Then to get the value of r for the 1st solution with c = .5 and `q = 1.5:



rVal[1][.5, 1.5]



-0.485047




To get all the value of r for all seven solutions at c = .25 and `q = 1.



Through[Array[rVal, 7][.25, 1.]]



-0.524898, 0.519435, 0.792331, 0.986375, 1.19151, 
0.0176256 - 0.0442494 I, 0.0176256 + 0.0442494 I



Notice that in this last case there are five real solutions.



2nd Edit



Here is an exploration of the seven root objects. Note that two of plots are are empty because two the root objects are never real in the specified domain. Also, you should be able to see where the icicles come from.



Column[
Table[
Plot3D[Evaluate[rVal[i][c, q]], c, 0, 1, q, 1, 2, PlotPoints -> 50],
i, 7]]


plots






share|improve this answer











$endgroup$












  • $begingroup$
    Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
    $endgroup$
    – ppp
    2 hours ago











  • $begingroup$
    @ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
    $endgroup$
    – m_goldberg
    2 hours ago










  • $begingroup$
    @ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
    $endgroup$
    – m_goldberg
    1 hour ago










  • $begingroup$
    Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
    $endgroup$
    – ppp
    1 hour ago










  • $begingroup$
    One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
    $endgroup$
    – ppp
    1 hour ago













3












3








3





$begingroup$

Mathematica's kernel seemed to slow down to a crawl when constraints on c and q were given along with your equation. However, the constraints are enforced in the s plot domain specification, so I thought it would worth trying to solve the equation without them which succeeded. Rationalization of the parameters suppressed all warnings.



Making the plot was a little slow because there seem some regions with very steep gradients that Plot3D has trouble with.



eq = (1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))));

Module[eqn, soln,
eqn = With[d = 4/5, s = 2, t = 2/100, Evaluate @ eq];
soln = Solve[Evaluate @ eqn == 0, r];
Plot3D[Evaluate[r /. soln], c, 0, 1, q, 1, 2, PlotPoints -> 50]]


plot



Edit



This addresses a issue raised by the OP in a comment below.



Here is way to get numerical values of r.



soln =
Module[eqn,
eqn = With[d = 0.8, s = 2, t = 0.02, Evaluate@eq];
NSolve[Evaluate @ eqn == 0, r]];
Do[rVal[i][c_, q_] = soln[[i, 1, 2]], i, Length[soln]];


Then to get the value of r for the 1st solution with c = .5 and `q = 1.5:



rVal[1][.5, 1.5]



-0.485047




To get all the value of r for all seven solutions at c = .25 and `q = 1.



Through[Array[rVal, 7][.25, 1.]]



-0.524898, 0.519435, 0.792331, 0.986375, 1.19151, 
0.0176256 - 0.0442494 I, 0.0176256 + 0.0442494 I



Notice that in this last case there are five real solutions.



2nd Edit



Here is an exploration of the seven root objects. Note that two of plots are are empty because two the root objects are never real in the specified domain. Also, you should be able to see where the icicles come from.



Column[
Table[
Plot3D[Evaluate[rVal[i][c, q]], c, 0, 1, q, 1, 2, PlotPoints -> 50],
i, 7]]


plots






share|improve this answer











$endgroup$



Mathematica's kernel seemed to slow down to a crawl when constraints on c and q were given along with your equation. However, the constraints are enforced in the s plot domain specification, so I thought it would worth trying to solve the equation without them which succeeded. Rationalization of the parameters suppressed all warnings.



Making the plot was a little slow because there seem some regions with very steep gradients that Plot3D has trouble with.



eq = (1/(24 (-1 + r)^3 r^4 s) (3 d^6 (-3 + 5 r) (-1 + t) + 6 d^5 (3 - 5 r + c (-3 + r (9 - 4 q + 8 (-1 + q) r))) (-1 + t) + 12 d (-1 + r)^3 r^2 (s + 2 t) - 4 (-1 + r)^3 r^2 s ((-1 + r^2) s (-1 + t) + 3 t) + 12 d^3 (-1 + r) r (2 + c (-2 + r (5 - 2 q + 4 (-1 + q) r)) + s - r (3 + s + r s) + 2 t - 3 r t + (-1 + r + r^2) s t) + 3 d^4 (3 + c^2 (1 + 2 (-1 + q) r) (-3 + r (7 - 2 q + 6 (-1 + q) r)) (-1 + t) - 2 c (-3 + r (9 - 4 q + 8 (-1 + q) r)) (-1 + t) - 3 t + r (3 + 4 r (-5 + 3 r) + 5 t)) + 12 d^2 (-1 + r) r (-(-1 + c) (s (-1 + t) - 2 t) + r^3 (-1 + c (-1 + q) s (-1 + t) - t) + r^2 (2 + (-1 + c (-1 + 2 q)) s (-1 + t) + (2 + 4 c - 4 c q) t) + r (-1 - (1 + c (-2 + q)) s (-1 + t) + (2 - 5 c + 2 c q) t))));

Module[eqn, soln,
eqn = With[d = 4/5, s = 2, t = 2/100, Evaluate @ eq];
soln = Solve[Evaluate @ eqn == 0, r];
Plot3D[Evaluate[r /. soln], c, 0, 1, q, 1, 2, PlotPoints -> 50]]


plot



Edit



This addresses a issue raised by the OP in a comment below.



Here is way to get numerical values of r.



soln =
Module[eqn,
eqn = With[d = 0.8, s = 2, t = 0.02, Evaluate@eq];
NSolve[Evaluate @ eqn == 0, r]];
Do[rVal[i][c_, q_] = soln[[i, 1, 2]], i, Length[soln]];


Then to get the value of r for the 1st solution with c = .5 and `q = 1.5:



rVal[1][.5, 1.5]



-0.485047




To get all the value of r for all seven solutions at c = .25 and `q = 1.



Through[Array[rVal, 7][.25, 1.]]



-0.524898, 0.519435, 0.792331, 0.986375, 1.19151, 
0.0176256 - 0.0442494 I, 0.0176256 + 0.0442494 I



Notice that in this last case there are five real solutions.



2nd Edit



Here is an exploration of the seven root objects. Note that two of plots are are empty because two the root objects are never real in the specified domain. Also, you should be able to see where the icicles come from.



Column[
Table[
Plot3D[Evaluate[rVal[i][c, q]], c, 0, 1, q, 1, 2, PlotPoints -> 50],
i, 7]]


plots







share|improve this answer














share|improve this answer



share|improve this answer








edited 16 mins ago

























answered 2 hours ago









m_goldbergm_goldberg

88.1k872199




88.1k872199











  • $begingroup$
    Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
    $endgroup$
    – ppp
    2 hours ago











  • $begingroup$
    @ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
    $endgroup$
    – m_goldberg
    2 hours ago










  • $begingroup$
    @ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
    $endgroup$
    – m_goldberg
    1 hour ago










  • $begingroup$
    Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
    $endgroup$
    – ppp
    1 hour ago










  • $begingroup$
    One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
    $endgroup$
    – ppp
    1 hour ago
















  • $begingroup$
    Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
    $endgroup$
    – ppp
    2 hours ago











  • $begingroup$
    @ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
    $endgroup$
    – m_goldberg
    2 hours ago










  • $begingroup$
    @ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
    $endgroup$
    – m_goldberg
    1 hour ago










  • $begingroup$
    Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
    $endgroup$
    – ppp
    1 hour ago










  • $begingroup$
    One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
    $endgroup$
    – ppp
    1 hour ago















$begingroup$
Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
$endgroup$
– ppp
2 hours ago





$begingroup$
Thanks so much! Can you please give me some comments on how to interpret the graphical result? It looks very weird... Why it has layers? What are those icicle-like things? Why change in the color, say from green to purple? Etc... From the graph, what is the approximate value of r for, say, c=0.2 & q=1.2? Isn't it impossible to tell since there are three corresponding values of r? Green, blue, and orange?
$endgroup$
– ppp
2 hours ago













$begingroup$
@ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
$endgroup$
– m_goldberg
2 hours ago




$begingroup$
@ppp. Solve finds seven root functions for your equation, so at any patch of the plot domain, up to seven surfaces could show up. But it seems that only three of the root functions actually evaluate to real values over any significant part of the domain. The icicle-like things are areas where the solver found small patches, evidently with very steep gradients, of a real surface for some of the other possible surfaces.There is also a smallish patch of a red surface which can only be seen when the plot is rotated to different view.
$endgroup$
– m_goldberg
2 hours ago












$begingroup$
@ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
$endgroup$
– m_goldberg
1 hour ago




$begingroup$
@ppp. You can can explore each of the seven root objects returned from Solve individually to find the one that best represents a solution to your problem. Further, you can get a value of r for any pair of c and q from the individual root objects.
$endgroup$
– m_goldberg
1 hour ago












$begingroup$
Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
$endgroup$
– ppp
1 hour ago




$begingroup$
Do you know how to 'explore each of the root objects returned from Solve individually to find the one that best represents a solution to my problem'? By the way, I have restricted the roots to 0<=r<=1 by adding PlotRange->0,1. It would be desirable to get a unique positive root lying between 0 and 1.
$endgroup$
– ppp
1 hour ago












$begingroup$
One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
$endgroup$
– ppp
1 hour ago




$begingroup$
One more question: I'm sorry but I'm not quite understanding your comments about the icicle-like things. May I ask your explanation once more? I really appreciate!
$endgroup$
– ppp
1 hour ago

















draft saved

draft discarded
















































Thanks for contributing an answer to Mathematica Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid


  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.

Use MathJax to format equations. MathJax reference.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f194463%2fsolving-an-equation-with-constraints%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Can not update quote_id field of “quote_item” table magento 2Magento 2.1 - We can't remove the item. (Shopping Cart doesnt allow us to remove items before becomes empty)Add value for custom quote item attribute using REST apiREST API endpoint v1/carts/cartId/items always returns error messageCorrect way to save entries to databaseHow to remove all associated quote objects of a customer completelyMagento 2 - Save value from custom input field to quote_itemGet quote_item data using quote id and product id filter in Magento 2How to set additional data to quote_item table from controller in Magento 2?What is the purpose of additional_data column in quote_item table in magento2Set Custom Price to Quote item magento2 from controller

How to solve knockout JS error in Magento 2 Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?(Magento2) knockout.js:3012 Uncaught ReferenceError: Unable to process bindingUnable to process binding Knockout.js magento 2Cannot read property `scopeLabel` of undefined on Product Detail PageCan't get Customer Data on frontend in Magento 2Magento2 Order Summary - unable to process bindingKO templates are not loading in Magento 2.1 applicationgetting knockout js error magento 2Product grid not load -— Unable to process binding Knockout.js magento 2Product form not loaded in magento2Uncaught ReferenceError: Unable to process binding “if: function()return (isShowLegend()) ” magento 2

Nissan Patrol Зміст Перше покоління — 4W60 (1951-1960) | Друге покоління — 60 series (1960-1980) | Третє покоління (1980–2002) | Четверте покоління — Y60 (1987–1998) | П'яте покоління — Y61 (1997–2013) | Шосте покоління — Y62 (2010- ) | Посилання | Зноски | Навігаційне менюОфіційний український сайтТест-драйв Nissan Patrol 2010 7-го поколінняNissan PatrolКак мы тестировали Nissan Patrol 2016рвиправивши або дописавши її