diff scripts/general/quadv.m @ 8507:cadc73247d65

style fixes
author John W. Eaton <jwe@octave.org>
date Tue, 13 Jan 2009 14:08:36 -0500
parents df9519e9990c
children eb63fbe60fab
line wrap: on
line diff
--- a/scripts/general/quadv.m
+++ b/scripts/general/quadv.m
@@ -42,7 +42,7 @@
 ## @seealso{triplequad, dblquad, quad, quadl, quadgk, trapz}
 ## @end deftypefn
 
-function [Q, fcnt] = quadv (f, a, b, tol, trace, varargin)
+function [q, fcnt] = quadv (f, a, b, tol, trace, varargin)
   if (nargin < 3)
     print_usage ();
   endif
@@ -81,49 +81,49 @@
   endif
 
   h = (b - a) / 2;
-  Q = (b - a) / 6 * (fa + 4 * fc + fb);
+  q = (b - a) / 6 * (fa + 4 * fc + fb);
  
-  [Q, fcnt, hmin] = simpsonstp (f, a, b, c, fa, fb, fc, Q, fcnt, abs (b - a), 
+  [q, fcnt, hmin] = simpsonstp (f, a, b, c, fa, fb, fc, q, fcnt, abs (b - a), 
 				tol, trace, varargin{:});
 
   if (fcnt > 10000)
-    warning("Maximum iteration count reached");
-  elseif (isnan(Q) || isinf (Q))
-    warning ("Infinite or NaN function evaluations were returned");
+    warning ("maximum iteration count reached");
+  elseif (isnan (q) || isinf (q))
+    warning ("infinite or NaN function evaluations were returned");
   elseif (hmin < (b - a) * myeps)
-    warning ("Minimum step size reached. Possibly singular integral");
+    warning ("minimum step size reached -- possibly singular integral");
   endif
 endfunction
 
-function [Q, fcnt, hmin] = simpsonstp (f, a, b, c, fa, fb, fc, Q0, 
+function [q, fcnt, hmin] = simpsonstp (f, a, b, c, fa, fb, fc, q0, 
 				       fcnt, hmin, tol, trace, varargin)
   if (fcnt > 10000)
-    Q = Q0;
+    q = q0;
   else
     d = (a + c) / 2;
     e = (c + b) / 2;
     fd = feval (f, d, varargin{:});
     fe = feval (f, e, varargin{:});
     fcnt += 2;
-    Q1 = (c - a) / 6 * (fa + 4 * fd + fc);
-    Q2 = (b - c) / 6 * (fc + 4 * fe + fb);
-    Q = Q1 + Q2;
+    q1 = (c - a) / 6 * (fa + 4 * fd + fc);
+    q2 = (b - c) / 6 * (fc + 4 * fe + fb);
+    q = q1 + q2;
 
     if (abs(a -  c) < hmin)
       hmin = abs (a - c);
     endif
 
     if (trace)
-      disp ([fcnt, a, b-a, Q]);
+      disp ([fcnt, a, b-a, q]);
     endif
 
-    ## Force at least one adpative step
-    if (fcnt == 5 || abs (Q - Q0) > tol)
-      [Q1, fcnt, hmin] = simpsonstp (f, a, c, d, fa, fc, fd, Q1, fcnt, hmin,
+    ## Force at least one adpative step.
+    if (fcnt == 5 || abs (q - q0) > tol)
+      [q1, fcnt, hmin] = simpsonstp (f, a, c, d, fa, fc, fd, q1, fcnt, hmin,
 				    tol, trace, varargin{:});
-      [Q2, fcnt, hmin] = simpsonstp (f, c, b, e, fc, fb, fe, Q2, fcnt, hmin,
+      [q2, fcnt, hmin] = simpsonstp (f, c, b, e, fc, fb, fe, q2, fcnt, hmin,
 				     tol, trace, varargin{:});
-      Q = Q1 + Q2;
+      q = q1 + q2;
     endif
   endif
 endfunction