PROGRAM d10r1(input,output);
(* driver for routine MNBRAK *)
VAR
   ax,bx,cx,fa,fb,fc : real;
   i : integer;

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

FUNCTION func(x: real): real;
BEGIN
   func := bessj0(x)
END;

(*$I MNBRAK.PAS *)

BEGIN
   FOR i := 1 to 10 DO BEGIN
      ax := i*0.5;
      bx := (i+1.0)*0.5;
      mnbrak(ax,bx,cx,fa,fb,fc);
      writeln('a':14,'b':12,'c':12);
      writeln('x':3,'  ',ax:12:6,bx:12:6,cx:12:6);
      writeln('f':3,'  ',fa:12:6,fb:12:6,fc:12:6)
   END
END.