Ошибка в программе решения (Pascal)

[Surik]

Приветствую всех
Дали задание в лабораторной работе по
Алгоритмизации и программированию
Нужно написать программу для этого примера:

program bla;
var x,y,z,a,b,c,d,e,f:real;
function Arcsin (x :real ) : real;
begin
if abs(x) = 1 then Arcsin:= 0
else  Arcsin:= Arctan(x/Sqrt(1-sqr(x)));
end;
begin
readln(x);
a:=(x*x*x-1)/(2*sqrt(3));
b:=exp(x*x+1);
c:=abs(x-1)/abs(x-2);
c:=ln(c);
y:=ln(2.3);
d:=(x-1)/2;
d:=Arcsin(d);
e:=sin(x/2);
e:=e*e*e;
f:=sqr(2*x/3);
f:=sin(f)/cos(f);
z:=a*b*c*y-d+e+f;
writeln(z:8:1);
readln;
end.

Компилируется нормально
но при вычислении "z"
TurboPascal ругается вот так:

Причем в двух местах
Первое место то что на скрине
А второе это

c:=ln(c);

Причем вторая ошибка не появляется
если вводимое число больше "3"
А Если ввести число "3" то вычисление проходит на ура!
Пожалуйста подскажите где я ошибся
Изучаю паскаль всего 3 дня по колледжной программе
До этого были два самостоятельных урока по дельфи
поэтому надеюсь на местных знатоков и профессионалов!

[JUDAS]

Проверил Ваш код и нашёл 3 уязвимости
1. b:=exp(x*x+1); - вылетает при больших Х
2. c:=ln(c); вылетает при Х=2
Кусок кода c:=abs(x-1)/abs(x-2); написан не правильно, нужно так:
c:=abs((x-1)/(x-2));
В формуле же написано - модуль от результата деления
3. d:=Arcsin(d); вылетает при если вводимое число больше "3"
тут ошибка понятна слону: d:=(x-1)/2; при Х>3 число D>1
при проверке if abs(x) = 1 then Arcsin:= 0 ошибка в условии, нужно так:
if abs(x) >= 1 then Arcsin:= 0
-------------------
Я тут "накидал" Ваш пример с условиями и проверками.

program Project1;
var z, x  : double;
    v1,v2,v3,v4,v5,v6 : double;
begin
 readln(x);
 v1 := (x*x*x-1)/2*sqrt(3);
 v1 := arctan(v1);
 v2 := exp(x*x+1);
 if (x=2) then
 begin
  writeln('bezkonechost');
  exit;
 end;
 v3 := ln(abs((x-1)/(x-2)))*ln(2.3);
 if abs(x)>3 then
 begin
  writeln('net resheniy');
  exit;
 end;
 v4 := (x-1)/2;
 if v4=1 then v4:=pi/2 else
 v4 := arctan(v4/Sqrt(1-sqr(v4)));

 v5:=sin(x/2);
 v5:=v5*v5*v5;

 v6:=sqr(2*x/3);
 v6:=sin(v6)/cos(v6);

 Z:=v1*v2*v3 - v4+v5+v6;
 writeln(z:8:1);
end.

Посмотрите, может поможет в обнаружении ошибок в Вашем коде

[Surik]

Спасибо большое
Сейчас буду пробывать