Re: ODE system solving

[Geraint Paul Bevan]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

address@hidden wrote:
| i found this for solving a BVP problem:
|
| you have a condition like this:
| equations
| y(0)=2
| y(L)=4
|
| Using ODE function, you can't do this, but you have to write something
like this:
|
| equations
| y(0)=2
| Y'(0)=alfa
|
| then you write this new equation: F(alfa)=y(L)-4
|
| and you need to find alfa for wich F(alfa)=0 (so y(L)=4)
|
| my question is: how to solve this equation? can i use something like
the MATLAB's "fzero" function? (i think it's ok in matlab, but i don't
konw really)

Have a look at the Differential-Algebraic Equations and Nonlinear
equations sections in the manual - you should be able to use daspk to
solve the equations for a given alfa. Perhaps you could then wrap the
entire integration loop within another function that you pass to fsolve.

- --
Geraint Bevan
http://www.mech.gla.ac.uk/~gbevan

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.5 (GNU/Linux)
Comment: Using GnuPG with Thunderbird - http://enigmail.mozdev.org

iEYEARECAAYFAkHj5hoACgkQcXV3N50QmNP88gCgiJhZgIkKfGyMDtcm2j09ym4L
5FMAoIG5Zxk52sRyuKlOa27YiZB+z1AI
=Skwj
-----END PGP SIGNATURE-----



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------