# Luento 17.11.05
# gamma
# x^(a-1)exp(-x)

x<-seq(0, 20, 0.01)
a<-2.5
plot(x,x^(a-1)*exp(-x),type="line")
gamma(a)


x<-seq(0, 20, 0.01)
a<-5.5
plot(x,x^(a-1)*exp(-x),type="line")
gamma(a)

# Jakaumien kuvaajia
x<-seq(-3, 3, 0.01)
plot(x,dnorm(x),type="line")

mean<-0; sd<-sqrt(2)
plot(x,dnorm(x,mean,sd),type="line")