powered by
dr = alpha(beta-r)dt + sigma sqrt(r) dW,
with market price of risk q = q1/sqrt(r) +q2 sqrt(r). The maturities are 1,3,6 and 12 months.
bond.cir(alpha, beta, sigma, q1, q2, r0, n, maturities, days = 360)
out = bond.cir(0.5,2.55,0.365,0.3,0,3.55,1080,c(1/12, 3/12, 6/12, 1),365)
Run the code above in your browser using DataLab