niter = 1e5
proba = rep(0,niter)
probb = rep(0,niter)
for (i in 1:niter)
{
r = rnorm(5,.001,.015)
logPrice = log(1000) + cumsum(r)
proba[i] = as.numeric(logPrice[1]<log(990))
probb[i] = as.numeric(logPrice[5]<log(990))
}
mean(proba)
## [1] 0.23088
Also the same problem(a) can be solved as P(x< 990). Given, the daily log returns as normal iid. so, log(Pf/P0) as x log(990/1000) = -0.01005034 P(x<-0.01005034) = cdf(-.01005034) = pnorm(-0.01005034,.001,.015)
pnorm(log(990/1000),.001,.015)
## [1] 0.2306557
mean(probb)
## [1] 0.3273