function taud=root(a,b,tau)
    taud = tau*nthroot(a,b)
end