function rs = blood_pressure_resistance ( )

%*****************************************************************************80
%
%% blood_pressure_resistance() returns the resistance of a blood network.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 March 2022
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real RS(:), the resistance.
%
  [ pdia, psys, pulse, q, told, t0, y0, tstop ] = blood_pressure_parameters ( );

  T = 1.0 / pulse;
  ca = blood_pressure_compliance ( );
  rs = T / ( ca * ( log ( psys ) - log ( pdia) ) );

  return
end
