powered by
rybifr(x, y, np, nthet = np + 1, itype = 2, icoll = 0, isigma = 1, ch = 1.345, ck = 1.05 * sqrt(nthet), bm = 1.05 * sqrt(nthet), tol = 0.001, tau = 1e-06, maxitt = 50, maxitw = 80)