Numerical simulation of magmatic system evolution

This MATLAB code simulates the evolution of the magma system as presented in article "The Bimodal Road to Supereruptions". Beginning with the temperature and pressure conditions characteristic of a smaller-scale eruption, the model incorporates key processes (including crystal-melt separation and mafic magma injection) to reproduce the conditions leading to a supereruption. The code generates the results corresponding to Figure 3 of the article.