Mathematica routine package calculating shape of bended rod under the force load