GRAPHS_SRAVNENIYE

clear all;
close all;
F0 = 0.01;
dg = 0:0.1:5;
Dd = zeros(51:1);
% Ddexp = zeros(51:1);
Dd1 = 0.8;
Ps1 = 0.5*(1-Dd1+F0);
for i = 1:51
    Dd(i) = func_D(dg(i),F0);
    Ddexp(i) = func2_exp_const_thr_calc(F0,dg(i));
end
dg1 = 3.168;
% dg2 = 19.6377;


figure(228);
hold on;
p = plot(dg1, Dd1);
p(1).Marker = '.';p(1).MarkerSize = 30;
% p = plot(dg2, Dd1);
% p(1).Marker = '.';p(1).MarkerSize = 30;

p = plot(dg, Dd);
p(1).LineWidth = 2;

Dp = [0.0537	0.0979	0.2024	0.3045	0.3907	0.4948	0.5910	0.7047	0.793 0.9075 0.9554	0.991];
dp = [0.6815 1.0448	1.4847	1.8019	2.0730	2.3263	2.5797	2.8507	3.168 3.6079 3.9712	4.6527];
p = plot(dp,Dp);
p(1).LineWidth = 1; p(1).Marker = 'o';p(1).MarkerSize = 6;p(1).Color = "blue";

 p = plot(dg, Ddexp);
 p(1).LineWidth = 2;

% Dpexp = [0.0554	0.0958	0.2262	0.3049	0.3907	0.4977	0.5973	0.7001	0.7988	0.9059];
% dpexp = [0.5372	1	1.8614	2.825	4.0259	5.6439	8.0152	11.9114	19.64	42.7087];
% p = plot(dpexp,Dpexp);
% p(1).LineWidth = 1; p(1).Marker = 'o';p(1).MarkerSize = 6;p(1).Color = "blue";


grid on;
xlabel("d");
ylabel("D(d)");





Ps = 0.5*(1-Dd+F0);
Psp = 0.5*(1-Dp+F0);
figure(1337);
hold on;
p = plot(dg1, Ps1);
p(1).Marker = '.';p(1).MarkerSize = 30;

p = plot(dg, Ps);
p(1).LineWidth = 2;p(1).Color = "red";
p = plot(dp, Psp);
p(1).LineWidth = 1; p(1).Marker = 'o';p(1).MarkerSize = 6;p(1).Color = "blue";
grid on;
xlabel("d");
ylabel("Ps(d)");