## Example 1:
## dX(t) = 2*(3-X(t)) *dt + dW(t)
f <- expression( 2*(3-x) )
g <- expression( 1 )
res <- rsde1d(drift=f,diffusion=g,M=100,N=1000,tau=0.5412)
kurtosis(res)
skewness(res)
bconfint(res,level = 0.95)
moment(res,order=c(2,3,4,5))
Run the code above in your browser using DataLab