st_prima[i,j] = st_prima[i,j-1]*exp((r/360-rf/360-(volatilidad^2)/2)*dt+volatilidad*sqrt(dt)*rnorm(1)) compensacionesCall[i] <- max(st_prima[i,T+1]-k2,0)*exp(-(r-rf ...