Вычислить значение выражения, учитывая ОДЗ - Turbo Pascal

Формулировка задачи:

Даны числа x, y, z. Вычислить N, если N=(8y)^(1/3)-pi*ln(arctan(1/(zx))

Решение задачи: «Вычислить значение выражения, учитывая ОДЗ»

Листинг программы

var x,y,z,n,s:real;
begin
writeln('Введите значения x,y,z');
readln(x,y,z);
if(z*x=0)or(arctan(1/z/x)<=0) then
writeln('Значение N не определено')
else
 begin
  if y=0 then s:=0
  else if y>0 then s:=exp(ln(8*y)/3)
  else s:=-exp(ln(-8*y)/3);
  n:=s-pi*ln(arctan(1/z/x));
  writeln('N=',n:0:4);
 end;
readln;
end.

Объяснение кода листинга программы

Создается переменная x типа real, и ей присваивается значение, введенное пользователем.
Создается переменная y типа real, и ей присваивается значение, введенное пользователем.
Создается переменная z типа real, и ей присваивается значение, введенное пользователем.
Создается переменная n типа real.
Если z*x равно нулю или arctan(1/z/x) меньше или равно нулю, то выводится сообщение Значение N не определено.
Иначе, если y равно нулю, переменной s присваивается значение 0.
Если y больше нуля, переменной s присваивается значение exp(ln(8*y)/3).
Если y меньше нуля, переменной s присваивается значение -exp(ln(-8*y)/3).
Переменной n присваивается значение s-pi*ln(arctan(1/z/x)).
Выводится значение переменной n с точностью до четырех знаков после запятой.
Программа завершается вызовом функции readln().