lnHeyde<-function(a,x){
ff<-(1/x/((2*pi)^0.5))*exp((-(log(x))^2)/2)*(1+a*sin(2*pi*log(x)))}
par(mfrow=c(2,1))
plot(function(x)lnHeyde(0,x),0.01,2,ylim=c(0,1),main = "Lognormalvariante: Dichten mit gleicher MGF (a=-0.5,-0.1,0,0.1,0.5)",ylab="")
plot(function(x)lnHeyde(0.1,x),0.01,2,add=TRUE)
plot(function(x)lnHeyde(-0.1,x),0.01,2,add=TRUE)
plot(function(x)lnHeyde(0.5,x),0.01,2,add=TRUE)
plot(function(x)lnHeyde(-0.5,x),0.01,2,add=TRUE)
plot(function(x)lnHeyde(0,x),0.01,5,ylim=c(0,1),ylab="")
plot(function(x)lnHeyde(0.1,x),0.01,5,add=TRUE)
plot(function(x)lnHeyde(-0.1,x),0.01,5,add=TRUE)
plot(function(x)lnHeyde(0.5,x),0.01,5,add=TRUE)
plot(function(x)lnHeyde(-0.5,x),0.01,5,add=TRUE)