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)");