-- triquad.ads    compute integration coordinates and weights for order 1..6
--                Gauss Legendre over triangle
--                1 <= n <=6 at this time
--                nn returned is n*n
--                integral = sum i=0 to nn-1 W(i)*f(X(i),Y(i))

with Real_Arrays; -- types real and real_vector and real_matrix
use Real_Arrays;

package Triquad is

  procedure Triquad_Msg;

  procedure tint(n:Integer;  T:Real_Vector; nn:in out Integer;
                 X:out Real_Vector; Y:out Real_Vector; W:out Real_Vector);
end Triquad;