N=100;jmax=100;delt=.05;sig=.12;mu=.05-((sig^2)/2);ylim0=0:1*3;x0=0:N
ctr1<-ctr2<-0
par(new=T)
for(j in 1:jmax){
w=2;for(i in x0[-1])w=c(w,w[i]+mu*delt+sig*sqrt(delt)*rnorm(1))
if(w[N+1]>2.5)ctr1<-ctr1+1
if(max(w)>2.5)ctr2=ctr2+1
par(new=T)
plot(x0*delt*2,w,type='l',ylim=ylim0,col=j)
}
abline(2.5, 0)
ctr1
ctr2

dev.off()

pdf('question4d.pdf')