Finding the solution of this long equation using Maple

Click For Summary

Discussion Overview

The discussion revolves around solving a complex equation using Maple, specifically for very small values of ##x## (##x \ll 1##). Participants explore commands and methods for finding solutions, including potential reformulations of the equation and numerical approximations.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant requests assistance with a lengthy equation in Maple, emphasizing the need for a command to find small values of ##x##.
  • Another participant reformats the equation using symbols for lengthy decimal numbers, suggesting this may improve readability and workability.
  • There is mention of a potential root for ##x \ll 1## lying between ##8.071385266 \times 10^{-19}## and ##8.071385267 \times 10^{-19}##, though the accuracy of this reformulation is questioned.
  • A suggestion is made to use a Taylor expansion in Maple to simplify the equation to a polynomial form, which may facilitate finding solutions.
  • Further exploration indicates the possibility of two roots for ##x \ll 1##, with one root near ##3.0381613366 \times 10^{-37}## and another near ##6.8635202356 \times 10^{-28}##, but caution is advised due to the delicate nature of the problem.
  • Some participants express uncertainty about the correctness of the reformulations and calculations, urging careful verification of each step.

Areas of Agreement / Disagreement

Participants express differing views on the accuracy of the reformulations and the existence of roots, indicating that the discussion remains unresolved with multiple competing perspectives on the approach to solving the equation.

Contextual Notes

There are limitations regarding the assumptions made about the equation's structure and the potential for errors in reformulation. The discussion highlights the complexity of the equation and the sensitivity of solutions to small changes in parameters.

jamesbrazil
Messages
1
Reaction score
0
TL;DR
Command to solve a strange equation in Maple for ##x \ll 1##
I need help solving an equation. I started using Maple, but had no success. Could someone explain to me which command to use? I need to find a very small value of ##x##, that is, ##x \ll 1##. The equation is

$$434972871000000000.0+{\frac {\sqrt {6} \left( { 1.488388992\times 10^
{-36}}\,\ln \left( 1/12\,{\frac {12\,\sqrt {6}{x}^{2}+{
1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {6\,{x}^{2}+{
1.488388992\times 10^{-36}}}x}{x}} \right) \sqrt {6}-12\,
\sqrt {6\,{x}^{2}+{ 1.488388992\times 10^{-36}}}x \right) }{
72\,{x}^{3}}}-{\frac { \left( 0.001704000000\,{x}^{2}+{
1.488388992\times 10^{-36}} \right) \left( { 1.488388992\times 10^{
-36}}\,\ln \left( 1/12\,{\frac { 0.003408000000\,\sqrt {6}{x}^{2
}+{ 1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}x
}{x}} \right) \sqrt {6}-12\,\sqrt { 0.0000004839360000\,{x}^
{2}+{ 4.227024737\times 10^{-40}}}x \right) \sqrt {6}}{72\,\sqrt
{ 6.000000001\,{x}^{2}+{ 5.240806311\times 10^{-33}}}\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}{
x}^{3}}}=0$$
 
Last edited by a moderator:
Physics news on Phys.org
jamesbrazil said:
TL;DR Summary: Command to solve a strange equation in Maple for ##x \ll 1##

I need help solving an equation. I started using Maple, but had no success. Could someone explain to me which command to use? I need to find a very small value of ##x##, that is, ##x \ll 1##. The equation is

$$434972871000000000.0+{\frac {\sqrt {6} \left( { 1.488388992\times 10^
{-36}}\,\ln \left( 1/12\,{\frac {12\,\sqrt {6}{x}^{2}+{
1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {6\,{x}^{2}+{
1.488388992\times 10^{-36}}}x}{x}} \right) \sqrt {6}-12\,
\sqrt {6\,{x}^{2}+{ 1.488388992\times 10^{-36}}}x \right) }{
72\,{x}^{3}}}-{\frac { \left( 0.001704000000\,{x}^{2}+{
1.488388992\times 10^{-36}} \right) \left( { 1.488388992\times 10^{
-36}}\,\ln \left( 1/12\,{\frac { 0.003408000000\,\sqrt {6}{x}^{2
}+{ 1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}x
}{x}} \right) \sqrt {6}-12\,\sqrt { 0.0000004839360000\,{x}^
{2}+{ 4.227024737\times 10^{-40}}}x \right) \sqrt {6}}{72\,\sqrt
{ 6.000000001\,{x}^{2}+{ 5.240806311\times 10^{-33}}}\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}{
x}^{3}}}=0$$
Can you repost this equation, using symbols for all of the lengthy decimal numbers, followed with a separate list of all the symbols and their respective decimal values? That would make everything much easier to read and work with!
 
If anyone might still be interested in this then I've tried to reformat it as requested by jamesbrazil

a=434972871000000000.0;
b=1.488388992*10^-36;
c=0.001704000000;
d=0.003408000000;
e=0.0000004839360000;
f=4.227024737*10^-40;
g=5.240806311*10^-33;
a+
(sqrt(6)*(b*ln(1/12*
12*sqrt(6)*x^2+b*sqrt(6)+
12*sqrt(6*x^2+b)*x)/x))*sqrt(6)-
12*sqrt(6*x^2+b)*x)/(72*x^3)-
((c*x^2+b)*
(b*ln(1/12*(d*sqrt(6)*x^2+b*sqrt(6)+
12*sqrt(e*x^2+f)*x)/x)*sqrt(6)-
12*sqrt(e*x^2+f)*x)*sqrt(6))/
(72*sqrt(6*x^2+g)*
sqrt(e*x^2+f)*x^3)

I believe there is a root for x<<1 that lies between
8.071385266*10^-19 and 8.071385267*10^-19

It is certainly possible that I've made a mistake when trying to reformat this.
And I suppose it might be possible that the OP made a mistake when trying to format this.

Please check this carefully to try to find any of my mistakes before depending on this.

When I look at a plot of this it seems that the expression is approximately 4.34973*10^17
for a range of modest positive x.
It also seems that it is approximately that same value for many small negative x values,
but it seems there may be a number of points which are indeterminant.
 
Last edited:
jamesbrazil said:
TL;DR Summary: Command to solve a strange equation in Maple for ##x \ll 1##

I need help solving an equation. I started using Maple, but had no success. Could someone explain to me which command to use? I need to find a very small value of ##x##, that is, ##x \ll 1##. The equation is

$$434972871000000000.0+{\frac {\sqrt {6} \left( { 1.488388992\times 10^
{-36}}\,\ln \left( 1/12\,{\frac {12\,\sqrt {6}{x}^{2}+{
1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {6\,{x}^{2}+{
1.488388992\times 10^{-36}}}x}{x}} \right) \sqrt {6}-12\,
\sqrt {6\,{x}^{2}+{ 1.488388992\times 10^{-36}}}x \right) }{
72\,{x}^{3}}}-{\frac { \left( 0.001704000000\,{x}^{2}+{
1.488388992\times 10^{-36}} \right) \left( { 1.488388992\times 10^{
-36}}\,\ln \left( 1/12\,{\frac { 0.003408000000\,\sqrt {6}{x}^{2
}+{ 1.488388992\times 10^{-36}}\,\sqrt {6}+12\,\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}x
}{x}} \right) \sqrt {6}-12\,\sqrt { 0.0000004839360000\,{x}^
{2}+{ 4.227024737\times 10^{-40}}}x \right) \sqrt {6}}{72\,\sqrt
{ 6.000000001\,{x}^{2}+{ 5.240806311\times 10^{-33}}}\sqrt {
0.0000004839360000\,{x}^{2}+{ 4.227024737\times 10^{-40}}}{
x}^{3}}}=0$$
If you know that x is small, I'd suggest a Taylor expansion first (use mtaylor in Maple), then solve it. That way it's just a polynomial equation.

-Dan
 
  • Like
Likes   Reactions: Greg Bernhardt
If anyone might still be interested in this then I've tried to reformat it as requested by jamesbrazil.

In experimenting with that I see that I dropped a couple of ( ) in my post above.

I apologize for that and I can't seem to edit that now to correct those mistakes.

This will hopefully correct my ( ) mistakes.

a=434972871000000000.0;
b=1.488388992*10^-36;
c=0.001704000000;
d=0.003408000000;
e=0.0000004839360000;
f=4.227024737*10^-40;
g=5.240806311*10^-33;
a+
(sqrt(6)*(b*ln(1/12*
(12*sqrt(6)*x^2+b*sqrt(6)+
12*sqrt(6*x^2+b)*x)/x))*sqrt(6)-
12*sqrt(6*x^2+b)*x)/(72*x^3)-
((c*x^2+b)*
(b*ln(1/12*(d*sqrt(6)*x^2+b*sqrt(6)+
12*sqrt(e*x^2+f)*x)/x)*sqrt(6)-
12*sqrt(e*x^2+f)*x)*sqrt(6))/
(72*sqrt(6*x^2+g)*
sqrt(e*x^2+f)*x^3)

Strangely enough, I now believe there may be two roots for x<<1

One may lie between 10^-37 and 10^-36 and may lie near 3.0381613366*10^-37

One may lie between 10^-28 and 10^-27 and may lie near 6.8635202356*10^-28

But with the size of the numbers involved this problem is very delicate and needs some care.

If anyone wants to try to check this then I would urge you to begin with checking that I have reformatted the original expression correctly and then begin carefully looking for the roots. I would greatly appreciate if someone would do that independently, I am always happy when someone finds and points out any of my mistakes.

Thank you and again I apologize for any and all of my errors.
 
Last edited:
Bill Simpson said:
If anyone might still be interested in this then I've tried to reformat it as requested by jamesbrazil.
Frankly, until @jamesbrazil returns to show that he's still interested (and to explain the problem background and show the derivation his equation) I'm don't think this is worth pursuing.
 
  • Like
Likes   Reactions: topsquark

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 0 ·
Replies
0
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
10
Views
3K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 11 ·
Replies
11
Views
8K
  • · Replies 8 ·
Replies
8
Views
2K