function y = eqlbrm(x,a,epsk,c)
% y = eqlbrm(x,a,epsk,c)
% This is the function for the malate reaction equations.
y = [ x(1)*(1.0 + c*x(3)/epsk) - a(1);
      x(2)*(1.0 + x(3)/epsk^2) - a(2);
      x(3)*(1.0 + c*x(1)/epsk + x(2)/epsk^2) - a(3)];