[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: ODE system solving
From: |
Geraint Paul Bevan |
Subject: |
Re: ODE system solving |
Date: |
Tue, 11 Jan 2005 14:43:38 +0000 |
User-agent: |
Mozilla Thunderbird 0.7.1 (X11/20040715) |
-----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
-------------------------------------------------------------
- ODE system solving, address@hidden, 2005/01/07
- Re:Re: ODE system solving, address@hidden, 2005/01/07
- Re: ODE system solving, address@hidden, 2005/01/10
- Re: ODE system solving, address@hidden, 2005/01/11
- Re: ODE system solving, address@hidden, 2005/01/11
- Re: ODE system solving,
Geraint Paul Bevan <=
- Re: ODE system solving, address@hidden, 2005/01/11