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

<p dir="ltr">%This code implements an ODE system for anthocyanin monomers A^-, A, AH^+, B and the self-association of A^-, A, AH^+. </p><p dir="ltr">%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.</p><p dir="ltr">%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.</p>