function Int=FEintergrate2D(CtrlVar,MUA,f)

Int=FEintergrate2D(CtrlVar,MUA,f)
calculates the integral of f over each of the elements of the FE mesh given in the structure MUA
CtrlVar can be entered as an empty variable

Example. Calculate volume of ice:

Int=FEintergrate2D([],MUA,h);
TotalIceVolume=sum(Int)

ndim=2;

fnod=reshape(f(MUA.connectivity,1),MUA.Nele,MUA.nod);

[points,weights]=sample('triangle',MUA.nip,ndim);

Int=zeros(MUA.Nele,1);

for Iint=1:MUA.nip

    fun=shape_fun(Iint,ndim,MUA.nod,points) ;
    detJ=MUA.DetJ(:,Iint);
    detJw=detJ*weights(Iint);
    fint=fnod*fun;
    Int=Int+fint.*detJw;

end

end

Published with MATLAB® R2016b