my.Psi <-
function(x, my.pi){
# calculates Psi
exp(log(1-my.pi) + dnorm(x, 0, 1, log=T) - log(my.pi + (1 - my.pi) * pnorm(x, 0, 1) ))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.