--- a/src/DLD-FUNCTIONS/quad.cc	Thu Mar 06 01:56:55 2008 -0500
+++ b/src/DLD-FUNCTIONS/quad.cc	Thu Mar 06 02:27:55 2008 -0500
@@ -320,6 +320,33 @@
 }
 
 /*
+
+%!function y = f (x) 
+%! y = x + 1;
+%!test
+%! [v, ier, nfun, err] = quad ("f", 0, 5);
+%! assert(ier == 0 && abs (v - 17.5) < sqrt (eps) && nfun > 0 && 
+%!        err < sqrt (eps))
+
+%!function y = f (x)
+%!  y = x .* sin (1 ./ x) .* sqrt (abs (1 - x));
+%!test
+%!  [v, ier, nfun, err] = quad ("f", 0.001, 3);
+%! assert((ier == 0 || ier == 1) && abs (v - 1.98194120273598) < sqrt (eps) && nfun > 0);
+
+%!error <Invalid call to quad.*> quad ();
+
+%!error <Invalid call to quad.*> quad ("f", 1, 2, 3, 4, 5);
+
+%!test
+%! quad_options ("absolute tolerance", eps);
+%! assert(quad_options ("absolute tolerance") == eps);
+
+%!error <Invalid call to quad_options.*> quad_options (1, 2, 3);
+
+*/
+
+/*
 ;;; Local Variables: ***
 ;;; mode: C++ ***
 ;;; End: ***