function [weight,point]=quad_3()
% 3 point quadrature on the reference triangle (2nd order exactly)
w1=1/3; weight=[w1,w1,w1];
point(:,1)=[0;1/2]; point(:,2)=[1/2;0]; point(:,3)=[1/2;1/2];