[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: fsolve features
From: |
Francesco Potorti` |
Subject: |
Re: fsolve features |
Date: |
Mon, 5 Jul 1999 12:35:12 +0200 (CEST) |
(2) a call to fsolve like : [x,info]=fsolve('f',x0) is adequate for
function (vector) Y = f ( vector X )
>>> where Y and X have same length <<<<
It's supposed to handle equally sized problems only!
It is easy to work around this. Maybe you are suggesting that a workaround
is nt what you would like to have?
[x,RC]= fsolve(f='x^2 +1',randn(1)); #>V2.0.14
#NOTE: prior to V2.0.14 , you have to define f explicitely:
#function y=f(x) y=x^2+1; end
#[x,RC]=fsolve('f',randn(1))
!!! Where is this documented?
Because of all these pitfalls of 'fsolve', I decided to write my my
personal solver 'cfsolve' , returning correct complex roots for this (or
similar) problem
Where is your cfsolve?
--
Francesco Potortì (researcher) Voice: +39-050-593 203 (op. 211)
Computer Networks Group Fax: +39-050-904052
CNUCE-CNR, Via Santa Maria 36 Email: address@hidden
56126 Pisa - Italy Web: http://fly.cnuce.cnr.it/
---------------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL. To ensure
that development continues, see www.che.wisc.edu/octave/giftform.html
Instructions for unsubscribing: www.che.wisc.edu/octave/archive.html
---------------------------------------------------------------------