PROGRAM d4r4(input,output);
(* driver for routine QROMB *)
CONST
   pio2=1.5707963;
   n=5;
TYPE
   glnarray = ARRAY [1..n] OF real;
VAR
   glit : integer;
   a,b,s : real;

FUNCTION func(x: real): real;
(* Test function *)
BEGIN
   func := sqr(x)*(sqr(x)-2.0)*sin(x)
END;

FUNCTION fint(x: real): real;
(* Integral of test function func *)
BEGIN
   fint := 4.0*x*(sqr(x)-7.0)*sin(x)
      -(sqr(sqr(x))-14.0*sqr(x)+28.0)*cos(x)
END;

(*$I MODFILE.PAS *)
(*$I TRAPZD.PAS *)

(*$I POLINT.PAS *)

(*$I QROMB.PAS *)

BEGIN
   a := 0.0;
   b := pio2;
   writeln ('Integral of func computed with QROMB');
   writeln;
   writeln ('Actual value of integral is',fint(b)-fint(a):12:6);
   qromb(a,b,s);
   writeln ('Result from routine QROMB is',s:12:6);
END.