help-octave
[Top][All Lists]

## Re: Using quad() for multidimensional integration

 From: Ted Harding Subject: Re: Using quad() for multidimensional integration Date: Thu, 17 Aug 1995 23:21:11 +0200 (BST)

```( 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!