formula by akapitanov | SimScale