MathGroup Archive 2010

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

Search the Archive

Re: Find the solution of a system of two nonlinear, complicated equations

  • To: mathgroup at smc.vnet.net
  • Subject: [mg108824] Re: Find the solution of a system of two nonlinear, complicated equations
  • From: Alois Steindl <Alois.Steindl at tuwien.ac.at>
  • Date: Fri, 2 Apr 2010 05:19:47 -0500 (EST)
  • References: <hp1ufs$6d4$1@smc.vnet.net>

Christian Schneider schrieb:
> Dear members of the mathgroup email group,
> 
> I have serious difficulties in solving the following problem in Mathematica
> v.7 and would very grateful if anyone could give me a helping hand with
> that.
> 
> I want to solve a system of two nonlinear equations. Unfortunately the
> equations are rather bulky. I've tried to use FindRoot for this problem, but
> I get three error messages:
> 
> General::ovfl : Overflow occurred in computation
> 
> General::ovfl : Overflow occurred in computation
> 
> FindRoot::nlnum : The function value {Overflow[], Overflow[]} is not a list
> of numbers with dimensions {2} at {sigmabare, dmu}={0.2, 4}.
> 
> I think mathematica is telling me that the iterations or the starting points
> need to be changed. I tried to do that, but always end up with that error
> messages.
> 
> I have a set of two equations with two variables and two parameters each,
> and I want to solve the two variables of the set.
> 
> The variables of the two equations are: sigmabare, dmu
> 
> The two parameters per equation are
> 
> for equation1: epsilon1, c01 
> 
> for equation2: epsilon2, c02
> 
>  
> 
> The complete syntax as plain text is:
> 
>  
> 
> FindRoot[{c01*6.022*10^(23) == 
> 
>        Abs[sigmabare/(2*(437*10^(-12))*3*(1.602*10^(-19)))]*
> 
>         Exp[(-(4.11641*10^(-21))*(1.65*(1/(4*(4.11641*10^(-21))*
> 
>                     epsilon1*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^\
> 
> 3*sigmabare/3.141592654)^(1/2)) - 
> 
>               2.61*(1/(4*(4.11641*10^(-21))*
> 
>                     epsilon1*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^\
> 
> 3*sigmabare/3.141592654)^(1/2))^(1/4) + 
> 
>               0.26*Log[
> 
>                 1/(4*(4.11641*10^(-21))*
> 
>                     epsilon1*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^\
> 
> 3*sigmabare/3.141592654)^(1/2)] + 1.95))/((4.11641*10^(-21)))]*
> 
>         Exp[dmu/((4.11641*10^(-21)))]/(6.022*10^(23)) /. 
> 
>       c01 -> {0.5*10^(-5)} /. epsilon1 -> {54} , 
> 
>     c02*6.022*10^(23) == 
> 
>      Abs[sigmabare/(2*(437*10^(-12))*3*(1.602*10^(-19)))]*
> 
>       Exp[(-(4.11641*10^(-21))*(1.65*(1/(4*(4.11641*10^(-21))*
> 
>                   epsilon2*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^3*
> 
>                   sigmabare/3.141592654)^(1/2)) - 
> 
>             2.61*(1/(4*(4.11641*10^(-21))*
> 
>                    epsilon2*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^3*
> 
>                    sigmabare/3.141592654)^(1/2))^(1/4) + 0.26*Log[
> 
>               1/(4*(4.11641*10^(-21))*
> 
>                   epsilon2*(8.854*10^(-12)))*((1.602*10^(-19))^3*3^3*
> 
>                   sigmabare/3.141592654)^(1/2)] + 
> 
>             1.95))/((4.11641*10^(-21)))]*
> 
>       Exp[dmu/((4.11641*10^(-21)))]/(6.022*10^(23))} /. 
> 
>    c02 -> {1*10^(-5)} /. 
> 
>   epsilon2 -> {78}, {{sigmabare, 0.2}, {dmu, 4}}]
> 
>  
> 
> Sorry for this mess of equations.
> 
>  
> 
> Has somebody an idea where I am going wrong and how to solve this?
> 
>  
> 
> Thanks a lot in advance for your help,
> 
>  
> 
> Chris
> 
Hello,
this looks really like a mess of equations!
I do not want to play around with your equations, instead I want to give 
you some general advices:
First you use several constants (like some numeric appoximation of pi 
and Avogadro's constant) throughout your equations without assigning 
them to symbols.
Second you use very large and small quantities simultaneously. You 
surely are aware of the limits of floating point calculations? I would 
advice to rescale your system before trying to solve it. You should use 
non-dimensional quantities and scale the system, such that all 
expressions in the equations are of order 1; also the expected solution 
should be of that order.
Alois


  • Prev by Date: Re: 0/1 knapsack-like minimalization problem and file
  • Next by Date: Re: Find the solution of a system of two nonlinear, complicated
  • Previous by thread: Re: Find the solution of a system of two nonlinear, complicated equations
  • Next by thread: Re: Adding data to excel sheet (v7.0)