function l1l6=GetNodalFunctions(pentaelem,gauss_coord)
%GETNODALFUNCTIONS - return the values of nodal functions at gaussian point on the reference element
%
%   Usage:
%      l1l6=GetNodalFunctions(pentaelem,gauss_coord)
%
%   See also GETNODALFUNCTIONSSTOKES, GETNODALFUNCTIONSDERIVATIVESPARAMS, GETNODALFUNCTIONSDERIVATIVESBASIC

l1l6=[gauss_coord(1,1)*(1- gauss_coord(1,4))/2,  gauss_coord(1,2)*(1- gauss_coord(1,4))/2,  gauss_coord(1,3)*(1- gauss_coord(1,4))/2, ...
    gauss_coord(1,1)*(1+ gauss_coord(1,4))/2,  gauss_coord(1,2)*(1+gauss_coord(1,4))/2,  gauss_coord(1,3)*(1+ gauss_coord(1,4))/2];
