changeset 24914:d7293106945c stable

cosint fails for large x along negative real axis. cosint: add xtest.
author Colin Macdonald <cbm@m.fsf.org>
date Tue, 23 Jan 2018 23:55:12 -0800
parents 9da779d2f029
children c6cdacb1794e
files scripts/specfun/cosint.m
diffstat 1 files changed, 9 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/specfun/cosint.m	Tue Jan 23 23:14:50 2018 -0800
+++ b/scripts/specfun/cosint.m	Tue Jan 23 23:55:12 2018 -0800
@@ -226,3 +226,12 @@
 %! assert (A, B, -3*eps)
 %! B = cosint (single (x));
 %! assert (A, B, -3*eps ("single"))
+
+## fails along negative real axis
+%!xtest
+%! x = [-25 -100 -1000];
+%! yex = [-0.0068485971797025909189 + pi*1i
+%!        -0.0051488251426104921444 + pi*1i
+%!        0.000826315511090682282 + pi*1i];
+%! y = cosint (x);
+%! assert (y, yex, -5*eps)