lg(x_):=log(x_)/log(10);