FLOAD C:\WIN32FOR\JVN_PROGS\SIMP_INT.F

: f1   fdup  f0=  if  exit  then  fln  fcos  ;  ok
use( f1  0e0 1e0  1e-3 )integral cr f. 33  function calls
.50036707479007  ok

use( f1  0e0 1e0  1e-5 )integral cr f. 89  function calls
.49998854561883  ok

use( f1  0e0 1e0  1e-6 )integral cr f. 161  function calls
.50000066355684  ok

: f2   fdup  f^2  fnegate  1e0  f+  fsqrt  ;  ok
use( f2  -1e0 1e0  1e-4 )integral cr f.
Error: )integral Too many subdivisions!

use( f2  0e0 1e0  1e-4 )integral f2* cr f. 73  function calls
1.5707947345751  ok

use( f2  0e0 1e0  1e-6 )integral f2* cr f.
Error: )integral Too many subdivisions!

use( f2  0e0 1e0  1e-5 )integral f2* cr f.
Error: )integral Too many subdivisions!

: f3   f2  f^3  ;  ok
use( f3  0e0 1e0  1e-5 )integral f2* cr f. 41  function calls
1.1780972890983  ok

use( f3  -1e0 1e0  1e-5 )integral  cr f. 81  function calls
1.1780981124834  ok

use( f3  -1e0 1e0  1e-6 )integral  cr f. 125  function calls
1.1780973967987  ok

use( f3  -1e0 1e0  1e-7 )integral  cr f. 221  function calls
1.1780972500625  ok