symulacje

##########################################################################################
###                                                                     Ocena prawdopodobieństwa zdarzenia rzadkiego                                                           ###
##########################################################################################

N=10^3
y=rexp(N)+4.8
waga=dnorm(y)/dexp(y-4.8)
sum(waga)/N
### Graficzna prezentacja wyników
plot(cumsum(waga)/1:N, type='l')
abline(a=pnorm(-4.8),b=0, col='red')