как найти состояние равновесия аналитическим способом

function Payk_model
function Payk_model
%модель Qd-модель функции спроса 
Pt=[0.5:0.01:20]; 
a = 50; 
b = 1.2;
Qd = a-b*(Pt-1); 
%модель Qs-модель функции предложения
Pt=[0.5:0.01:20]; 
c = 23;
d = 2.8;
Qs = c+d*Pt-1;
V1 = Qd; V2 = Qs; 
% рисуем график(угол наклона спроса меньше угла наклона предложения)
figure;
plot(V1,Pt,V2,Pt,'LineWidth',2); grid on; 
title('Паутинообразная модель равновесия (b<d)'); 
xlabel('Объём предложения Q'); ylabel('фактическая цена экономического блага P');
legend('Qs','Qd');
 
 
%модель Qd-модель функции спроса
Pt=[0:0.01:20]; 
a = 50; 
b = 1.5;
Qd = a-b*(Pt-1); 
%модель Qs-модель функции предложения
Pt=[0:0.01:20]; 
c = 23;
d = 1.5;
Qs = c+d*Pt-1;
V3 = Qd; V4 = Qs; 
% рисуем график(b=d)
figure;
plot(V3,Pt,V4,Pt,'LineWidth',2); grid on; 
title('Паутинообразная модель равновесия, (b=d)'); 
xlabel('Объём предложения Q'); ylabel('фактическая цена экономического блага P');
legend('Qd','Qs');
 
%модель Qd-модель функции спроса
Pt=[0.5:0.01:20]; 
a = 50; 
b = 2.6;
Qd = a-b*(Pt-1); 
%модель Qs-модель функции предложения
Pt=[0.5:0.01:20]; 
c = 23;
d = 1.2;
Qs = c+d*Pt-1;
V5 = Qd; V6 = Qs; 
% рисуем график (b>d)
figure;
plot(V5,Pt,V6,Pt,'LineWidth',2); grid on; 
title('Паутинообразная модель равновесия, (b>d)'); 
xlabel('Объём предложения Q'); ylabel('фактическая цена экономического блага P');
legend('Qd','Qs');