Mercurial > hg > octave-nkf
diff scripts/polynomial/residue.m @ 6998:cc049a392a97
[project @ 2007-10-10 16:51:51 by jwe]
author | jwe |
---|---|
date | Wed, 10 Oct 2007 16:53:23 +0000 |
parents | 3fbf49bde6d4 |
children | 4a682c7b2bd6 |
line wrap: on
line diff
--- a/scripts/polynomial/residue.m +++ b/scripts/polynomial/residue.m @@ -302,14 +302,14 @@ %! b = [1, 1, 1]; %! a = [1, -5, 8, -4]; %! [r, p, k, e] = residue (b, a); -%! assert ((abs (r - [-2; 7; 3]) < 1e-5 -%! && abs (p - [2; 2; 1]) < 1e-7 -%! && isempty (k) -%! && e == [1; 2; 1])); +%! assert (abs (r - [-2; 7; 3]) < 1e-5 +%! && abs (p - [2; 2; 1]) < 1e-7 +%! && isempty (k) +%! && e == [1; 2; 1]); %! k = [1 0]; %! [b, a] = residue (r, p, k); %! assert ((abs (b - [1, -5, 9, -3, 1]) < 1e-12 -%! && abs (a - [1, -5, 8, -4]) < 1e-12)); +%! && abs (a - [1, -5, 8, -4]) < 1e-12)); %!test %! b = [1, 0, 1]; @@ -318,11 +318,9 @@ %! r1 = [-5i; 12; +5i; 12]/54; %! r2 = conj(r1); %! p1 = [+3i; +3i; -3i; -3i]; -%! p2 = conj(p1); -%! assert ((((abs (r - r1) < 1e-7) && (abs (p - p1) < 1e-7)) -%! || ((abs (r - r2) < 1e-7) && (abs (p - p2) < 1e-7))) -%! && isempty (k) -%! && (e == [1; 2; 1; 2])); +%! assert (abs (r - r1) < 1e-7 && abs (p - p1) < 1e-7 +%! && isempty (k) +%! && e == [1; 2; 1; 2]); %! [br, ar] = residue (r, p, k); %! assert ((abs (br - b) < 1e-12 -%! && abs (ar - a) < 1e-12)); +%! && abs (ar - a) < 1e-12));