Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Are some equations unsolvable?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg128079] Re: Are some equations unsolvable?
  • From: "Alexei Boulbitch" <boulbitch at gmx.de>
  • Date: Fri, 14 Sep 2012 00:22:55 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net

Hi,
This is what I have:

f = (1/(sigma*Sqrt[2 Pi]))*Sqrt[(Pi/(-1/(2*sigma^2)))]*
  Exp[((m/sigma^2)^2)/4*(-1/2*sigma^2) + (m^2/(2*sigma^2))]

Solve[f == 216, sigma]

And I get this message: "This system cannot be solved with the methods available to Solve" 

Is it because there is no way to isolate sigma? Or am I doing something wrong?

Thanks

Some are unsolvable indeed. And not only due to limitations of Mathematica, but just due to their internal nature.
In addition come limitations are imposed also by Mathematica. If you look into Menu/Help/Solve/More Information you find the following:

Solve deals primarily with linear and polynomial equations. 
The option InverseFunctions specifies whether Solve should use inverse functions to try and find solutions to more general equations. The default is InverseFunctions->Automatic. In this case, Solve can use inverse functions, but prints a warning message. See notes on InverseFunctions. 

That is, some non-polynomial equation it can solve, but not all of them.
It is a good idea to help Mathematica understand, of what do you need to achieve. In your case, for example, you may have specified, if the variables are real or complex. However, 
This is your left hand part:

In[19]:= Clear[f];
f = (1/(sigma*Sqrt[2 Pi]))*Sqrt[(Pi/(-1/(2*sigma^2)))]*
  Exp[((m/sigma^2)^2)/4*(-1/2*sigma^2) + (m^2/(2*sigma^2))]

ff = Simplify[f /. sigma -> x*m, {m > 0}]

Out[20]= (\[ExponentialE]^((3 m^2)/(8 sigma^2)) Sqrt[-sigma^2])/sigma

Out[21]= (\[ExponentialE]^(3/(8 x^2)) Sqrt[-x^2])/x


Now it is clear that if solution exists, it is imaginary:

In[17]:= FindRoot[(\[ExponentialE]^(3/(8 x^2)) Sqrt[-x^2])/x, {x, 
  0.3*\[ImaginaryI]}]

Out[17]= {x -> 0.+ 0.0592883 \[ImaginaryI]}

Have fun, Alexei
 

-- 
Alexei BOULBITCH,
Dr. , habil. 
IEE S.A.
Department for Material Development
ZAE Weiergewan
11, rue Edmond Reuter
L-5326 Contern, Luxembourg

Tel.          ++352-2454-2566
Fax.:        ++352 424737201
mobile:       +49 (0) 151 524 066 44  
E-mail:      alexei.boulbitch at iee.lu 
                boulbitch at gmx.de





  • Prev by Date: Re: Kernel exits
  • Next by Date: Re: Eigenvalue and eigenvectors of a 10x10 matrix
  • Previous by thread: Re: Are some equations unsolvable?
  • Next by thread: Eigenvalue and eigenvectors of a 10x10 matrix