added Niels' formula