[Octave-bug-tracker] [bug #33863] Two test fails in quadl.m (dev. source MinGW)

[Tatsuro MATSUOKA]

URL:
  <http://savannah.gnu.org/bugs/?33863>

                 Summary: Two test fails in quadl.m (dev. source MinGW)
                 Project: GNU Octave
            Submitted by: tmacchant
            Submitted on: 2011年07月26日 04時29分56秒
                Category: Interpreter
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Incorrect Result
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

Hello

I have encountered two test fails in quadl.m


  ***** assert ( quadl (@(x) sin (3*x).*cosh (x).*sinh (x),10,15),
         2.588424538641647e+10, -9e-15)
!!!!! test failed
assert (quadl (@(x) sin (3 * x) .* cosh (x) .* sinh (x), 10,
15),2.588424538641647e+10,-9e-15) expected
2.5884e+010
but got
2.5884e+010
maximum relative error 1.04636e-014 exceeds tolerance 9e-015



  ***** assert ( quadl (@(x) sin (2 + 3*x).^2, 0, 10, 0.3, []),
        (60 + sin(4) - sin(64))/12, -0.3)
!!!!! test failed
assert (quadl (@(x) sin (2 + 3 * x) .^ 2, 0, 10, 0.3, []),(60 + sin (4) - sin
(64)) / 12,-0.3) expected
 4.8603
but got
 8.1963
maximum relative error 0.686392 exceeds tolerance 0.3


For the first fail,
> quadl (@(x) sin (3 * x) .* cosh (x) .* sinh (x), 10, 15) -
2.588424538641647e+10
ans = -2.70843505859375e-004

Error seems large.

For second fail, error is apparently large.
I tried


octave:8> quad (@(x) sin (2 + 3 * x) .^ 2, 0, 10, 0.3, [])
warning: implicit conversion from null_matrix to real vector
ans =  4.86026428887461


quadl does not work correctly in some cases.

Regards

Tatsuro




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?33863>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/