From bug-octave-request at che dot utexas dot edu  Sat Oct 30 16:25:29 1993
Subject: Crash and burn report
From: Tony Mullins <jamull at che dot utexas dot edu>
To: bug-octave
Date: Sat, 30 Oct 93 16:25:20 EDT


The following input causes Octave to crash on my DEC 5000/200.  I know
that there are bugs in this file, and I have fixed them.  However, I
thought you might want samples of input that cause Octave to crash.

Tony Mullins

******************** Cut Here ********************
alpha = 3;
nc = 5;

global k1=1;
global De=1;

[u,A,B] = colloc (nc,"right");

# First we transform the roots to r^2.

r = sqrt(u);

[m,n] = size(u);

for j=1:n
  a(j) = alpha-1;
endfor

# Here we are going to transform the 2nd-derivative weights to
# even polynomials.

global Br2;
Br2 = diag((2*u).^2)*B + 6*A;

function rhs = f(C)
  rhs = Br2*C + (k1/De)*C;
endfunction

C1 = ones(n);

[C,info] = fsolve("f",C1)