Using MATLAB R2017b:

Write a program, plotharmonic.m, that plots harmonic motion of x in time, as described by the following equation: x(t) = A cos(t + ) In this equation, A is the amplitude of the motion, is the angular frequency, and is the phase shift, all of which are constants in the equation. Run the program with the values A = 2, and = 2 . Pick a value for the period T that is related to by the relation = 2/T . Also pick a good value for the final time Tf that sets the range of t = [0, ... , Tf ].