%
%   GFSERIES1B.M   Called by quadl to determine ym*yn integrals
%
%   File prepared by J. R. White, UMass-Lowell (Nov. 2003)
%
      function f = ifile(x)
      global h L ALPM ALPN
      f = sin(ALPM*x).*sin(ALPN*x);
%
%   end of function
