#
# File Read Set  fixation2.vary2.set
# Init Conds  2 Par Range
#  change the :1 in Crv(1) Array(2) to :2
# screen will turn red/orange. 
# when done click on 
# nUmerics Color Code  Another quan
# Color via: rho
# Click on Optimize
# click [Esc] 
# click Restore
#  will see rho vs w 9<w<.3 and N < 200
#
par TT=5,RR=3,PP=1,SS=0;
par m=10
!a=m*RR
!b=SS+(m-1)*PP
!c=TT+(m-1)*PP
!d=m*PP
f(i,w,N)=1-w+w*(a*(i-1)+b*(N-i))/(N-1)
g(i,w,N)=1-w+w*(c*i+d*(N-i-1))/(N-1)
gam(i,w,N)=g(i,w,N)/f(i,w,N)
init p=1,s=1
@ meth=discrete,total=1000
par w,nn=1000
aux ww=w
p'=p*gam(t+1,w,nn)
s'=s+p
z'=z+1/nn
init z=0
@ bound=10000
aux rho=nn/s
done