Refix the non-negativity of denom in step_length_positive_definite().

[octave.git] / c_inner_product.m

diff --git a/c_inner_product.m b/c_inner_product.m
index 0ff7f9631dee488c7365d8d4df7676fca131c0e4..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>.
-  ##
-  integrand = @(x) w(x)*v1(x)*v2(x)
-  ip = quad(integrand, a, b)
+  %
+  % 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