Mercurial > hg > octave-max

changeset 12802:412882f498b4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
codesprint: Wrote 5 tests for quadl.m
author David Wells <drwells@vt.edu>
date Sat, 16 Jul 2011 12:52:47 -0500
parents 2837514cf178
children b7a6a3644f3b
files scripts/general/quadl.m
diffstat 1 files changed, 22 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/general/quadl.m
+++ b/scripts/general/quadl.m
@@ -194,3 +194,25 @@
     w = v;
   endif
 endfunction
+
+
+## basic functionality
+%!assert( quadl (@(x) sin (x), 0, pi, [], []), 2, -3e-16)
+
+## the values here are very high so it may be unavoidable that this fails
+%!assert ( quadl (@(x) sin (3*x).*cosh (x).*sinh (x),10,15),
+%!         2.588424538641647e+10, -9e-15)
+
+## extra parameters
+%!assert (quadl (@(x,a,b) sin (a + b*x), 0, 1, [], [], 2, 3),
+%!        cos(2)/3 - cos(5)/3, - 3e-16)
+
+## test different tolerances. This test currently fails for a very high
+## tolerances.
+%!assert ( quadl (@(x) sin (2 + 3*x).^2, 0, 10, 0.3, []),
+%!        (60 + sin(4) - sin(64))/12, -0.3)
+
+
+## for lower tolerances the test passes.
+%!assert ( quadl (@(x) sin (2 + 3*x).^2, 0, 10, 0.1, []),
+%!        (60 + sin(4) - sin(64))/12, -0.1)
\ No newline at end of file