betaf<-function(const, k, sigma, t){
gamma=sqrt(k^2+2*const*sigma^2)
betaf=-2*const*(1-exp(-gamma*t))/(2*gamma-(gamma-k)*(1-exp(-gamma*t)))
return(betaf)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.