Mercurial > octave-nkf

comparison scripts/polynomial/residue.m @ 6998:cc049a392a97

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-10-10 16:51:51 by jwe]
author jwe
date Wed, 10 Oct 2007 16:53:23 +0000
parents 3fbf49bde6d4
children 4a682c7b2bd6
comparison
equal deleted inserted replaced
6997:6d0d8d621532 6998:cc049a392a97
300 300
301 %!test 301 %!test
302 %! b = [1, 1, 1]; 302 %! b = [1, 1, 1];
303 %! a = [1, -5, 8, -4]; 303 %! a = [1, -5, 8, -4];
304 %! [r, p, k, e] = residue (b, a); 304 %! [r, p, k, e] = residue (b, a);
305 %! assert ((abs (r - [-2; 7; 3]) < 1e-5 305 %! assert (abs (r - [-2; 7; 3]) < 1e-5
306 %! && abs (p - [2; 2; 1]) < 1e-7 306 %! && abs (p - [2; 2; 1]) < 1e-7
307 %! && isempty (k) 307 %! && isempty (k)
308 %! && e == [1; 2; 1])); 308 %! && e == [1; 2; 1]);
309 %! k = [1 0]; 309 %! k = [1 0];
310 %! [b, a] = residue (r, p, k); 310 %! [b, a] = residue (r, p, k);
311 %! assert ((abs (b - [1, -5, 9, -3, 1]) < 1e-12 311 %! assert ((abs (b - [1, -5, 9, -3, 1]) < 1e-12
312 %! && abs (a - [1, -5, 8, -4]) < 1e-12)); 312 %! && abs (a - [1, -5, 8, -4]) < 1e-12));
313 313
314 %!test 314 %!test
315 %! b = [1, 0, 1]; 315 %! b = [1, 0, 1];
316 %! a = [1, 0, 18, 0, 81]; 316 %! a = [1, 0, 18, 0, 81];
317 %! [r, p, k, e] = residue(b, a); 317 %! [r, p, k, e] = residue(b, a);
318 %! r1 = [-5i; 12; +5i; 12]/54; 318 %! r1 = [-5i; 12; +5i; 12]/54;
319 %! r2 = conj(r1); 319 %! r2 = conj(r1);
320 %! p1 = [+3i; +3i; -3i; -3i]; 320 %! p1 = [+3i; +3i; -3i; -3i];
321 %! p2 = conj(p1); 321 %! assert (abs (r - r1) < 1e-7 && abs (p - p1) < 1e-7
322 %! assert ((((abs (r - r1) < 1e-7) && (abs (p - p1) < 1e-7)) 322 %! && isempty (k)
323 %! || ((abs (r - r2) < 1e-7) && (abs (p - p2) < 1e-7))) 323 %! && e == [1; 2; 1; 2]);
324 %! && isempty (k)
325 %! && (e == [1; 2; 1; 2]));
326 %! [br, ar] = residue (r, p, k); 324 %! [br, ar] = residue (r, p, k);
327 %! assert ((abs (br - b) < 1e-12 325 %! assert ((abs (br - b) < 1e-12
328 %! && abs (ar - a) < 1e-12)); 326 %! && abs (ar - a) < 1e-12));