fs = function(x,epsilon,delta) dnorm(sinh(delta*asinh(x)-epsilon))*delta*cosh(delta*asinh(x)-epsilon)/sqrt(1+x^2)
vec = seq(-5,5,0.001)
plot(vec,fs(vec,0,1),type="l",ylim=c(0,0.5), xlab="Rendement espéré",ylab="Densité de probabilité")
points(vec,fs(vec,-0.5,1.1),type="l",col="red")