forced oscillations code snippets in matlab

phase as a function of frequency for forced oscillations in matlab