Hi,
I'm trying to do something very simple. I just want to write down a standard bond pricing model. Specifically, suppose my model block is
r=mu*r(-1) = eps_r;
p = exp((1+r)*eps_r(+1) + (1+r)^2)
where eps_r has sd=1.
We know that p=1 always since eps_r(+1) is normal. But when I do an approximation, no matter how many orders, I find a nonzero and potentially large coefficient on r (it must be zero). is there any way to get around this problem?
Thanks in advance!