本スクリプトは代表的なワールの応力修正係数を採用しましたが、図からでも割り出すことができます。

ここでは次の計算式を採用してます。

κ = (4c-1)/(4c-4)+0.615/c

で割り出せます。