Here you find a Bloch-McConnell simulation introduced in the publications
It was initially based on the code provided by Shanrong Zhang published in Woessner, Zhang, Merritt and Sherry: Numerical solution of the Bloch equations provides insights into the optimum design of PARACEST agents for MRI, MRM 2005. It was extended for multiple pools, including semi-solid MT, and also extended for pulsed RF saturation.
Dowload zipped Matlab implementations here or find the package on https://github.com/cest-sources/BM_sim_fit/
This code allows for simulation of BM system upon RF presaturation
Read more about usage on github in BM_Documentation.docx and BM_tutorial.pptx