Вычисление второй производной функции

Если вычислить производную от производной, то получается вторая производная функции. Эта программа вычисляет вторую производную функции с заданной точностью (можете менять). Кроме того в программе можно поменять и функцию и точку в которой производная вычисляется.

program proizvodnaya;
 uses crt;
 const epsmin=1e-6;
type
 func=function(x:real):real;
var
 L1,L2,L3,epsilon:real; i:integer;
 {$F+}
 function ftest(x:real):real;
 begin
 ftest:=sin(x)
 end;
 function df(f:func; a,eps:real):real;
 var pme:integer; y1, y2, dx:real;
 begin
 pme:=-1;y2:=(f(a+5)-f(a))/5;
 i:=1;
 repeat
 y1:=y2;
 pme:=pme*(-1);
 dx:=pme/i;
 y2:=(f(a+dx)-f(a))/dx;
 i:=i+1;
 until (abs (y2-y1) < eps ) ;
 df:=y2;
 end;
Begin
 clrscr;
 {Производная в точке pi/2}
 L1:=df(ftest, pi/2, epsmin);
 {Производная в точке (pi/2+приращение)}
 L2:=df(ftest, pi/2+0.01, epsmin);
 {Отношение приращения первой производной к приращению аргумента}
 L3:=(L2-L1)/0.01;
 writeln('Вторая производная sin(x) при x=pi/2 ',L3:8:4);
 readln;
End.