Значение определенного интеграла

[Aerowalk]

Здравствуйте, уважаемые программисты. Мне достался код программы (вычисление интеграла используя метод двойного пересчета по формуле средних прямоугольников), по идее он должен работать, но что-то не то. Я слабоват в паскале, так что просмотрите, пожалуйста, прогу и найдите ошибки. Сам паскаль вот на это ругается:
write('Znachenie integral = ', d_a(a,b,f));

PHP код:

program calc_int;

 type func=function (x:real):real;
 var a,b: real;
 
 function f(x:real):real;
  begin f:=Exp((-0.5)*Ln(3*x*x-1));
  end;

  
 {function v(x:real):real;}
 {begin v:=}
 {end;}

 
 function d_a(a,b:real; f:func):real;
  const N0=10; eps=0.001;
  var new,old,delta:real;
   N1:integer;

    function integral(a,b:real; N2:integer; f:func):real;
      var x,h,sum:real;
      i:integer;
      begin h:=(b-a)/N2;
        x:=a-h/2;
        for i:=1 to N2 do begin x:=x+h;
        sum:=sum+f(x) end;
        integral:=sum*h;
      end;

 begin {body of d_a}
  N1:=N0;
  old:=integral(a,b,N1,f);
  repeat
   N1:=2*N1;
   new:=integral(a,b,N1,f);
   delta:=abs(new-old);
   old:=new;
  until delta<eps;
  d_a:=new;
 end; {of d_a}

begin {main}
 write('a,b = ');
 read (a,b);
 write('Znachenie integral = ', d_a(a,b,f));
readln
end. 


Алгоритм правильный, если что.

[Vanta11a]

Вытащите integral из d_a и посмотрите, что получится.

[Aerowalk]

я что-то недопонял, а зачем вытаскивать integral ведь все логично там... меня волнует другой вопрос, а как вообще в теле правильно вызвать функцию d_a? может не надо в write писать? а как тогда написать? почему паскалю не нравится write('', d_a(a,b,f)); вообще правильно ли так вызывать?

a,b - граничные значения интеграла, f-сама подынтегральная функция

[Vanta11a]

Так. F должна вызываться со значением, т.к. в коде прописано

function f(x:real):real;

Т.е. надо как-то так:

write('Znachenie integral = ', d_a(a,b,f(c)));

Просто привычка, что все процедуры расположены отдельно.

И вообще, зачем Вы вызываете с f()?

[Aerowalk]

О точно насчет f(x), вместо c x должен стоять же, но почему-то паскалю это не нравится все равно. Да и сказали мне так делать, до этой программы я вообще думал что можно только write('Znachenie integral = ', просто переменная); Вы, Vanta11a, можете исправить тело (если хотите то и функции вначале растусовать (алгоритмы правильные, не вникайте)) программы так, как вам видится правильнее, у меня у самого что-то не получается, сил нет видеть эти ошибки при запуске

А а можно вызывать без f()? паскаль ругается на это

может в паскале дело, он 7.1 версии?

[Aerowalk]

аа все разобрался спасибо Vanta11a