
这个是流程图。结果要模拟K。图中R<0时K=K+1
参数λ=0.01,c=1,β=0.001,u=1000,t0=180,要求循环1000次
我写的长这样:
workfile my u 1 10
!N=1000
series u=1000
series tt= 180
series m=1
series β=1/1000
series λ=0.01
series K=0
for !i=1 to 1000
series t=@rpoisson(λ)
series R=u
if t<=tt then
R=R+m*t-@rexp(β)
t=t+@rpoisson(λ)
if R<0 then K=K+1
endif
endif
next
show K
运行不能啊。。结果明显不对劲。。求大神助攻!