c_inner_product.m

   1 function ip = c_inner_product(w, a, b, v1, v2)
   2   %
   3   % The usual inner product defined on the space of continuous
   4   % functions over the interval [a,b].
   5   %
   6   % INPUT:
   7   %
   8   %   * ``w`` -- The weight function.
   9   %
  10   %   * ``a`` -- The left endpoint of the interval.
  11   %
  12   %   * ``b`` -- The right endpoint of the interval.
  13   %
  14   %   * ``v1`` -- The first vector.
  15   %
  16   %   * ``v2`` -- The second vector.
  17   %
  18   % OUTPUT:
  19   %
  20   % The inner product <v1, v2>.
  21   %
  22   integrand = @(x) w(x)*v1(x)*v2(x);
  23   ip = quad(integrand, a, b);
  24 end