erho.bt.lim.p <-
function(p,c1)
# derivative of erho.bt.lim wrt c1
return(chi.int.p(p,2,c1)+c1^2*chi.int2.p(p,0,c1)+2*c1*chi.int2(p,0,c1))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.