Mercurial > hg > octave-nkf

comparison scripts/polynomial/residue.m @ 7188:fdd7cd70dc14

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-11-26 20:31:25 by jwe]
author jwe
date Mon, 26 Nov 2007 20:31:25 +0000
parents c0be321eb472
children 4f403ccc4e30
comparison
equal deleted inserted replaced
7187:1124c84fc63d 7188:fdd7cd70dc14
376 %! assert ((abs (rr - r(n)) < 1e-5 376 %! assert ((abs (rr - r(n)) < 1e-5
377 %! && abs (pr - p(n)) < 1e-7 377 %! && abs (pr - p(n)) < 1e-7
378 %! && abs (kr - k) < 1e-12 378 %! && abs (kr - k) < 1e-12
379 %! && abs (er - e(n)) < 1e-12)); 379 %! && abs (er - e(n)) < 1e-12));
380 380
381 381 %!test
382 %! b = [1];
383 %! a = [1, 10, 25];
384 %! [r, p, k, e] = residue(b, a);
385 %! r1 = [0; 1];
386 %! p1 = [-5; -5];
387 %! assert (abs (r - r1) < 1e-7 && abs (p - p1) < 1e-7
388 %! && isempty (k)
389 %! && e == [1; 2]);
390 %! [br, ar] = residue (r, p, k);
391 %! assert ((abs (br - b) < 1e-12
392 %! && abs (ar - a) < 1e-12));