Index: /issm/trunk/src/c/shared/Numerics/GaussPoints.cpp
===================================================================
--- /issm/trunk/src/c/shared/Numerics/GaussPoints.cpp	(revision 3465)
+++ /issm/trunk/src/c/shared/Numerics/GaussPoints.cpp	(revision 3466)
@@ -2041,4 +2041,38 @@
 		
 	}
+	else if (num_gauss==4){
+		//order=3, num_gauss=4. Can integrate polynomials of degree 0 to 5
+
+		gauss_weights=(double*)xmalloc(num_gauss*sizeof(double));
+		gauss_weights[0]=0.347854845137454;
+		gauss_weights[1]=0.652145154862546;
+		gauss_weights[2]=0.652145154862546;
+		gauss_weights[3]=0.347854845137454;
+		
+		segment_coord=(double*)xmalloc(num_gauss*sizeof(double));
+		segment_coord[0]=-0.861136311594053;
+		segment_coord[1]=-0.339981043584856;
+		segment_coord[2]=0.339981043584856;
+		segment_coord[3]=0.861136311594053;
+		
+	}
+	else if (num_gauss==5){
+		//order=3, num_gauss=4. Can integrate polynomials of degree 0 to 5
+
+		gauss_weights=(double*)xmalloc(num_gauss*sizeof(double));
+		gauss_weights[0]=0.236926885056189;
+		gauss_weights[1]=0.478628670499366;
+		gauss_weights[2]=0.568888888888889;
+		gauss_weights[3]=0.478628670499366;
+		gauss_weights[4]=0.236926885056189;
+		
+		segment_coord=(double*)xmalloc(num_gauss*sizeof(double));
+		segment_coord[0]=-0.906179845938664;
+		segment_coord[1]=-0.538469310105683;
+		segment_coord[2]=0.0;
+		segment_coord[3]=0.538469310105683;
+		segment_coord[4]=0.906179845938664;
+		
+	}
 	else{
 		_printf_("GaussSegment error message: order not supported yet");