%
%   SOV_EX2B.M   Called by quadl to determine yn*yn integrals - which
%                are the norm of the eigenfunctions (used in SOV_EX2.M)
%
%   File prepared by J. R. White, UMass-Lowell (Nov. 2004)
%
      function f = ifile(x)
      global hok L ALPN uL uo CC1
      f = sin(ALPN*x).*sin(ALPN*x);
%
%   end of function
