function intB = GaussIntB(h)
intB = zeros(6,24);
x0 = 0; y0 = 0; z0=0;
numgp=2;
gpx = getgp(x0,x0+h,numgp);
gpy = getgp(y0,y0+h,numgp);
gpz = getgp(z0,z0+h,numgp);
for ix = 1:numgp*2
    x=gpx(ix);
    for iy = 1:numgp*2 
        y=gpy(iy);
        for iz = 1:numgp*2
            z=gpz(iz);
            intB = intB + LocalB(x,y,z,x0,y0,z0,h);
        end
    end
end
intB = intB.*h^3/(numgp*2)^3;
end