function E = anmEMtx(e)
% Generate material constituive matrix for a given element.
% The material constitutive matrix is assembled using the 
% element modulus of elasticity and poisson ratio.
% The matrix size and layout depends on the type of analysis.
% For plane stress and plane strain analysis, the size is (3,3).
% For axisymmetric analysis, the size is (4,4).
% Input arguments:
%  e:   index of element

 include_gblrefs;

 E_coeff = Elasticity(emat(e));          % get element elasticity coeff.
 P_ratio = Poisson(emat(e));             % get element poisson ration.

 if(analysis_type == AXISYMMETRIC)

  E = zeros(4,4);

  %%% COMPLETE HERE 14 %%%

 elseif(analysis_type == PLANE_STRAIN)

  E = zeros(3,3);

  %%% COMPLETE HERE 15 %%%

 else                                    % plane stress analysis

  E = zeros(3,3);

  %%% COMPLETE HERE 16 %%%

 end

end
