S-Lindley~CDF
pSL = function(x,p)
{
p = 0.22660
cos((pi/2)*(1+((x*p)/(1+p)))*exp(-p*x))
}
ks1=ad1=cvm1=NULL
ks1=ks.test(y,pSL)
ad1=ad.test(y,pSL)
cvm1=cvm.test(y,pSL)
result1=c(ks1$statistic,ks1$p.value,ad1$statistic,ad1$p.value,
```

```
cvm1$statistic,cvm1$p.value)