How can I define the angular velocity vectors for a mechanism to solve velocity loop closure equation on MATLAB?

Question

We have defined the thetas and solved the loop closure equations with "vpasolve" command, and to solve the velocity loop closure equations we have tried the simila