#pexp(2,1/3)
x<-seq(0,20,length=1000)
dens<-dexp(x,1)
distr<-pexp(x,1)
barva<-rainbow(10)

x<-seq(0,6,length=1000)
plot(x,dexp(x,dens),type='n',lwd=2,ylab='hustota',
     main='Exponencialni rozlozeni - hustota')
lines(x,dexp(x,1/2),col=barva[2],lwd=2)
lines(x,dexp(x,1),col='black',lwd=2)
lines(x,dexp(x,2),col=barva[4],lwd=2)
lines(x,dexp(x,4),col=barva[6],lwd=2)
lines(x,dexp(x,6),col=barva[8],lwd=2)
lines(x,dexp(x,8),col=barva[1],lwd=2)


plot(x,distr,type='n',lwd=2,
     ylab='pravdepodobnost',main='Exponencialni rozlozeni - distr.fce')
lines(x,pexp(x,1/2),col=barva[2],lwd=2)
lines(x,pexp(x,1),col='black',lwd=2)
lines(x,pexp(x,2),col=barva[4],lwd=2)
lines(x,pexp(x,4),col=barva[6],lwd=2)
lines(x,pexp(x,6),col=barva[8],lwd=2)
lines(x,pexp(x,8),col=barva[1],lwd=2)


#-----------------------------------------------------------------------
#-----------------------------------------------------------------------
#normalni rozlozeni
x<-seq(-30,30,length=1000)
dens<-dnorm(x,0,1/2)
distr<-pnorm(x,0,1/2)
barva<-rainbow(10)

par(mfrow=c(1,2))
plot(x,dens,type='n',col='red',lwd=2,ylab='hustota',
     main='Normalni rozlozeni - hustota')
lines(x,dnorm(x,0,1/2),col=barva[1],lwd=2)
lines(x,dnorm(x,0,1),col='black',lwd=2)
lines(x,dnorm(x,0,2),col=barva[2],lwd=2)
lines(x,dnorm(x,0,4),col=barva[4],lwd=2)
lines(x,dnorm(x,0,8),col=barva[8],lwd=2)
legend(10,0.8,fill=c(barva[1],'black',barva[2],barva[4],barva[8]),
       legend=c('N(0;1/2)','N(0;1)','N(0;4)','N(0;16)','N(0,64)'))

plot(x,distr,type='n',lwd=2,
     ylab='pravdepodobnost',main='Normalni rozlozeni - distr.fce')
lines(x,pnorm(x,0,1/2),col=barva[1],lwd=2)
lines(x,pnorm(x,0,1),col='black',lwd=2)
lines(x,pnorm(x,0,2),col=barva[2],lwd=2)
lines(x,pnorm(x,0,4),col=barva[4],lwd=2)
lines(x,pnorm(x,0,8),col=barva[8],lwd=2)
legend(10,0.6,fill=c(barva[1],'black',barva[2],barva[4],barva[8]),
       legend=c('N(0;1/2)','N(0;1)','N(0;4)','N(0;16)','N(0,64)'))

#------------------------------------------------------------------------------------

#-----------------------------------------------------------------------
#normalni rozlozeni
x<-seq(-18,18,length=1000)
dens<-dnorm(x,0,1)
distr<-pnorm(x,0,1)

#par(mfrow=c(1,2))
plot(x,dens,type='n',col='red',lwd=2,ylab='hustota',
     main='Normalni rozlozeni - hustota')
lines(x,dnorm(x,-7,1),col=barva[1],lwd=2)
lines(x,dnorm(x,0,1),col='black',lwd=2)
lines(x,dnorm(x,7,1),col=barva[2],lwd=1)
legend(10,0.4,fill=c(barva[1],'black',barva[2]),
       legend=c('N(-7;1)','N(0;1)','N(7;1)'),cex=0.75)

plot(x,distr,type='n',lwd=2,
     ylab='pravdepodobnost',
     main='Normalni rozlozeni - distr.fce')
lines(x,pnorm(x,-7,1),col=barva[1],lwd=2)
lines(x,pnorm(x,0,1),col='black',lwd=2)
lines(x,pnorm(x,7,1),col=barva[2],lwd=1)
legend(10,0.4,fill=c(barva[1],'black',barva[2]),
       legend=c('N(-7;1)','N(0;1)','N(7;1)'),cex=0.75)