Possible bug in Solve function?












15














Bug introduced in 10.0 and persisting through 11.3 or later




In 11.3.0 for Microsoft Windows (64-bit) (March 7, 2018) writing:



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)



eqn = {D[f[w, x, y, z], w] == 0, 

       D[f[w, x, y, z], x] == 0, 

       D[f[w, x, y, z], y] == 0, 

       D[f[w, x, y, z], z] == 0};



sol = Solve[eqn];



Table[eqn /. sol[[n]], {n, Length[sol]}]


I get:




{{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True}}




from which there are four wrong solutions.



Am I wrong or is it a Solve bug?




EDIT: through the email address support@wolfram.com I contacted Wolfram Technical Support who in less than three working days have confirmed that it is a bug and have already proceeded to report to their developers.






equation-solving bugs







share|improve this question




















  • 1




    You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
    – Szabolcs
    Jan 5 at 22:30










  • Select[sol, And @@ eqn /. # &]
    – Bob Hanlon
    Jan 5 at 23:53






  • 5




    Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
    – J. M. is computer-less♦
    Jan 6 at 2:45
















15














Bug introduced in 10.0 and persisting through 11.3 or later




In 11.3.0 for Microsoft Windows (64-bit) (March 7, 2018) writing:



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)



eqn = {D[f[w, x, y, z], w] == 0, 

       D[f[w, x, y, z], x] == 0, 

       D[f[w, x, y, z], y] == 0, 

       D[f[w, x, y, z], z] == 0};



sol = Solve[eqn];



Table[eqn /. sol[[n]], {n, Length[sol]}]


I get:




{{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True}}




from which there are four wrong solutions.



Am I wrong or is it a Solve bug?




EDIT: through the email address support@wolfram.com I contacted Wolfram Technical Support who in less than three working days have confirmed that it is a bug and have already proceeded to report to their developers.






equation-solving bugs







share|improve this question




















  • 1




    You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
    – Szabolcs
    Jan 5 at 22:30










  • Select[sol, And @@ eqn /. # &]
    – Bob Hanlon
    Jan 5 at 23:53






  • 5




    Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
    – J. M. is computer-less♦
    Jan 6 at 2:45














15












15








15


2





Bug introduced in 10.0 and persisting through 11.3 or later




In 11.3.0 for Microsoft Windows (64-bit) (March 7, 2018) writing:



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)



eqn = {D[f[w, x, y, z], w] == 0, 

       D[f[w, x, y, z], x] == 0, 

       D[f[w, x, y, z], y] == 0, 

       D[f[w, x, y, z], z] == 0};



sol = Solve[eqn];



Table[eqn /. sol[[n]], {n, Length[sol]}]


I get:




{{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True}}




from which there are four wrong solutions.



Am I wrong or is it a Solve bug?




EDIT: through the email address support@wolfram.com I contacted Wolfram Technical Support who in less than three working days have confirmed that it is a bug and have already proceeded to report to their developers.






equation-solving bugs







share|improve this question















Bug introduced in 10.0 and persisting through 11.3 or later




In 11.3.0 for Microsoft Windows (64-bit) (March 7, 2018) writing:



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)



eqn = {D[f[w, x, y, z], w] == 0, 

       D[f[w, x, y, z], x] == 0, 

       D[f[w, x, y, z], y] == 0, 

       D[f[w, x, y, z], z] == 0};



sol = Solve[eqn];



Table[eqn /. sol[[n]], {n, Length[sol]}]


I get:




{{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True},
{False, True, True, False},
{True, True, True, True},
{True, True, True, True}}




from which there are four wrong solutions.



Am I wrong or is it a Solve bug?




EDIT: through the email address support@wolfram.com I contacted Wolfram Technical Support who in less than three working days have confirmed that it is a bug and have already proceeded to report to their developers.






equation-solving bugs




equation-solving bugs






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited 15 hours ago









Szabolcs

158k13432928




158k13432928










asked Jan 5 at 21:50









TeMTeM

1,990621




1,990621








  • 1




    You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
    – Szabolcs
    Jan 5 at 22:30










  • Select[sol, And @@ eqn /. # &]
    – Bob Hanlon
    Jan 5 at 23:53






  • 5




    Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
    – J. M. is computer-less♦
    Jan 6 at 2:45














  • 1




    You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
    – Szabolcs
    Jan 5 at 22:30










  • Select[sol, And @@ eqn /. # &]
    – Bob Hanlon
    Jan 5 at 23:53






  • 5




    Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
    – J. M. is computer-less♦
    Jan 6 at 2:45








1




1




You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
– Szabolcs
Jan 5 at 22:30




You could use List@ToRules@Reduce[eqn, {x, y, z, w}] to get all valid solutions. Filter for those that only have numeric values on the RHS of ->.
– Szabolcs
Jan 5 at 22:30












Select[sol, And @@ eqn /. # &]
– Bob Hanlon
Jan 5 at 23:53




Select[sol, And @@ eqn /. # &]
– Bob Hanlon
Jan 5 at 23:53




5




5




Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
– J. M. is computer-less♦
Jan 6 at 2:45




Next time, please do not add the bugs tag yourself on a question. The tag is only supposed to be added after your observations have been confirmed by other users.
– J. M. is computer-less♦
Jan 6 at 2:45










2 Answers
2






active

oldest

votes


















5














Thanks for asking! In version 9.0 only 16 solutions are returned and they are all valid. In version 10.2 there are 20 solutions, with the extra 4 all being invalid. Contragulations! I think you found a bug. You may want to click "Help", then "Give Feedback...", and then fill out the form in your browser to report.



As another answer notes, you can always try Reduce instead which may give better results in some cases, but Solve is usually what you want.






share|improve this answer























  • Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
    – TeM
    Jan 5 at 22:11






  • 2




    @TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
    – Szabolcs
    Jan 5 at 22:31










  • @Szabolcs: Could you direct me where I can do it correctly?
    – TeM
    Jan 5 at 22:44






  • 4




    @TeM wolfram.com/support/contact
    – Szabolcs
    Jan 5 at 23:49



















4














You can use Reduce



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)

eqn = {D[f[w, x, y, z], w] == 0, D[f[w, x, y, z], x] == 0, 

   D[f[w, x, y, z], y] == 0, D[f[w, x, y, z], z] == 0};

red = Reduce[eqn, Backsubstitution -> True]



$left(z=0land x=0land w=-sqrt{1-y^2}right)lor left(z=0land x=0land w=sqrt{1-y^2}right)lor left(z=0land y=0land
w=-sqrt{1-x^2}right)lor left(z=0land y=0land w=sqrt{1-x^2}right)lor (z=0land y=-1land x=0land w=0)lor (z=0land y=0land x=0land
w=-1)lor (z=0land y=0land x=0land w=1)lor (z=0land y=1land x=0land w=0)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land
x=-frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)\
lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)$




First@eqn //. {ToRules[red]}



{True, True, True, True, True, True, True, True, True, True, True,
True, True, True, True, True}







share|improve this answer





















  • Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
    – TeM
    Jan 6 at 1:14






  • 4




    I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
    – Okkes Dulgerci
    Jan 6 at 1:18











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%2f188900%2fpossible-bug-in-solve-function%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























2 Answers
2






active

oldest

votes








2 Answers
2






active

oldest

votes









active

oldest

votes






active

oldest

votes









5














Thanks for asking! In version 9.0 only 16 solutions are returned and they are all valid. In version 10.2 there are 20 solutions, with the extra 4 all being invalid. Contragulations! I think you found a bug. You may want to click "Help", then "Give Feedback...", and then fill out the form in your browser to report.



As another answer notes, you can always try Reduce instead which may give better results in some cases, but Solve is usually what you want.






share|improve this answer























  • Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
    – TeM
    Jan 5 at 22:11






  • 2




    @TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
    – Szabolcs
    Jan 5 at 22:31










  • @Szabolcs: Could you direct me where I can do it correctly?
    – TeM
    Jan 5 at 22:44






  • 4




    @TeM wolfram.com/support/contact
    – Szabolcs
    Jan 5 at 23:49
















5














Thanks for asking! In version 9.0 only 16 solutions are returned and they are all valid. In version 10.2 there are 20 solutions, with the extra 4 all being invalid. Contragulations! I think you found a bug. You may want to click "Help", then "Give Feedback...", and then fill out the form in your browser to report.



As another answer notes, you can always try Reduce instead which may give better results in some cases, but Solve is usually what you want.






share|improve this answer























  • Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
    – TeM
    Jan 5 at 22:11






  • 2




    @TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
    – Szabolcs
    Jan 5 at 22:31










  • @Szabolcs: Could you direct me where I can do it correctly?
    – TeM
    Jan 5 at 22:44






  • 4




    @TeM wolfram.com/support/contact
    – Szabolcs
    Jan 5 at 23:49














5












5








5






Thanks for asking! In version 9.0 only 16 solutions are returned and they are all valid. In version 10.2 there are 20 solutions, with the extra 4 all being invalid. Contragulations! I think you found a bug. You may want to click "Help", then "Give Feedback...", and then fill out the form in your browser to report.



As another answer notes, you can always try Reduce instead which may give better results in some cases, but Solve is usually what you want.






share|improve this answer














Thanks for asking! In version 9.0 only 16 solutions are returned and they are all valid. In version 10.2 there are 20 solutions, with the extra 4 all being invalid. Contragulations! I think you found a bug. You may want to click "Help", then "Give Feedback...", and then fill out the form in your browser to report.



As another answer notes, you can always try Reduce instead which may give better results in some cases, but Solve is usually what you want.







share|improve this answer














share|improve this answer



share|improve this answer








edited Jan 6 at 1:16

























answered Jan 5 at 22:05









SomosSomos

3628




3628












  • Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
    – TeM
    Jan 5 at 22:11






  • 2




    @TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
    – Szabolcs
    Jan 5 at 22:31










  • @Szabolcs: Could you direct me where I can do it correctly?
    – TeM
    Jan 5 at 22:44






  • 4




    @TeM wolfram.com/support/contact
    – Szabolcs
    Jan 5 at 23:49


















  • Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
    – TeM
    Jan 5 at 22:11






  • 2




    @TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
    – Szabolcs
    Jan 5 at 22:31










  • @Szabolcs: Could you direct me where I can do it correctly?
    – TeM
    Jan 5 at 22:44






  • 4




    @TeM wolfram.com/support/contact
    – Szabolcs
    Jan 5 at 23:49
















Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
– TeM
Jan 5 at 22:11




Well, it will mean that the next version will be the first calculation I will perform. Thank you! ^_^
– TeM
Jan 5 at 22:11




2




2




@TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
– Szabolcs
Jan 5 at 22:31




@TeM You should report it to Wolfram first, otherwise there's no chance for it to get fixed.
– Szabolcs
Jan 5 at 22:31












@Szabolcs: Could you direct me where I can do it correctly?
– TeM
Jan 5 at 22:44




@Szabolcs: Could you direct me where I can do it correctly?
– TeM
Jan 5 at 22:44




4




4




@TeM wolfram.com/support/contact
– Szabolcs
Jan 5 at 23:49




@TeM wolfram.com/support/contact
– Szabolcs
Jan 5 at 23:49











4














You can use Reduce



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)

eqn = {D[f[w, x, y, z], w] == 0, D[f[w, x, y, z], x] == 0, 

   D[f[w, x, y, z], y] == 0, D[f[w, x, y, z], z] == 0};

red = Reduce[eqn, Backsubstitution -> True]



$left(z=0land x=0land w=-sqrt{1-y^2}right)lor left(z=0land x=0land w=sqrt{1-y^2}right)lor left(z=0land y=0land
w=-sqrt{1-x^2}right)lor left(z=0land y=0land w=sqrt{1-x^2}right)lor (z=0land y=-1land x=0land w=0)lor (z=0land y=0land x=0land
w=-1)lor (z=0land y=0land x=0land w=1)lor (z=0land y=1land x=0land w=0)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land
x=-frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)\
lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)$




First@eqn //. {ToRules[red]}



{True, True, True, True, True, True, True, True, True, True, True,
True, True, True, True, True}







share|improve this answer





















  • Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
    – TeM
    Jan 6 at 1:14






  • 4




    I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
    – Okkes Dulgerci
    Jan 6 at 1:18
















4














You can use Reduce



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)

eqn = {D[f[w, x, y, z], w] == 0, D[f[w, x, y, z], x] == 0, 

   D[f[w, x, y, z], y] == 0, D[f[w, x, y, z], z] == 0};

red = Reduce[eqn, Backsubstitution -> True]



$left(z=0land x=0land w=-sqrt{1-y^2}right)lor left(z=0land x=0land w=sqrt{1-y^2}right)lor left(z=0land y=0land
w=-sqrt{1-x^2}right)lor left(z=0land y=0land w=sqrt{1-x^2}right)lor (z=0land y=-1land x=0land w=0)lor (z=0land y=0land x=0land
w=-1)lor (z=0land y=0land x=0land w=1)lor (z=0land y=1land x=0land w=0)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land
x=-frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)\
lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)$




First@eqn //. {ToRules[red]}



{True, True, True, True, True, True, True, True, True, True, True,
True, True, True, True, True}







share|improve this answer





















  • Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
    – TeM
    Jan 6 at 1:14






  • 4




    I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
    – Okkes Dulgerci
    Jan 6 at 1:18














4












4








4






You can use Reduce



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)

eqn = {D[f[w, x, y, z], w] == 0, D[f[w, x, y, z], x] == 0, 

   D[f[w, x, y, z], y] == 0, D[f[w, x, y, z], z] == 0};

red = Reduce[eqn, Backsubstitution -> True]



$left(z=0land x=0land w=-sqrt{1-y^2}right)lor left(z=0land x=0land w=sqrt{1-y^2}right)lor left(z=0land y=0land
w=-sqrt{1-x^2}right)lor left(z=0land y=0land w=sqrt{1-x^2}right)lor (z=0land y=-1land x=0land w=0)lor (z=0land y=0land x=0land
w=-1)lor (z=0land y=0land x=0land w=1)lor (z=0land y=1land x=0land w=0)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land
x=-frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)\
lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)$




First@eqn //. {ToRules[red]}



{True, True, True, True, True, True, True, True, True, True, True,
True, True, True, True, True}







share|improve this answer












You can use Reduce



f[w_, x_, y_, z_] := w*x^2*y^3 - z*(w^2 + x^2 + y^2 - 1)

eqn = {D[f[w, x, y, z], w] == 0, D[f[w, x, y, z], x] == 0, 

   D[f[w, x, y, z], y] == 0, D[f[w, x, y, z], z] == 0};

red = Reduce[eqn, Backsubstitution -> True]



$left(z=0land x=0land w=-sqrt{1-y^2}right)lor left(z=0land x=0land w=sqrt{1-y^2}right)lor left(z=0land y=0land
w=-sqrt{1-x^2}right)lor left(z=0land y=0land w=sqrt{1-x^2}right)lor (z=0land y=-1land x=0land w=0)lor (z=0land y=0land x=0land
w=-1)lor (z=0land y=0land x=0land w=1)lor (z=0land y=1land x=0land w=0)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land
x=-frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=-frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=-frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land
w=-frac{1}{sqrt{6}}right)lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=-frac{1}{sqrt{3}}land
w=frac{1}{sqrt{6}}right)\
lor left(z=frac{1}{4 sqrt{3}}land y=frac{1}{sqrt{2}}land x=frac{1}{sqrt{3}}land w=frac{1}{sqrt{6}}right)$




First@eqn //. {ToRules[red]}



{True, True, True, True, True, True, True, True, True, True, True,
True, True, True, True, True}








share|improve this answer












share|improve this answer



share|improve this answer










answered Jan 6 at 1:09









Okkes DulgerciOkkes Dulgerci

4,1851816




4,1851816












  • Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
    – TeM
    Jan 6 at 1:14






  • 4




    I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
    – Okkes Dulgerci
    Jan 6 at 1:18


















  • Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
    – TeM
    Jan 6 at 1:14






  • 4




    I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
    – Okkes Dulgerci
    Jan 6 at 1:18
















Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
– TeM
Jan 6 at 1:14




Yes, of course, I had already tried. I was almost certain it was a bug from Solve, so I pointed out. Thank you very much anyway, always very kind!
– TeM
Jan 6 at 1:14




4




4




I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
– Okkes Dulgerci
Jan 6 at 1:18




I believe Solve use the function Reduce under the hood. when you remove Backsubstitution -> True, you'll find implicit solution, somehow Solve messes up somewhere and it is definitely a bug..
– Okkes Dulgerci
Jan 6 at 1:18


















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.





Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


Please pay close attention to the following guidance:


  • 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.


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%2f188900%2fpossible-bug-in-solve-function%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

If I really need a card on my start hand, how many mulligans make sense? [duplicate]

Image
5 This question already has an answer here: How do you calculate the likelihood of drawing certain cards in your opening hand? 3 answers There is a funny modern deck where you play Treasure Hunt, but you really need to have it on the starting hand. How many mulligans should I take, if my deck contains 60 cards and contains 4 copies of Treasure Hunt? magic-the-gathering probability share | improve this question asked Nov 23 '18 at 12:12 dan dan 152 4
Read more

Alcedinidae

Image
Alcedinidae Martin-pêcheur à dos bleu ( Alcedo azurea ) Classification (COI) Règne Animalia Embranchement Chordata Sous-embr. Vertebrata Classe Aves Ordre Coraciiformes Famille Alcedinidae Rafinesque, 1815 Les Alcédinidés ou Alcedinidae forment une famille d'oiseaux plus connus sous les noms de martins-pêcheurs et martins-chasseurs. Sommaire 1 Description 2 Répartition et habitat 3 Systématique 3.1 Liste alphabétique des genres 3.2 Liste des espèces 4 Chez les Anciens : l'alcyon 4.1 Étymologie (Littré) 4.2 Oiseau « fabuleux » ? 5 Notes 6 Articles connexes 7 Liens externes 7.1 Bibliographie Description | Ce sont des oiseaux compacts de taille petite à moyenne (10 à 46 cm), au bec droit et long, en forme de poignard. Leurs pattes sont courtes, et ils portent un plumage aux couleurs vives. Répartition et habitat | Cosmopolites, ils fréquentent surtout
Read more

Can an atomic nucleus contain both particles and antiparticles? [duplicate]

Image
9 $begingroup$ This question already has an answer here: What happens if we put together a proton and an antineutron? 2 answers Is it theoretically possible to make a "deuterium" atom containing a proton and an antineutron in its nucleus? Would the strong nuclear force cause attraction between a proton and an antineutron? Would such a nucleus be stable, or would the proton somehow annihilate the antineutron when close enough? nuclear-physics antimatter share | cite | improve this question asked yesterday
Read more