Re: ODE system solving

[Geraint Paul Bevan]

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

address@hidden wrote:
| i've tryed to write an m-file containing all the instructions for
solving my equations. this is the file:
|
| # dovrebbe calcolare la T(z) e w(z) per una bava di GR, processo Neumag
| y0=[57 0.186 520];
| z=linspace(0,1000,2000)';
| y=lsode(SysO,y0,z);
| plot(z,y(:,3)-273)
| # inizio della descrizione del sistema
| function SysO=f(y,z)
...
| endfunction
|
|
| it doesn't work!!! octave crush!!!
| what can i do?

Do you have all of that in a single file? SysO.m? If so, you should
separate the function definition from the script:

run.m:
y0=[57 0.186 520];
z=linspace(0,1000,2000)';
y=lsode(SysO,y0,z);
plot(z,y(:,3)-273)

SysO.m:
function ...
endfunction


- --
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

iEYEARECAAYFAkHjuIYACgkQcXV3N50QmNM+cACeK/gnWYIWXWJhoTLbs4OA9CFb
t4MAnRDBs8L7dzasWKJOiSGMhHa3aQ1m
=hzsj
-----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
-------------------------------------------------------------