Не получается использовать solve().

close all
clear all
 
%% известные величины
G(1)=-15*pi/180; %% начальный угол качки, рад
dg=0.001; %%приращение угла G, рад
i=1; %%начало отсчёта
n=510; %%окончание отсчёта
R1=sym(3.2); %%длинна выносной балки, м
R2=sym(0.79); %%высота стойки, м
g0=sym(15*pi/180); %% базовый угол наклона относительно центра качания, рад
f=sym(10*pi/180); %% угол возвышения балки, рад
L=sym(5); %% расстояние от центра качания до точки установки устройства, м
 
%% решение системы уравнений
syms a b g
[a,b]=solve('L*(sin(g0+g)-sin(g0))-R1*(cos(f)+cos(a))-R2*sin(b)=0', 'L*(cos(g0)-cos(g+g0))-R1*(sin(a)-sin(f))+R2*(1-cos(b))=0');
 
%% поиск углов поворота
g=G(1);
f1(1)=a+b;
f2(1)=g+b;
for i=2:n
    g=G(i-1);
    f1(i)=a+b;
    f2(i)=g+b;
end
 
%%построение графика
figure;
plot(G,f1,'--r',G,f2,'-.b','linewidth',2)
grid on;
legend('f1, град','f2, град') 
xlabel('угол гамма, град')
ylabel('углы f1,f2, град')