function f = GetEquivalentNodalLoad(w,L,dir)
% f = GetEquivalentNodalLoad(w,L,dir)
%
% GetEquivalentNodalLoad
%
%   f = equivalent nodal load;
%   w = distributed load;
%   L = beam length;
%   dir = {'2' for y direction|'3' for z direction}.
%

% Walter Ponge-Ferreira
% Cotia, 17.06.2018

if nargin<3,
  dir = 2;
endif
if nargin<2,
  L = 1;
endif

if dir == 2,
  f = [w*L/2; w*L^2/12; w*L/2; -w*L^2/12];
elseif dir == 3,
  f = [-w; 0; w; 0];      
else
  error('Incorrect load direction');
endif

endfunction