function Simple(varargin)
% Create and configure a ChemPlugin instance
cp = ChemPlugin();
Config(cp, 'Na+ = 1 mmol/kg; Cl- = 1 mmol/kg; pH = 3');
Config(cp, 'react 2 mmol NaOH; delxi = 0.1');

% Calculate initial conditions
Initialize(cp);
fprintf('pH = %f\n', Report1(cp, 'pH'));

% March forward in time.
while 1
    deltat = ReportTimeStep(cp);
    if AdvanceTimeStep(cp, deltat)
        break;
    end
    if AdvanceChemical(cp)
        break;
    end
    fprintf('pH = %f\n', Report1(cp, 'pH'));
end

end
