Fixt bug u = h/ hu now is the right way u = hu/h