Найти радиус описанной окружности по трем сторонам. не могу понять где ошибка! - Pascal

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

program raduys;
function rad(var a,b,c:real):real;
var r,p,s:real;
begin
s:=a*b*c;
p:=(a+b+c)/2;
r:=s/sqrt(sqrt(p*(p-a)-(p-b)*(p-c)));
rad:=r;
end;
var
r,a,b,c:real;
 
begin
writeln('vvedite a=,b=,c,');readln(a,b,c);
writeln(r);readln;
end.

вводишь любые значения сторон, а результат все равно равен 0

Решение задачи: «Найти радиус описанной окружности по трем сторонам. не могу понять где ошибка!»

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

function rad(a,b,c:real):real;
 var r,p,s:real;
 begin
 s:=a*b*c;
 p:=(a+b+c)/2;
 r:=s/(4*sqrt(p*(p-a)*(p-b)*(p-c)));
 rad:=r;
 end;
 
 var
 r,a,b,c:real;
 
 begin
 writeln('vvedite a=,b=,c,');
 readln(a,b,c);
 writeln(rad(a,b,c));
 readln;
 end.

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

Функция rad(a,b,c) принимает три аргумента типа real (a, b и c) и возвращает значение типа real (радиус описанной окружности).
Внутри функции объявлены три переменные r, p и s типа real.
Значение переменной s вычисляется как произведение a, b и c.
Значение переменной p вычисляется как среднее арифметическое a, b и c.
Значение переменной r вычисляется как значение s, деленное на 4 и возведенное в квадрат, деленное на значение p и умноженное на (p-a), (p-b) и (p-c) соответственно.
Значение функции rad(a,b,c) присваивается переменной rad.
В основной программе объявлены четыре переменные r, a, b и c типа real.
В основной программе выводится приглашение для ввода значений a, b и c.
В основной программе вызывается функция rad(a,b,c) с введенными значениями a, b и c.
В основной программе выводится результат вызова функции rad(a,b,c).
В основной программе выводится приглашение для ввода дополнительной информации.
Программа завершается после чтения дополнительной информации.