Avoid divide-by-zero in step_length_positive_definite.m.

[octave.git] / c_inner_product.m

diff --git a/c_inner_product.m b/c_inner_product.m
index 1d7bfd1e67122823724ebafb8ee919fe58c9980e..b53deef48adeba65439e81961ddf865dce4d7e60 100644 (file)
--- a/c_inner_product.m
+++ b/c_inner_product.m
@@ -1,24 +1,24 @@
 function ip = c_inner_product(w, a, b, v1, v2)
-  ##
-  ## The usual inner product defined on the space of continuous
-  ## functions over the interval [a,b].
-  ##
-  ## INPUT:
-  ##
-  ##   * ``w`` -- The weight function.
-  ##
-  ##   * ``a`` -- The left endpoint of the interval.
-  ##
-  ##   * ``b`` -- The right endpoint of the interval.
-  ##
-  ##   * ``v1`` -- The first vector.
-  ##
-  ##   * ``v2`` -- The second vector.
-  ##
-  ## OUTPUT:
-  ##
-  ## The inner product <v1, v2>.
-  ##
+  %
+  % The usual inner product defined on the space of continuous
+  % functions over the interval [a,b].
+  %
+  % INPUT:
+  %
+  %   * ``w`` -- The weight function.
+  %
+  %   * ``a`` -- The left endpoint of the interval.
+  %
+  %   * ``b`` -- The right endpoint of the interval.
+  %
+  %   * ``v1`` -- The first vector.
+  %
+  %   * ``v2`` -- The second vector.
+  %
+  % OUTPUT:
+  %
+  % The inner product <v1, v2>.
+  %
   integrand = @(x) w(x)*v1(x)*v2(x);
   ip = quad(integrand, a, b);
 end