dl and fortran

[Stef Pillaert BK]

Hello again,

I would like to use the LAPACK-function DSYSVX to solve symmetric (but not
positive definite) matrices.

Right now, I'm saving my matrices to a file, I call a fortran
program, where I call DSYSVX, I save
the results to files, and read the results in octave.

I would prefer to do all this with dynamic linking, but since I'm not a
c++-specialist (neither fortran...), I have a few little problems:

* dynamic linking works (well, not in 2.0.7, but I should wait for 2.0.8,
am I well informed?). But I have no idea how to call my fortran-program 
from within a file like "oregonator.cc" (see examples). 

* How about precision? Does octave use "double precision"-equivalent (I'm
talking Fortran now...) And how do I pass variables\values?

Is it possible to give me a small example how to do all this?

Thanks,

address@hidden