<t>How can I simulate the path of a CIR process?I need to simulate the whole path,dr=k(h-r)dt+s*sqrt(r)dW(t)do i just simulate withr(t)=r(t-1)+k(h-r(t-1))*dt+s*sqrt(r(t-1))*I*sqrt(dt)where I is a N(0,1)..basically, i am doing that at the moment, but the volatility bit always drags r(t) into negative...