x <- c(1,2,3,4,5,6,7,8,9)

mean(x)

x<-4
dnorm(x, 0, 1)

x<- ((-4000):4000)/1000
y<-dnorm(x, 0, 1)
plot(x,y)

s<-2
x<- (0:4000)/1000
y<-dweibull(x, s, 1)
plot(x,y)

x<- (0:4000)/1000
y<-dgamma(x, 3, 3)
plot(x,y)

x<- (0:4000)/1000
yw<-dweibull(x, 2, 1)
yg<-dgamma(x, 3, 3)
plot(x,yw)
lines(x,yg,type = "l", col = "red", lwd=5)

1-pweibull(356,10,500)