That is an error, as the code on octave-forge is not part of Octave. The
copyright strings were recently updated in octave-forge and I believe
roots is an octave core function and not from Octave forge, so send the
proposed patches to roots to address@hidden. As for octave-forge
The correct implementation of roots.m is
function z = roots(p)
if(any(isnan(p) | isinf(p)))
error('Input to ROOTS must not contain NaN or Inf.');
end
while(any(isinf(p./p(1))))
p=p(2:end);
end
p = vec(p);
n = numel(p)-1;
A = diag(ones(n-1,1),-1);
A(1,:) = -p(2:n+1)./p(1);
z = eig(A);
changes, those are easier to deal with as the acceptance process of
patches is easy.. Just send to the address@hidden list
and in generally you then commit the code yorself..
OK, I'll see what can be done for that.