From bug-octave-request at bevo dot che dot wisc dot edu  Fri Oct 10 14:48:29 1997
Subject: Re: archive retrieval: welcome.txt
From: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
To: "Raymond E. Rogers" <rrogers at voyager dot net>
Cc: "bug-octave at bevo dot che dot wisc dot edu" <bug-octave@bevo.che.wisc.edu>
Date: Fri, 10 Oct 1997 14:46:36 -0500

On 10-Oct-1997, Raymond E. Rogers <rrogers at voyager dot net> wrote:

| The following is a cut&paste of a gammai failure.  Could I get others to
| check if it is a octave failure or some failure in my Linux machine? 
| Notice that it is corrected with the "whitspace" command.

| *************session***********
| /root# octave --braindead
| Octave, version 2.0.9 (i586-pc-linux-gnu).
| Copyright (C) 1996, 1997 John W. Eaton.
| This is free software with ABSOLUTELY NO WARRANTY.
| For details, type `warranty'.
| 
| >> gammai(7,8)
| parse error near line 78 of file
| /usr/local/share/octave/2.0.9/m/specfun/gammai.m
| 
| >>>     u   = [zeros (1, len); ones (1, len)];
|                          ^
| 
| error: `gammai' undefined near line 1 column 1
| error: evaluating index expression near line 1, column 1

Please try the following patch.

Thanks,

jwe


Fri Oct 10 11:18:10 1997  John W. Eaton  <jwe at bevo dot che dot wisc dot edu>

	* specfun/gammai.m: Avoid problems with whitespace when
	constructing matrices.


Index: specfun/gammai.m
===================================================================
RCS file: /home/jwe/src/master/octave/scripts/specfun/gammai.m,v
retrieving revision 1.14
diff -c -r1.14 gammai.m
*** gammai.m	1996/11/21 01:43:06	1.14
--- gammai.m	1997/10/10 19:43:01
***************
*** 75,82 ****
    k = find ((x >= a + 1) & (x < Inf) & (a > 0));
    if any (k)
      len = length (k);
!     u   = [zeros (1, len); ones (1, len)];
!     v   = [ones (1, len); x(k)];
      c_old = 0;
      c_new = v(1, :) ./ v(2, :);
      n   = 1;
--- 75,84 ----
    k = find ((x >= a + 1) & (x < Inf) & (a > 0));
    if any (k)
      len = length (k);
!     t0  = zeros (1, len);
!     t1  = ones (1, len);
!     u   = [t0; t1];
!     v   = [t1; x(k)];
      c_old = 0;
      c_new = v(1, :) ./ v(2, :);
      n   = 1;