g.r

g <- function(Z) {
    prix <- apply(grid,2,cumsum)
    prix <- grid
    for (i in 1:n) {
        prix[i,] <- S0*exp(sigma*prix[i,]+(r-sigma^2/2)*i/n*Time)
    }
    prix <- colMeans(prix) - K
    prix[prix<0] <- 0
    return(prix)
}