Code to call your function $(3)$

$\%$ Define an array theta between $0$ and $2**$ pi with $181$ elements using linear spacing

theta $=$linspace$(0,2**,181)$;

$\%$ Define a variable $R$ with value $3.$ This is the magnitude of the vector $r$

$R=3$;

$\%$ Define a row array of three columns for angular velocity. Set the last element to $3$

$W=\left[\begin{array}{ccc}0& 0& 3\end{array}\right]$;

$\%$ Define a row array of three columns for angular acceleration. Set all the elements equal to zero.

$A=\left[\begin{array}{ccc}0& 0& 0\end{array}\right]$;

$\%$ Define a row array of three columns for reference velocity. Set the first element to $-1,$ and second element to $3.$

v$\_$ref $=\left[\begin{array}{ccc}-1& 3& 0\end{array}\right]$;

$\%$ Define a row array of three columns for reference acceleration. Leave all values at zero.

a$\_$ref $=\left[\begin{array}{ccc}0& 0& 0\end{array}\right]$;

$\%$ Initiate a zeros array for vector $r.$ Number of rows: length of theta, Number of columns: $3$

length$($theta$),3)$;

$\%$ Initiate a zeros array for vector $v$ with the same dimensions as $r.$

$v=$ zeros$($length$($theta$),3)$;

$\%$ Initiate a zeros array for vector a with the same dimensions as $r.$

$a=$ zeros $($length$($theta$),3)$;

$\%$ In a for loop, determine the three components of $r$ for each theta.

$\%$ No need to create plots for this problem.