psip.bt <-
function(x,c1,M)
{
x1 <- (x-M)/c1
ivec1 <- (x1 < 0)
ivec2 <- (x1 > 1)
return(ivec1+(1-ivec1-ivec2)*((1-x1^2)^2+4*x*x1*(1-x1^2)/c1))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.