forced oscillation code snippets in matlab

phase as a function of frequency for forced oscillationz in matlab