function Rtot = U_RL_model_numeric(params, constants )
%------------------- LineshapeKin Simulation -----------------------------
% Function for numeric solving of equilibrium thermodynamic equations
%
% 5/12/09
% Evgenii Kovrigin
% Medical College of Wisconsin
%
% U-RL model developed in 
% /Users/kovrigin/Documents/Workspace/Data/Data.XV/EKM16.Analysis_of_multis
% tep_kinetic_mechanisms/Equilibria/LRIM_U_RL.mn




%--- this is what I need for K-matrix
Leq=params(1);


%-- these are constants 
K_a_A=constants(1);
K_a_B=constants(2);
Ltot=constants(3);



%-- equation

Rtot = - (Leq - Ltot)/(K_a_B + 1) - (K_a_B*(Leq - Ltot))/(K_a_B + 1) - (Leq - Ltot)/(K_a_A*Leq*(K_a_B + 1));
