Линейная интерполяция, программа на TP

[marishaliza]

Добрый день!

Пожалуйста, помогите с решением

имеется программа на ТР, которая выводит таблицу значений функции с рассчитанным шагом h на интервале [c, c+30h]
необходимо с помощью линейной интерполяции вычисляет значения функции в точках xi=c+i*h+(i mod 4+1)/5*h по таблице значений функции с шагом h и вывесли значения xi, приближенные и точные значения функции в точках xi (i = 0,1,...,29).

uses crt;
const c=9;
var h,x,xi:real;
var n,i,k:integer;
{vichislenie funczii}
function f(x:real):real;
begin
f:=2*c*c*c*sin(x/c);
end;
begin
clrscr;
writeln('vvedite shag h');
readln(h);
writeln;
writeln('tablica znachenija funkcii na intervale [',c,', ',(c+30*h):3:2,']');
x:=c;
while x<c+30*h do
begin
write(f(x):9:4,' ');
x:=x+h;
end;
writeln;
writeln;
writeln('znachenie funkcii v tochkah xi');
For i:=0 to 29 do
begin
xi:=c+i*h+(i mod 4+1)/5*h;
write('x':2,i,'=',xi:5:3,' f(x',i,')=',f(xi):8:3);

Коммент педагога: Не найдены приближенные значения f(xi) с помощью линейной интерполяции.


end;
writeln;
readln;
end.


подскажите, пожалуйста, как вычислить приближенные значения f(xi) с помощью линейной интерполяции, какая формула?
в условиях задачи вот так написано:
Рассчитать шаг таблицы значений функции f(x), по которой с помощью линейной интерполяции можно было бы найти промежуточные значения функции с точностью 0.0001, если табличные значения функции округлены до 4-х знаков после запятой.