有限元高斯积分点MATLAB

function [r w] = Gauss（l）
% RETURN COORDINATES OF GAUSS QUADRATURE POINTS AND THEIR WEIGHTS
%
% SYNTAX: [r w] = Gauss（l）
%
% INPUT PARAMETERS
%    l - Order of Gauss quadrature （ 1 =< r <= 10 ）
%
% OUTPUT PARAMETERS
%    r - Coordinates of Gauss quadrature points
%    w - Weights at Gauss quadrature points
%
xg = [
               0.0              0.0              0.0              0.0              0.0              0.0              0.0              0.0              0.0              0.0
  -0.5773502691896  0.5773502691896              0.0              0.0              0.0	            0.0              0.0              0.0              0.0              0.0
  -0.7745966692415  0.0000000000000  0.7745966692415              0.0              0.0	            0.0              0.0              0.0              0.0              0.0 
  -0.8611363115941 -0.3399810435849  0.3399810435849  0.8611363115941              0.0	            0.0              0.0              0.0              0.0              0.0
  -0.9061798459387 -0.5384693101057              0.0  0.5384693101057  0.9061798459387              0.0              0.0              0.0              0.0              0.0 
  -0.9324695142032 -0.6612093864663 -0.2386191860832  0.2386191860832  0.6612093864663  0.9324695142032              0.0              0.0              0.0              0.0 
  -0.9491079123428 -0.7415311855994 -0.4058451513774              0.0  0.4058451514774	0.7415311855994  0.9491079123428              0.0              0.0              0.0
  -0.9602898564975 -0.7966664774136 -0.5255324099163 -0.1834346424957  0.1834346424957	0.5255324099163  0.7966664774136  0.9602898564975              0.0              0.0 
  -0.9681602395076 -0.8360311073266 -0.6133714327006 -0.3242534234038              0.0	0.3242534234038  0.6133714327006  0.8360311073266  0.9681602395076              0.0
  -0.9739065285172 -0.8650633666890 -0.6794095682990 -0.433395394129