T_dospeca=20/252*5
r_god<- 0.05
Akcije <- c(49,48.12,47.37,50.25,51.75,53.12,53,51.87,51.38,53,
49.88,48.5,49.88,50.37,52.13,51.88,52.87,54.87,54.62,55.87,57.25)
#r_ned<-r_god*5/252
Sigma=0.2
delta<-function(T,S,K,Sigma,r)
{
d1<-(log(S/K) + (r + Sigma^2/2)*T) / (Sigma*sqrt(T))

delta<- pnorm(d1)

return(delta)

}

Tdo=T_dospeca-(0:20)*5/252

Delta=delta(Tdo,Akcije,50,Sigma,r_god)
kupljene_akcije <-(Delta-c(0,Delta[1:20]))*100000
cena_kupljenih_akcija <- Akcije*kupljene_akcije

troskovi <-c()
troskovi[1] <- cena_kupljenih_akcija[1]
for(i in 2:21)
troskovi[i] <- troskovi[i-1]*(1+r_god*5/252)+cena_kupljenih_akcija[i]

delta_hedging <- cbind(Akcije, Delta, kupljene_akcije, cena_kupljenih_akcija,troskovi)
delta_hedging