Re: Using quad() for multidimensional integration

[Ted Harding]

( Re Message From: Vinayak Dutt )
>  I have source code and had a look at it. I am now coming down to the view
> that the problem is due to global variables used in Octave. It seems that
> Octave does not support recursion, as Mike Rilee pointed to me earlier in a
> private e-mail. I myself checked that with a factorial function and it bombed.
> On the other hand, MATLAB does support recusion. But it still fails
> in quad(). So I don't know whats the problem with MATLAB code. But as far
> as I can see Octave's problem stem from no recursion support.

Octave does support recursion. Look at this:--

    octave:78> function  y = fact(x)
    > if x<1, y=gamma(x+1);
    > else y = x.*fact(x-1);
    > end
    > end

    octave:79> fact(5)
    ans = 120
    octave:80> fact(6)
    ans = 720
    octave:81> fact(7)
    ans = 5040
    octave:82> fact(20)
    ans =  2.4329e+18
    octave:83> fact(30)
    ans =  2.6525e+32

Looks OK to me!

Ted.                                    (address@hidden)