помогите с паскалем. вычислить значения функции sin(x^2+1) с точностью до 0.001 на интервале от 0 до 1

Не достаточно данных!
Не сказано про шаг…
Пусть шаг табулирования 0,1.
Тогда так:

h := 0.1;
x := 0;
repeat
WriteLn('F(', x, ') = ', sin(x*x+1):0:3);
x := x + h;
until x > 1;