function lbubble=GetNodalFunctionsBubbleStokes(pentaelem,gauss_coord)
%GETNODALFUNCTIONSBUBBLESTOKES - return the value of the bubble nodal function at the gaussian point
%
%   Usage:
%      lbubble=GetNodalFunctionsBubbleStokes(pentaelem,gauss_coord)
%
%   See also GETNODALFUNCTIONSDERIVATIVESBASICSTOKES, GETNODALFUNCTIONSDERIVATIVESPARAMSTOKES

lbubble=27*gauss_coord(1,1)*gauss_coord(1,2)*gauss_coord(1,3)*(1-gauss_coord(1,4))*(1+gauss_coord(1,4));
