function [num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights, vert_gauss_coord,vert_gauss_weights]=GaussPenta(area_order,num_vert_gauss)
%GAUSSPENTA - get Gauss point for penta elements
%
%   This routine calls the routines GAUSSTRIA and GAUSSSEGMENT
%
%   Usage:
%      [num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights, vert_gauss_coord,vert_gauss_weights]=GaussPenta(area_order,num_vert_gauss)

%Get Tria area gauss points
[num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights]=GaussTria(area_order);

%Get vertical segment gauss points
[vert_gauss_coord vert_gauss_weights]=GaussSegment(num_vert_gauss);
