# HG changeset patch
# User jwe
# Date 1196109085 0
# Node ID fdd7cd70dc14abe8af9762cc7facd02afac652d4
# Parent  1124c84fc63d848490024894feb10a828cb11024
[project @ 2007-11-26 20:31:25 by jwe]

diff --git a/scripts/polynomial/residue.m b/scripts/polynomial/residue.m
--- a/scripts/polynomial/residue.m
+++ b/scripts/polynomial/residue.m
@@ -378,4 +378,15 @@
 %!   && abs (kr - k) < 1e-12
 %!   && abs (er - e(n)) < 1e-12));
 
-
+%!test
+%! b = [1];
+%! a = [1, 10, 25];
+%! [r, p, k, e] = residue(b, a);
+%! r1 = [0; 1];
+%! p1 = [-5; -5];
+%! assert (abs (r - r1) < 1e-7 && abs (p - p1) < 1e-7
+%!   && isempty (k)
+%!   && e == [1; 2]);
+%! [br, ar] = residue (r, p, k);
+%! assert ((abs (br - b) < 1e-12
+%!   && abs (ar - a) < 1e-12));