% Function to get Gauss points 
% and weight factors
 
function [w,gp] = gauss(ngp)

 if ngp == 1
        gp = 0;
        w  = 2;
 elseif ngp == 2   % if four gauss points are used (two in each direction)
        gp = [-0.57735027,  0.57735027 ];    
        w  = [1,            1          ];   %
 end