Here you find a Bloch-McConnell simulation introduced in the publications

• Schmitt et al. (2011), MRM 65: 1620–1629. DOI: 10.1002/mrm.22750,
• Roeloffs et al. (2014), NMR Biomed., 28, 40–53 doi: 10.1002/nbm.3192
• Zaiss et al. (2014), NMR Biomed., 28; 217–230 doi: 10.1002/nbm.3237,
• Meissner et al. (2015), NMR Biomed., 28, 1196–1208. doi: 10.1002/nbm.3362,

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

• multiple pools possible (water, 5 CEST, 1 semisolid MT)
• trains RF pulses with arbitratry pulse shapes possible (gauss, sinc, spin-lock already implemented)
• both analytical (if possible) and full numerical simulation possible
• transient and steady-state possible
• first order readout schemes possible
• Both full Z-spectra and single offset evaluation possible
• Ready for multi parameter simulations

Read more about usage on github in BM_Documentation.docx and BM_tutorial.pptx

This code was created by Moritz Zaiss based on the cw BM simulation and fit code of Shanrong Zhang. A revised version was created with the help of Patrick Schuenke, Christian Meyer, Christian David and Volkert Roeloffs.