Introduction - If you have any usage issues, please Google them yourself
Given N non-equidistant points on the function value, calculated on a specified interval with the specified three interpolating polynomial interpolation points on the function value