Numerical integration/Adaptive Simpson's method: Difference between revisions

Content added Content deleted

Inline

@@ Line 566: / Line 566: @@
 ruleval1 picture 9(5)V9(20).
 delta picture 9(5)V9(20).
+abs-delta picture 9(5)V9(20).
 tol0 picture 9(5)V9(20).
 tol1 picture 9(5)V9(20).
@@ Line 590: / Line 591: @@
                compute delta = ruleval0 + ruleval1 - ruleval
                if delta < 0 then
-                   compute delta = -delta
+                   compute abs-delta = -delta
+               else
+                   compute abs-delta = delta
                end-if
                compute tol0 = tol * (x4 - x0)
                compute tol1 = tol * (x2 - x0)
-               if (tol1 = tol0) or (delta <= tol0) then
+               if (tol1 = tol0) or (abs-delta <= tol0) then
                    compute delta15 = delta / 15
                    compute stepval = ruleval0 + ruleval1 + delta15