Интерполяционная формула Ньютона. (Pascal)

Листинг программы
Program IntNewton;

Uses
Crt;

Type
matr = array [1..20, 1..20] of real;
matr2 = array [1..20] of real;

Function fact(z: longint): longint;
..
Результаты выполнения

Please enter a grade of interval(a:b)
0
0.5
...

Цель: Используя интерполяционную формулу, найти значение функции в точке, не являющейся интерполяционной и оценить погрешность вычислений; подобрать аналитическую формулу, представить с некоторой точностью данные табличного значения функции.

Задание: Вариант3. Функция cos(x^2). Интервал: (0; 0,5).

Вывод. Используя интерполяционную формулу, нашел значение функции в точке, не являющейся интерполяционной и оценил погрешность вычислений; подобрал аналитическую формулу, представил с некоторой точностью данные табличного значения функции.

p.s. В работе также содежится исходный код, блок схема, результаты выполнения и теоретические сведения