Программа вычисления площади фигуры, ограниченной кривой

    Program dkr3; 
    uses crt;
    var n,num:integer; 
    a,b,h,x,xb,s,s1,pg:real;
    
    
    function f(x: real):real;
    begin
     f:=(x*x*x)-4*x; //сама функция
    end;
    
    
    function fb(x:real):real;
    begin
      fb:=((x*x*x*x)/4)-(2*(x*x)); //производная
    end;
    
    
    procedure proc1;
      begin
    write('Введите нижний предел интегрирования '); 
    readln(a);
    write('Введите верхний предел интегрирования ');
    readln(b);
    write('Введите количество отрезков ');
    readln(n);
    h:=(b-a)/n; //дальше по методу правых прямоугольников
    s:=0; 
    xb:=a;
    for var i:=1 to n do begin 
    x:=xb+i*h; 
    s:=s+fb(x)*h;
    end;
    s1:=fb(b)-fb(a);
    pg:=abs(s1-s);
    writeln('Интеграл равен ',s:12:10); 
    writeln('Погрешность равна ', pg)
    end;
    
    
    
    begin
      writeln('Реализовать программу вычисления площади фигуры, ограниченной кривой 1*x^3+(0)*x^2+(-4)*x+(0) и осью Ox');
      repeat
      writeln('1 - узнать производную этой функции');
      writeln('2 - приближенное значение');
      writeln('3 - выход из программы');
      readln(num);
      case num of
        1:writeln('(x^4/4)-(2x^2)');
        2:proc1;
        3:halt;
      end;
      until num=3;
    end.

    s:=s+fb(x)*h;

Program dkr3;
 
uses
  crt; 
 
function f(x: real):real := x * (x * x - 4); //сама функция
    
function fb(x: real):real := x * x * (x * x / 4 - 2); //производная
 
procedure proc1;
  begin
    var a := ReadLnReal('Введите нижний предел интегрирования:'); 
    var b := ReadLnReal('Введите верхний предел интегрирования:');
    var n := ReadLnInteger('Введите количество отрезков:');
    var h := (b - a) / n; //дальше по методу правых прямоугольников
    var s := 0.0; 
    for var i := 1 to n do s += f(a + i * h) * h;
    var s1 := fb(b) - fb(a);
    var pg := abs(s1 - s);
    writeln('Приближённое значение интеграла равно ', s); 
    writeln('Погрешность равна ', pg)
  end;
   
begin
  writeln('Программа вычисления площади фигуры, ограниченной кривой'); 
  writeln('1*x^3+(0)*x^2+(-4)*x+(0) и осью Ox'); 
  writeln('1 - узнать первообразную этой функции'); 
  writeln('2 - приближенное значение интегала'); 
  writeln('3 - выход из программы'); 
  var num: char; 
  repeat
    num := ReadKey;
    GoToXY(1, 6);
    for var i := 7 to 18 do writeln(' ' * 79);
    gotoXY(1, 7);
    case num of
      '1': write('x^4/4-2x^2');
      '2': proc1
    end;
    GoToXY(1, 6);
  until num = '3'
end.