Производные паскаль

[Igomax]

Нужно вычислить 1 и2 производные следующих функций -4*ctg^3(x) и (e^x-e^-x)/(e^x+e^-x)

Вот что у меня получилось

program gra;
const
 epsilon=0.0001;
type func=function(x:real):real;
var z,z2,y,y2,x:real;   {$f+}
function f1(x:real):real;
begin
f1:=-4*(cos(x)/sin(x))*(cos(x)/sin(x))*(cos(x)/sin(x));  тут пишет ошибку division by zero
end;
function f2(x:real):real;
begin
 f2:=(exp(x)-exp(-x))/(exp(x)+exp(-x));
end;

function fdx(f:func; a:real; eps:real):real;
var x1,x2,y1,y2:real;
    i,p:longint;
begin
 i:=1;
 p:=-1;
 x1:=a+5;
  y1:=f(x1);
  repeat
   y2:=y1;
   x2:=x1;
   p:=p*(-1);
   x1:=(p/i)+a;
   y1:=(f(x2)-f(x1))/(x2-x1);
   i:=i+1;
  until abs(y1-y2)<eps;
  fdx:=y2;
 end;

function fdx2(f:func; a:real; eps:real):real;
var x1,x2,y1,y2:real;
    i,p:longint;
begin
 i:=1;
 p:=-1;
 x1:=a+5;
 y1:=f(x1);
   repeat
   y2:=y1;
   x2:=f(x1)/x1;
   p:=p*(-1);
   x1:=(p/i)+a;
   y1:=(fdx(f,x2,eps)-fdx(f,x1,eps))/(x2-x1);
   i:=i+1;
  until abs(y1-y2)<eps;
  fdx2:=y2;
 end;
begin
z:=fdx(f1,1,epsilon);
y:=fdx(f2,5,epsilon);
z2:=fdx2(f1,1,epsilon);
y2:=fdx2(f2,5,epsilon);
write (z:10:5, z2:10:5, y:10:5, y2:10:5);
end.

Как её исправить?

[russian-stalker]

if sin(x)<>0 then
f1:=-4*(cos(x)/sin(x))*(cos(x)/sin(x))*(cos(x)/sin(x));

[Igomax]

Сообщение от russian-stalker

if sin(x)<>0 then
f1:=-4*(cos(x)/sin(x))*(cos(x)/sin(x))*(cos(x)/sin(x));

От этого программа зависает и всё

[russian-stalker]

Ну ошибку то исправляет, если зависает, значит условие цикла не выполняется, проверьте ещё раз циклы

[Igomax]

Не могу найти что неправильно

[raxp]

деление на ноль тут ошибка...

[netrino]

Сообщение от Igomax

until abs(y1-y2)<eps

Здесь не ошиблись? repeat until <условие> - цикл, который можно прочесть как "делать, пока <условие> не будет выполнено".
То есть
repeat
i := i + 1;
until i = 10;
Будет увеличивать i до тех пор, пока оно не станет равным 10

[Igomax]

Сообщение от netrino

Здесь не ошиблись? repeat until <условие> - цикл, который можно прочесть как "делать, пока <условие> не будет выполнено".
То есть
repeat
i := i + 1;
until i = 10;
Будет увеличивать i до тех пор, пока оно не станет равным 10

Тут нет ошибка т.к. мне нужно добиться точности вычисления 0.00001

[Igomax]

Сообщение от raxp

деление на ноль тут ошибка...

Это понятно, но что надо изменить чтобы её не было?

[russian-stalker]

Я во втором посту написал как избавится от ошибки деления на ноль,