Оценка оптимальной ставки налога на прибыль

#Определение оптимальной налоговой ставки
#Двухфакторный дискретный эксперимент, интервал исследования Т=10 
using Plots, PlotThemes

Ct=zeros(11)
Cpt=zeros(11)
Bd = zeros(11)
Bdsum=zeros(4,10)
K0=100  # исходный капитал
i=0
for Ren in (0.2:0.2:0.8) 
    j=0
    i=i+1
    for Snal in (0.1:0.1:1.0)
        j=j+1
        
       
        Cpt[1]=K0
        
        for t in (2:11)
           Ct[t]=Cpt[t-1]*Ren
           Cpt[t] =Ct[t]*(1-Snal)+Cpt[t-1]
           Bd[t] =Bd[t-1]+Ct[t]*Snal
           
        end
        Bdsum[i,j]=Bd[11]
    end
end

Ren = (0.2:0.2:0.8);
renlab1 = "Ren= "*string(Ren[1])
renlab2 = "Ren= " *string(Ren[2])
renlab3 = "Ren= " *string(Ren[3])
renlab4 = "Ren= " *string(Ren[4])

Bdsum

#plotly()
theme(:dao)
Snal = (0.1:0.1:1.0)
plot(Snal, Bdsum[1, :], label=renlab1, legend=:topright, wsize=(800, 500), linewidth=3)
plot!(Snal, Bdsum[2, :], label=renlab2, linewidth=3)
plot!(Snal, Bdsum[3, :], label=renlab3, linewidth=3)
plot!(Snal, Bdsum[4, :], label=renlab4, linewidth=3)
ylabel!("Поступление в бюджет")
xlabel!("Налоговая ставка")
title!("Определение оптимальной налоговой ставки (T=10)")