function I = ingegrandSolid1DExample(xi)
if nargin < 1
    xi = -1/sqrt(3);
end

global x0;
global x1;
global x2;
global E;
global A;
global L;
global xAve;
D = E * A;

B = [(2 * xi - 1) / 2  (2 * xi + 1) / 2 -2 * xi];
J = (xAve - x2) * (2 * xi) + L / 2.0;
I = 1 / J * D * transpose(B) * B;