function blood_pressure_exact_test ( )

%*****************************************************************************80
%
%% blood_pressure_exact_test() displays the exact solution of blood_pressure_ode().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 March 2022
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'blood_pressure_exact_test():\n' );
  fprintf ( 1, '  blood_pressure_exact() evaluates the exact solution\n' );
  fprintf ( 1, '  to blood_pressure_ode()\n' );
%
%  Get the parameters.
%
  [ pdia, psys, pulse, q, told, t0, y0, tstop ] = blood_pressure_parameters ( );
%
%  Evaluate the solution at 501 points.
%
  t = linspace ( t0, tstop, 501 );

  ye = blood_pressure_exact ( t );
%
%  Plot the exact solution.
%
  plot ( t, ye, 'r-', 'linewidth', 3 );

  title ( 'blood\_pressure\_ode():, exact' );
  grid ( 'on' );
  xlabel ( '<--- T --->' );
  ylabel ( '<--- Y(T) --->' );
%
%  Save a copy in a file.
%
  filename = 'blood_pressure_exact.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end

