function f = BeamElementForces(k,u)
% BeamElementForces

f = k * u;

endfunction