Напечатать таблицу значений функции - Pascal ABC

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

Напечатать таблицу значений функции f=(3x^2-13x+5)^(1/2) в интервале [-5;5] с шагом 0.5 с использованием подпрограмм. Как это сделать? Что только не пробовал вообще не как не получается. Заранее спасибо.

Решение задачи: «Напечатать таблицу значений функции»

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

program pr;
procedure ca;
var
  f, x: real;
begin
  x := -5;
  repeat 
    f := exp((1 / 2) * ln(3 * x * x - 13 * x + 5));
    writeln(x:4:1, ' | ', f:5:2);
    x := x + 0.5;
  until x > 5;
end;
begin
  writeln(' x   |   y ');
  writeln('------------');
  ca;
end.

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

В программе объявлены две переменные: f и x, оба типа real.
В процедуре ca инициализируется значение переменной x равное -5.
Затем, в цикле repeat, переменная f вычисляет значение функции exp(ln(3xx-13*x+5)) и записывает результат в переменную y.
После каждого вычисления значение переменной x увеличивается на 0.5.
Цикл продолжается до тех пор, пока значение переменной x не станет больше 5.
В конце программы выводится сообщение « x | y », где x - это текущее значение переменной x, а y - результат вычисления функции exp(ln(3xx-13*x+5)).
Выводится также сообщение «------------» для разделения результатов.