This thesis covers detailed analysis of bond pricing function. It focuses on connections between mathematical definitions and financial practice and it points out advantages and drawbacks of currently used function. Well known properties of this function are extended to negative internal rate of return values. This topic is further discussed with internal rate of return polynomial equations solving. Taylor series approximation is also shown regarding duration and convexity of bonds.