function value = stiff_solution ( t )

%*****************************************************************************80
%
%% stiff_solution() evaluates the exact solution of the stiff ODE.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    07 February 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real T: the evaluation time.
%
%  Output:
%
%    real VALUE, the solution at the given time.
%
  value = 50.0 * ( sin ( t ) + 50.0 * cos ( t ) - 50.0 * exp ( - 50.0 * t ) ) / 2501.0;

  return
end

