MathGroup Archive 2007

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

Search the Archive

FindRoot Precision Problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg73302] FindRoot Precision Problem
  • From: Scott Kasen <sdk5n at virginia.edu>
  • Date: Fri, 9 Feb 2007 23:45:25 -0500 (EST)

I'm running Mathematica 5.1 and have been struggling with this problem 
for the last few days...

I'm trying to use FindRoot to do a simultaneous solve on five equations 
involving five unknown variables.  Unfortunately for me, I continued 
getting the "Encountered a Singular Jacobian" error not far from my 
guess values.

In my attempts at debugging, I noticed the following peculiarity (I'm 
giving a very simplified example below):

If one of my equations is:
x1+x2 = = 0
then FindRoot is able to converge successfully.

However, if I change the equation to:
x1*100+x2*100 = = 0
I get the "Singular Jacobian" error.  It is this latter equation type 
that I need to be able to solve (i.e. the one with larger orders of 
magnitude.)

It is for this reason that I feel my problem is related to precision and 
accuracy.

I've tried playing around with the precision and accuracy, but to no 
avail.  A couple of things I've tried:  changing my PrecisionGoal, 
AccuracyGoal, and WorkingPrecision parameters, as well as globally 
increasing my precision using $MinPrecision.  My search of the math 
group archive pulled this up written by Daniel Lichtblau:
http://forums.wolfram.com/mathgroup/archive/1998/Mar/msg00336.html
This is hauntingly similar to my problem, and it is solved here using 
the GroebnerBasis function.  But, in my situation, my equations are 
nonlinear...I have some exponentials in my equations, and the 
GroebnerBasis function reminds me of that by producing the error, 
"Equations appear to be solved in an essentially non-algebraic way".

Any help would be immensely appreciated!  I've included the actual five 
equations in a Mathematica notebook file at the URL below, along with my 
FindRoot call and guesses.  I've tried to clean the equations up a 
little bit using Simplify before making this notebook file.

http://people.virginia.edu/~sdk5n/math_group.nb


-- 
______________________

Scott D. Kasen
Graduate Research Assistant

Department of Materials Science and Engineering
University of Virginia
P.O. Box 400745
116 Engineers Way
Charlottesville, VA 22904-4745

email:sdk5n at virginia.edu
voice: 434-982-5675
fax: 434-982-5660




  • Prev by Date: Re: Would someone confirm they also see this problem?
  • Next by Date: Re: Integrating SphericalHarmonicY
  • Previous by thread: Re: Array reference help please
  • Next by thread: Neural Networks add-on