MATLAB: How to derive an equation of Theta_a in terms of L_a using ‘solve’ command or any other way in Matlab

equationsolve

aIa=1;
bIa=aIa;
aIIa=1.25*aIa;
landaIa=pi/3;
landaIIa=acos((aIa*cos(landaIa))/aIIa);
KIa=1;
Kca=KIa;
KIIa=20*KIa;
s1a=sin(landaIa);
t2a=tan(landaIIa);
t1a=tan(landaIa);
syms L_a theta_a
eqn1_a=aIa*s1a*((sqrt(((t2a^2)/(t1a^2))-((cos(theta_a).^2))))-sin(theta_a))-L_a==0;
g=solve(eqn1_a,L_a)

Best Answer

There appear to be two solutions:
atan2(-(16*L_a^2-9)*sqrt(3)/(48*L_a), sqrt(-768*L_a^4+3168*L_a^2-243)/(48*L_a))
atan2(-(16*L_a^2-9)*sqrt(3)/(48*L_a), -sqrt(-768*L_a^4+3168*L_a^2-243)/(48*L_a))
Found with Maple. MATLAB does not appear to be able to find this.