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