R/derpsiBY3.R

derpsiBY3 <-
function(t,c3)
{
res=NULL

  for (i in 1:length(t))

{
if (t[i] <= c3)

 { res=rbind(res,0) }

else

{res=rbind(res,-exp(-sqrt(t[i]))/(2*sqrt(t[i]))) }

}
res

}
musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.