ElementRHS4ConPen.m

function fe = ElementRHS4ConPen(nnpe,nqptv,wtq,bigN,detj,data)
% computes elemental force matrix for body forces
%  (currently zeros)
%  Only essential BC's are applied in this code, so
% there are no integrations of surface tractions.

fe= zeros(nnpe,1);
feqp= zeros(nnpe,1);
bigNqp = zeros(nnpe,1);

  
for j=1:1:nqptv  
	
  xkfac=wtq(j)*detj(j);
  bigNqp(:) = bigN(:,j);
  feqp = xkfac*(bigNqp*data(j));
  fe = fe + feqp;

end