Cumulative distribution function of Vitality 2013 model
pvit13(x, par1, par2, par3, par4)
x
time
par1
r
par2
s
par3
lambda
par4
beta
cumulative probability