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