Решение системы уравнений

CE*cos(fi3)+DE*sin(fi5)=CD*cos(fi3+fi4)
CE*sin(fi3)+DE*cos(fi5)=CD*sin(fi3+fi4)

B=solve('CE*cos(fi3)+DE*sin(fi5)=CD*cos((fi4+fi3))','CE*sin(fi3)+DE*cos(fi5)=CD*sin((fi3+fi4))','fi3','fi5')

syms fi3 fi5
CE=22; DE=22; CD=40; 
fi4=acos((CE^2+CD^2-DE^2)/(2*CE*CD));
B=solve(CE*cos(fi3)+DE*sin(fi5)==CD*cos((fi4+fi3)),CE*sin( fi3)+DE*cos(fi5)==CD*sin((fi3+fi4)),fi3,fi5)
B.fi3
B.fi5

Warning: The solutions are parametrized by the symbols:
z = (Dom::ImageSet(PI +
arccos((20*sin(7740781025443129/18014398509481984))/11) +
2*PI*k, k, Z_) union Dom::ImageSet(PI -
arccos((20*sin(7740781025443129/18014398509481984))/11) +
2*PI*k, k, Z_)) intersect (Dom::ImageSet(arcsin(-
(20*cos(7740781025443129/18014398509481984))/11 + 1) +
2*PI*k, k, Z_) union Dom::ImageSet(PI - arcsin(-
(20*cos(7740781025443129/18014398509481984))/11 + 1) +
2*PI*k, k, Z_))
 
> In solve at 190
  In NonLinSys at 6 
B = 
    fi3: [1x1 sym]
    fi5: [1x1 sym]
ans =
pi
ans =
z

clear, clc
 
CE=22;
DE=22 ;
CD=40 ;
fi4 = acos((CE^2+CD^2-DE^2)/(2*CE*CD));
 
% u(1) -> fi3
% u(2) -> fi5
F = @(u) [CE*cos(u(1))+DE*sin(u(2))-CD*cos(u(1)+fi4);
          CE*sin(u(1))+DE*cos(u(2))-CD*sin(u(1)+fi4)];
u0 = [0.1 0.1]; % начальное приближение
R = fsolve(F, u0) % находим численное решение
 
fi = rad2deg(R) % ответ в градусах

clear, clc
 
CE=22;
DE=22 ;
CD=40 ;
fi4 = acos((CE^2+CD^2-DE^2)/(2*CE*CD));
 
% u(1) -> fi3
% u(2) -> fi5
F = @(u) [CE*cos(u(1))+DE*sin(u(2))-CD*cos(u(1)+fi4);
          deg2rad(4)+2*u(1)+u(2)-90];
u0 = [0.1 0.1]; % начальное приближение
R = fsolve(F, u0) % находим численное решение
 
fi = rad2deg(R) % ответ в градусах