levin [ar,sigma2,rc]=lev(r)