Write a function in Matlab that plots harmonic motion of x in time, as described by the following equation

x(t)=Acos(wt + phi)

In this equation, A is the amplitude of the motion, w is the angular frequency, and phi is the phase shift, all of which are defined by the user.

Write a second function when another harmonic motion in the y direction is also present:

y(t)= Bcos(2wt - phi)