function BeamElementShearDiagram(f, L)
% BeamElementShearDiagram

x = [0; L];
z = [f(1); -f(3)];

plot(x,z,x,zeros(2,1),'k');
title('Shear Force Diagram');

endfunction