%Definition von b_1/b_2
b_2 = (0:0.001:1.6)'.*B.*rho; 
p_2 = b_2./B;
b_1 = NaN(size(b_2)); 
p_1 = NaN(size(p_2));
%definition von differenz
differenz = NaN(size(b_2));
% Schleife
for i = 1:1:size(b_2,1)
    %Bestimme Hashrate und dazugehörige Sicherheitsfunktion
    H_b = ((m-1)/m)*delta/rho*((rho-1)/(2*kappa))*b_2(i);
    if H_b>A
        S = 1-((A/H_b)^k);         
    else 
        S=0;
    end
    %Bestimme Konsummenge q
    q = delta/rho*b_2(i)/n;
    %Fallunterscheidung: Wenn q>1 ist f*(u_ableitung(q)-1)=0
    if q > 1
        %Rechne b1 aus
        b_1(i) = delta/rho*b_2(i)*S;
    else
        %Rechne b2 aus
        b_1 (i) = (delta/rho)*S*b_2(i)*(1+f*(u_ableitung(q)-1));
    end
    %differenz aus b_2 und b_1 bilden
    differenz(i)=b_1(i)-b_2(i);
end
    %bl bestimmen: VZW geht von negativ auf positiv 
    if differenz(i-1)<=0 && differenz(i+1)>=0
        b_l=b_2(i);
    
    %b2 bestimmen: VZW geht von positiv auf negativ
    elseif differenz(i-1)>=0 && differenz(i+1)<=0
        b_h=b_2(i);
    end