Implement mupad Codes into Matlab Script

Can anyone help me to implement the following mupad codes into matlab script so that I can directly save the result into the working space?

l1 := 0.167;
l2 := 1.136037;
THETA:= -2.0;

argu:=sqrt(k)*l1;

num  := 1/sqrt(k)*sin(argu)+l2*cos(argu);
den := -l2*sqrt(k)*sin(argu)+cos(argu);

F:=-THETA+arctan(num/den)

numeric::fsolve(F, [k=-40..40], MultiSolutions);

Thanks a lot!

Answer 1

[Walter Roberson]

l1 = 0.167;
l2 = 1.136037;
THETA = -2.0;
syms k
argu = sqrt(k)*l1;

num  = 1/sqrt(k)*sin(argu)+l2*cos(argu);
den = -l2*sqrt(k)*sin(argu)+cos(argu);

F = -THETA+arctan(num/den);
vpasolve(F, k, [-40,40])

With regards to multiple solutions, see http://www.mathworks.com/help/symbolic/vpasolve.html#zmw57dd0e120083

However, there are no real-valued solutions to that equation.

[Karry Li]

Thanks a lot. Only one question: how do you come up with vpasolve instead of fzero and fsolve?