Re: definite numerical integration

[roberto]

On 9/21/05, David Bateman <address@hidden> wrote:
> roberto wrote:
> >3. doing this, evaluating the integral by varying r* in [0,R] (e.g. R = 110)
> >i'll have a function F(r*,a) where a should be obtained by imposing
> >some boundary condition;
> >
> >i looked for some help into Octave like
> >http://octave.sourceforge.net/index/Q.html
> >
> >but in libraries like:
> >[v, ier, nfun, err] = quad (f, a, b, tol, sing)
> >
> >the function to be integrated are to be inserted by their explicit 
> >formulation
> >which i don't have now, as already stated above;
> >
> >
> >
> I don't quite see the problem as in 1) you saw you have a function to
> numerically calculate the values in a range [0,r_max] and then in 3) you
> state you can use a explicit formulation
> 
sorry, but actually i do not have an explicit formulation of f(r) like 
function y = fun (x)
# Nasty integrable singularity to be unkind to quad and show octave's
   better here :-)
     y  = 1/x;
endfunction

but instead i have only values of f(r) in a finite set of points r,
maybe this is useful to solve the problem



-- 
roberto



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