> n<-c(1,2,1,1,3)
> nzp.simple(n,20)
 [1] 1.276966 1.153788 1.095297 1.064895 1.048378 1.039191 1.034015 1.031077
 [9] 1.029404 1.028448 1.027902 1.027589 1.027410 1.027308 1.027249 1.027215
[17] 1.027196 1.027185 1.027179 1.027175
> nzp.newton(n,20)
 [1] 0.7787627 0.9714403 1.0244661 1.0271638 1.0271701 1.0271701 1.0271701
 [8] 1.0271701 1.0271701 1.0271701 1.0271701 1.0271701 1.0271701 1.0271701
[15] 1.0271701 1.0271701 1.0271701 1.0271701 1.0271701 1.0271701
> nzp.mos(n,20)
 [1] 1.057311 1.027271 1.027170 1.027170 1.027170 1.027170 1.027170 1.027170
 [9] 1.027170 1.027170 1.027170 1.027170 1.027170 1.027170 1.027170 1.027170
[17] 1.027170 1.027170 1.027170 1.027170
> 
> n<-c(1,2,1,1,1)
> nzp.simple(n,20)
 [1] 0.8385669 0.6812044 0.5927914 0.5366615 0.4983636 0.4709712 0.4507253
 [8] 0.4354010 0.4235938 0.4143723 0.4070941 0.4013021 0.3966627 0.3929270
[15] 0.3899063 0.3874557 0.3854619 0.3838363 0.3825085 0.3814223
> nzp.newton(n,20)
 [1]   -0.7903745   -2.5059333   -6.6397071  -13.5648162  -27.1321084
 [6]  -54.2642168 -108.5284336 -217.0568672 -434.1137344 -868.2274687
[11]          NaN          NaN          NaN          NaN          NaN
[16]          NaN          NaN          NaN          NaN          NaN
> nzp.mos(n,20)
 [1] 0.4513577 0.3772372 0.3764381 0.3764380 0.3764380 0.3764380 0.3764380
 [8] 0.3764380 0.3764380 0.3764380 0.3764380 0.3764380 0.3764380 0.3764380
[15] 0.3764380 0.3764380 0.3764380 0.3764380 0.3764380 0.3764380