function vanderpol_period_test ( )

%*****************************************************************************80
%
%% vanderpol_period_test() tests vanderpol_period().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 December 2021
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'vanderpol_period_test():\n' );
  fprintf ( 1, '  vanderpol_period() returns values of\n' );
  fprintf ( 1, '  the period P of vanderpol_ode() as a function\n' );
  fprintf ( 1, '  of the parameter MU.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '      MU          P1        P2        P3        P4        P5\n' );
  fprintf ( 1, '\n' );

  for mu = 0 : 10

    [ p1, p2, p3, p4, p5 ] = vanderpol_period ( mu );

    fprintf ( 1, '  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f\n', ...
      mu, p1, p2, p3, p4, p5 );

  end

  for mu = 20 : 10 : 100

    [ p1, p2, p3, p4, p5 ] = vanderpol_period ( mu );

    fprintf ( 1, '  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f\n', ...
      mu, p1, p2, p3, p4, p5 );

  end

  for mu = 110 : 100 : 1000

    [ p1, p2, p3, p4, p5 ] = vanderpol_period ( mu );

    fprintf ( 1, '  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f\n', ...
      mu, p1, p2, p3, p4, p5 );

  end

  for mu = 2000 : 1000 : 10000

    [ p1, p2, p3, p4, p5 ] = vanderpol_period ( mu );

    fprintf ( 1, '  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f\n', ...
      mu, p1, p2, p3, p4, p5 );

  end

  return
end
