inst/book/ex2.07.R

# ex2.07.R
# Cox-Ingersoll-Ross (CIR-2)
d <- expression( 6-3*x ) 
s <- expression( 2*sqrt(x) ) 
s.x <- expression( 1/sqrt(x) )
require(sde)
set.seed(123)
sde.sim(X0=10, drift=d, sigma=s, sigma.x=s.x, 
    method="milstein") -> X
plot(X,main="Cox-Ingersoll-Ross")

Try the sde package in your browser

Any scripts or data that you put into this service are public.

sde documentation built on Sept. 9, 2022, 3:07 p.m.