formula by xungb1n | SimScale