0001 function [num_gauss,first_gauss_coord,second_gauss_coord,gauss_weights]=GaussQuad(order)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 if order==1,
0015 num_gauss=1;
0016 gauss_weights=4;
0017 first_gauss_coord= 0;
0018 second_gauss_coord= 0;
0019
0020
0021 elseif order==2,
0022 num_gauss=4;
0023 gauss_weights=[ 1, 1, 1, 1];
0024 first_gauss_coord=[ -0.57735026918962, 0.57735026918962, -0.57735026918962, 0.57735026918962];
0025 second_gauss_coord=[0.57735026918962, 0.57735026918962, -0.57735026918962, -0.57735026918962];
0026
0027
0028 elseif order==3,
0029 num_gauss=9;
0030 gauss_weights=[ 0.30864197530864, 0.49382716049383, 0.30864197530864, 0.49382716049383, 0.79012345679012, 0.49382716049383, 0.30864197530864, 0.49382716049383, 0.30864197530864];
0031 first_gauss_coord=[ -0.77459666924148, 0, 0.77459666924148, -0.77459666924148, 0, 0.77459666924148, -0.77459666924148, 0, 0.77459666924148];
0032 second_gauss_coord=[ -0.77459666924148, -0.77459666924148, -0.77459666924148, 0, 0, 0, 0.77459666924148, 0.77459666924148, 0.77459666924148];
0033
0034 else
0035 error('GaussQuad error message: order not supported yet');
0036 end