function driver

  fname = 'f16b';
  a = 0;
  b = 1;

  driverf (fname, a, b);

