You can put something along these lines in a script:
subplot(2,2,1)
text(0.02,0.98,'a','Units', 'Normalized', 'VerticalAlignment', 'Top')
subplot(2,2,2)
text(0.02,0.98,'b','Units', 'Normalized', 'VerticalAlignment', 'Top')
subplot(2,2,3)
text(0.02,0.98,'c','Units', 'Normalized', 'VerticalAlignment', 'Top')
subplot(2,2,4)
text(0.02,0.98,'d','Units', 'Normalized', 'VerticalAlignment', 'Top')
Please note that I am assuming that your subheadings are located in a 2x2 grid and that they are all 2D graphics. If these assumptions are not fulfilled, change the first two arguments of the subtask and / or add the z coordinate to the text.
source
share