new file mode 100644
--- /dev/null
+++ b/test/test_poly.m
@@ -0,0 +1,176 @@
+%% Automatically generated from DejaGNU files
+
+%% test/octave.test/poly/compan-1.m
+%!assert(all (all (compan ([1, 2, 3]) == [-2, -3; 1, 0])));
+
+%% test/octave.test/poly/compan-2.m
+%!assert(all (all (compan ([1; 2; 3]) == [-2, -3; 1, 0])));
+
+%% test/octave.test/poly/compan-3.m
+%!assert(isempty (compan (4)));
+
+%% test/octave.test/poly/compan-4.m
+%!assert(all (all (compan ([3, 2, 1]) == [-2/3, -1/3; 1, 0])));
+
+%% test/octave.test/poly/compan-5.m
+%!error compan ([1,2;3,4]);
+
+%% test/octave.test/poly/compan-6.m
+%!error compan ([]);
+
+%% test/octave.test/poly/conv-1.m
+%!assert(all (all (conv (ones (3, 1), ones (3, 1)) == [1, 2, 3, 2, 1])));
+
+%% test/octave.test/poly/conv-2.m
+%!assert(all (all (conv (ones (1, 3), ones (3, 1)) == [1, 2, 3, 2, 1])));
+
+%% test/octave.test/poly/conv-3.m
+%!assert(all (all (conv (3, [1, 2, 3]) == [3, 6, 9])));
+
+%% test/octave.test/poly/conv-4.m
+%!error conv ([1, 2; 3, 4], 3);
+
+%% test/octave.test/poly/conv-5.m
+%!assert(conv (2, 3),6);
+
+%% test/octave.test/poly/conv-6.m
+%!error conv (2, []);
+
+%% test/octave.test/poly/deconv-1.m
+%!test
+%! [b, r] = deconv ([3, 6, 9, 9], [1, 2, 3]);
+%! assert(all (all (b == [3, 0])) && all (all (r == [0, 0, 0, 9])));
+
+%% test/octave.test/poly/deconv-2.m
+%!test
+%! [b, r] = deconv ([3, 6], [1, 2, 3]);
+%! assert(b == 0 && all (all (r == [0, 3, 6])));
+
+%% test/octave.test/poly/deconv-3.m
+%!test
+%! [b, r] = deconv ([3, 6], [1; 2; 3]);
+%! assert(b == 0 && all (all (r == [0, 3, 6])));
+
+%% test/octave.test/poly/deconv-4.m
+%!error [b, r] = deconv ([3, 6], [1, 2; 3, 4]);;
+
+%% test/octave.test/poly/deconv-5.m
+%!error <number of rows must match> [b, r] = deconv ([3; 6], [1, 2, 3]);
+
+%% test/octave.test/poly/poly-1.m
+%!assert(all (all (poly ([1, 2, 3]) == [1, -6, 11, -6])));
+
+%% test/octave.test/poly/poly-2.m
+%!assert(all (all (abs (poly ([1, 2; 3, 4]) - [1, -5, -2]) < sqrt (eps))));
+
+%% test/octave.test/poly/poly-3.m
+%!error poly ([1, 2, 3; 4, 5, 6]);
+
+%% test/octave.test/poly/poly-4.m
+%!assert(poly ([]),1);
+
+%% test/octave.test/poly/polyderiv-1.m
+%!assert(all (all (polyderiv ([1, 2, 3]) == [2, 2])));
+
+%% test/octave.test/poly/polyderiv-2.m
+%!assert(polyderiv (13) == 0);
+
+%% test/octave.test/poly/polyderiv-3.m
+%!error polyderiv ([]);
+
+%% test/octave.test/poly/polyderiv-4.m
+%!error polyderiv ([1, 2; 3, 4]);
+
+%% test/octave.test/poly/polyfit-1.m
+%!test
+%! x = [-2, -1, 0, 1, 2];
+%! assert(all (all (abs (polyfit (x, x.^2+x+1, 2) - [1, 1, 1]) < sqrt (eps))));
+
+%% test/octave.test/poly/polyfit-2.m
+%!test
+%! x = [-2, -1, 0, 1, 2];
+%! assert(all (all (abs (polyfit (x, x.^2+x+1, 3) - [0, 1, 1, 1]) < sqrt (eps))));
+
+%% test/octave.test/poly/polyfit-3.m
+%!error polyfit ([1, 2; 3, 4], [1, 2; 3, 4], 4);
+
+%% test/octave.test/poly/polyfit-4.m
+%!test
+%! x = [-2, -1, 0, 1, 2];
+%! fail("polyfit (x, x.^2+x+1)");
+
+%% test/octave.test/poly/polyfit-5.m
+%!test
+%! x = [-2, -1, 0, 1, 2];
+%! fail("polyfit (x, x.^2+x+1, [])");
+
+%% test/octave.test/poly/polyinteg-1.m
+%!assert(all (all (polyinteg ([2, 2]) == [1, 2, 0])));
+
+%% test/octave.test/poly/polyinteg-2.m
+%!assert(isempty (polyinteg ([])));
+
+%% test/octave.test/poly/polyinteg-3.m
+%!assert(all (all (polyinteg (3) == [3, 0])));
+
+%% test/octave.test/poly/polyinteg-4.m
+%!error polyinteg ([1, 2; 3, 4]);
+
+%% test/octave.test/poly/polyreduce-1.m
+%!assert(all (all (polyreduce ([0, 0, 1, 2, 3]) == [1, 2, 3])));
+
+%% test/octave.test/poly/polyreduce-2.m
+%!assert(all (all (polyreduce ([1, 2, 3, 0, 0]) == [1, 2, 3, 0, 0])));
+
+%% test/octave.test/poly/polyreduce-3.m
+%!assert(all (all (polyreduce ([1, 0, 3]) == [1, 0, 3])));
+
+%% test/octave.test/poly/polyreduce-4.m
+%!assert(isempty (polyreduce ([])));
+
+%% test/octave.test/poly/polyreduce-5.m
+%!error polyreduce ([1, 2; 3, 4]);
+
+%% test/octave.test/poly/polyval-1.m
+%!assert(polyval ([1, 1, 1], 2) == 7);
+
+%% test/octave.test/poly/polyval-2.m
+%!assert(all (all (polyval ([1, 1, 1], [0; 1; 2]) == [1; 3; 7])));
+
+%% test/octave.test/poly/polyval-3.m
+%!assert(isempty (polyval ([1, 1, 1], [])));
+
+%% test/octave.test/poly/polyval-4.m
+%!assert(all (all (polyval ([1, 1, 1], [-1, 0; 1, 2]) == [1, 1; 3, 7])));
+
+%% test/octave.test/poly/polyval-5.m
+%!error polyval ([1, 2; 3, 4], [-1, 0; 1, 2]);
+
+%% test/octave.test/poly/polyval-6.m
+%!assert(isempty (polyval ([], [-1, 0; 1, 2])));
+
+%% test/octave.test/poly/polyvalm-1.m
+%!assert(isempty (polyvalm ([], [1, 2; 3, 4])));
+
+%% test/octave.test/poly/polyvalm-2.m
+%!error polyvalm ([1, 1, 1], [1, 2; 3, 4; 5, 6]);
+
+%% test/octave.test/poly/residue-1.m
+%!test
+%! b = [1, 1, 1];
+%! a = [1, -5, 8, -4];
+%! [r, p, k, e] = residue (b, a);
+%! assert((abs (r - [-2; 7; 3]) < 1e-6
+%! && abs (p - [2; 2; 1]) < 1e-7
+%! && isempty (k)
+%! && e == [1; 2; 1]));
+
+%% test/octave.test/poly/roots-1.m
+%!assert(all (all (abs (roots ([1, -6, 11, -6]) - [3; 2; 1]) < sqrt (eps))));
+
+%% test/octave.test/poly/roots-2.m
+%!assert(isempty (roots ([])));
+
+%% test/octave.test/poly/roots-3.m
+%!error roots ([1, 2; 3, 4]);
+