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.
ARC Training Centre for Innovative Wine Production
Australian Research CouncilFind out more...