Rheological properties of reconstituted skim milk coagulated with plant coagulants Cynara cardunculus L., Cynara humilis L. and chymosin was monitored by dynamic low amplitude oscillation. There are no published reports on the modelling of the gelation behaviour of milk by plant coagulants. Three mathematical models, Scott Blair, Douillard and Carlson, were fitted to the storage modulus (G′) as function of time curves. For all coagulants, Scott Blair model was the most efficient in modelling the gelation process, and gave both the smallest residuals and standard error of residuals, Se (P < 0·0001). Douillard model gave the poorest fit and in particular it was not able to predict the initial part of the gelation curves. Carlson model had an intermediate behaviour and, in the case of chymosin, it gave results that were quite comparable to Scott Blair model. The parameters of the Scott Blair model were different for plant coagulants and chymosin. Chymosin had the longest rate constant (τ) and the time shift coefficient (ts) was also different between vegetable coagulants and chymosin (P < 0·0005). These results are in agreement with the overall trends for gelation profiles obtained for vegetable coagulants and chymosin. In the beginning of gelation both plant coagulants produced gels with slightly higher G′ values than chymosin, but after longer incubation times chymosin gels had higher G′ values. It was concluded that the Scott Blair model was the best equation to follow the gelation of milk induced by both plant coagulants as well as chymosin. Modelling is an important and useful method for comparing the gelation process in gels formed by different types of coagulants.