function system2
    clc; clear;
    global s S D v e;
    %кин константы
    %изменяемая константа комментируется
    s = 1;
    S = 1;
%   D = 1;
    v = 0;
    e = 1;
    %литературное
    T_hl = 9; %время полураспада
    T_max = 2.5; %в минутах, время максимальной концентрации 

    %нач приближение
    A0 = 62; %мкг/мл %литературное
    B0 = 0;
    y0 = [A0 B0];
    
    rtol = 1.e-9;%6%11%12
    atol = 1.e-11;%7%12%15
    rfn = 5;
    options = odeset('BDF','ON','RelTol',rtol,'AbsTol',atol,'Refine',rfn);
     
    tspan = [0:0.1:6.7*T_hl]; %6.7

    figure(7);  hold on; grid on;
    ylabel('C[A], C[B], calculated');
	xlabel('t');
	title('C[A] - black, C[B] - red    —    t' );
    
    %цикл ведется по изменяемой константе
    for D = [0.1 0.5 1 10]
        %решающая функция
        [t,y] = ode15s(@f, tspan, y0, options); 
        plot(t,y(:, 1), 'k', 'LineWidth', 2); %Ca
        plot(t,y(:, 2), 'r', 'LineWidth', 2); %Cb
    end
    %графики
    axis([0 20 0 62]);
    
    Cb_max = max(y(:,2));
    Cb_hl = Cb_max/2;
end

function f = f(t,y)
    global s S D v e;

    Ca=y(1); Cb=y(2);
    f = zeros(2,1);
    
    f(1)= -s/S*D*(Ca - Cb) - v*Ca;
    f(2)=  s/S*D*(Ca - Cb) - e*Cb;
end