function [Kint,Mint,Fint]=splninteg(BX1,BX2,DB1,DB2,W,X)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% find the finite element stifness, mass, load coeff.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
global px qx fx
Kint=sum(W.*px(X).*DB1.*DB2);
Mint=sum(W.*qx(X).*BX1.*BX2);
Fint=sum(W.*fx(X).*BX1);