Re: Solving a trivial ODE with lsode?

[Geraint Paul Bevan]

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

Dirk Koschützki wrote on 02/12/04 12:09:
| Hello,
|
| I'm trying to solve a trivial ODE with octave (2.1.35) . This example is
| used just to get a feeling for octave.
|
| My problem is: I did not get the correct solution.
|

The correct solution of \int_0^{20} 2x dt" is not x^2, so your target
function 'quadrat' is not what the result of using lsode should show.

Your initial state is 0. At t=0, the derivative is 2*0=0, hence the next
~ value of the state is also 0.

Integrating a system that has no external inputs and is in a steady
state (dx/dt=x=0) will not give you a very interesting simulation.

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

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.2 (GNU/Linux)

iEYEARECAAYFAkArn7kACgkQcXV3N50QmNNwJACZAQPrlrdWjJcqIxAentiuPWBN
H/gAoIIWaD5Vs6BclLlv+Yfn7puopw/1
=fEnE
-----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
-------------------------------------------------------------