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.










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.










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.










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






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

"Incorrect syntax near the keyword 'ON'. (on update cascade, on delete cascade,)

Alcedinidae

Origin of the phrase “under your belt”?