%
%   CHIP2D_SOV1B.M   Called by quadl to determine In integrals
%                            (used with CHIP2D_SOV1.M)
%
%   File prepared by J. R. White, UMass-Lowell (Nov. 2003)
%
      function I = ifile(x)
      global hs kc Lc Q EIGN EIGM
      W = Q*Lc/(2*hs) + (Q/2/kc)*((Lc/2)^2 - x.^2);                                    
      I = W.*cos(EIGN*x);
%
%   end of function
