solving two of nonlinear equations and two unknwons.



Hi

I have a problem about solving two of nonlinear equations. And I only
have two unknown.
The function is run by some complex program and I think it's very hard
to solve it analytically.
Thus, I cannot have anallytic Jacobian.
I try numerically. I use IMSL and fortran.

First I use NEQBF(purpose:Solves a system of nonlinear equations using
factored secant update with a finite-difference approximation to the
Jacobian.). Then I find there is no constraint (reasonable bounds)
for this.
Does anybody know if there is constrainted nonlinear equations solver
for imsl?

Then I switch to use BCLSF(purpose:
Solves a nonlinear least squares problem subject to bounds on the
variables using a modified
Levenberg-Marquardt algorithm and a finite-difference Jacobian.)
It is quite sensitive to guess vectors. It traps to a wrong solution.
Since I have two unkown for two nonlinear equations. I doubt if I
need to use this (least squares).

Is there a better choice for IMSL fortran?

Mike
.



Relevant Pages

  • Re: solving two of nonlinear equations and two unknwons.
    ... I have a problem about solving two of nonlinear equations. ... I cannot have anallytic Jacobian. ...  I use IMSL and fortran. ... What does "Number of good digits in the function." ...
    (sci.math.num-analysis)
  • Re: solving two of nonlinear equations and two unknwons.
    ... I have a problem about solving two of nonlinear equations. ... I use IMSL and fortran. ... A nonlinear system of dimension 2 and 2 unknowns has ... If you add a constraint, ...
    (sci.math.num-analysis)
  • Re: solving two of nonlinear equations and two unknwons.
    ... I have a problem about solving two of nonlinear equations. ... have two unknown. ... I use IMSL and fortran. ... If you add a constraint, ...
    (sci.math.num-analysis)
  • Re: solving two of nonlinear equations and two unknwons.
    ... I have a problem about solving two of nonlinear equations. ... I cannot have anallytic Jacobian. ... Solves a nonlinear least squares problem subject to bounds on the ... you are using codes which try to approximate derivatives (Jacobians) ...
    (sci.math.num-analysis)
  • Re: out of memory message
    ... 300 nonlinear equations equations in the 100 grid points for each ... Therefore the Jacobian that fsolves needs to handle is ... 30000x30000 in dimension in my problem. ... function and pass that options structure into FSOLVE, ...
    (comp.soft-sys.matlab)