-- galke.adb for fem_checke_la.adb  Galerkin FEM
with Ada.Numerics.Long_Elementary_Functions;
use Ada.Numerics.Long_Elementary_Functions;
with Real_Arrays; -- types real, real_vector, real_matrix
use  Real_Arrays;
with laphi;       -- phi, phip, phipp, phippp, phipppp
use laphi;

  -- the following two functions must be external for aquad3 using 'access
function galke(x:Real; i,j,n:Integer; Xg:Real_vector) return real is
                                           -- Galerkin k stiffness function
begin                                      -- for this specific problem
  return phip(x,j,n,xg)*phi(x,i,n,xg);
end galke;