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