MATLAB code to model anthocyanin proton transfer, hydration and self-association

%This code implements an ODE system for anthocyanin monomers A^-, A, AH^+, B and the self-association of A^-, A, AH^+.

%Aggregation is implemented as an isodesmic %model; the aggregation rate vectors jp(1,2,3)v, jn(1,2,3)v have constant entries. The code can easily be modified for nonisodesmic models by choosing aggregation rate vectors with nonconstant entries.

%The constants used are for malvidin. The model can easily be modified to describe kinetic behaviour for all anthocyanins by assigning the rate constants kp, kn, jp and jn appropriately.