Численый методы для курсовик

program Integral;
uses
  Crt, Dos;
var
  dx,x1,x2,e,i:real;
function Fx(x:real):real;
begin
  Fx:=x*x; 
end;
procedure CountViaTrap;
var
  xx1,xx2,xx3:real;
  c:longint;
begin
  writeln('------------------------------------------------');
  writeln('-->metod trapecii.');
  writeln('vsego iteracii:',round(abs(x2-x1)/e));
  i:=0;
  for c:=1 to round(abs(x2-x1)/e) do begin
    write('iteracia ',c,chr(13));
    xx1:=Fx(x1+c*e);
    xx2:=Fx(x1+c*e+e);
    if xx2>xx1 then xx3:=xx1 else xx3:=xx2;
    i:=i+abs(xx2-xx1)*e+abs(xx3)*e;
  end;
  writeln('------------------------------------------------');
  writeln('integral=',i);
end;
begin
clrscr;
begin
  writeln('------------------------------------------------');
  writeln('Vvedite ishodnie znachenia:');
  write('Nachalnoe znachenie x (x1)=');Readln(x1);
  write('Konechnoe znachenie x (x2)=');Readln(x2);
  write('Tochnost vichislenia (e)=');Readln(e);
 
end;
  CountViaTrap;
  writeln('------------------------------------------------');
  readkey;
end.

unit Dec_Integral;
 
interface
 
uses crt;
 
type
  TDec_Integral=object
  function Fun(x:real):real;
  function Met_Simpsona(a,b,e:real):real;
  procedure Integration;
  constructor Init;
  destructor Done;
end;
 
implementation
 
  function TDec_Integral.Fun(x:real):real;
  begin
    fun:=x*x;
  end;
 
  function TDec_Integral.Met_Simpsona(a,b,e:real):real;
  var n,i:integer;
      s_ab,h,s_even,s_odd,s_res,s_pres:real;
  begin
    n:=4;
    s_ab:=Fun(a)+Fun(b);
    h:=(b-a)/n;
    s_even:=0; s_res:=0;
 
    i:=1;
    while(i<n)do
    begin
    s_even:=s_even+Fun(a+i * h);
    i:=i+2;
    end;
    repeat
    begin
      s_odd:=0;
      s_pres:=s_res;
 
      i:=1;
      while(i<n) do
      begin
        s_odd:=Fun(a+i*h);
        i:=i+2;
      end;
 
      s_res:=h/3*(s_ab+2*s_even+4*s_odd);
 
      if (s_res=0) then s_res:=e;
      s_even:=s_even+s_odd;
      n:=n*2;
      h:=h/2;
    end; until (abs((s_pres-s_res)/s_res)>e);
 
    Met_Simpsona:=abs(s_res);
  end;
 
  procedure TDec_Integral.Integration;
  var a,b,e,result:real;
  begin
    clrscr;
 
    writeln('?*⥣аЁа®ў**ЁҐ:');
    write('‚ўҐ¤ЁвҐ a:'); readln(a);
    write('‚ўҐ¤ЁвҐ b:'); readln(b);
    write('‚ўҐ¤ЁвҐ в®з*®бвм:'); readln(e);
    result:=Met_Simpsona(a,b,e);
    writeln;
    writeln('3        ');
    writeln('ф        ');
    writeln('іXэdX=',result:4:4);
    writeln('х        ');
    writeln('1        ');
    readln;
  end;
 
  constructor TDec_Integral.Init;
  begin
 
  end;
 
  destructor TDec_Integral.Done;
  begin
    clrscr;
  end;
 
end.

unit Dec_Integral;
interface
uses
  Crt, Dos;
type TDec_Integral=object
function Fun(x:real):real;
function M_Trapecii(a,b,e:real):real;
procedure Integration;
constructor Init;
destructor Done;
end;
 
implementation
 
function TDec_Integral.Fun(x:real):real;
begin
fun:=x*x;
end;
 
function TDec_Integral.M_trapecii(a,b,e:real):real;
var
  dx,i:real;
  xx1,xx2,xx3:real;
  c:longint;
begin
  writeln('------------------------------------------------');
 
  i:=0;
  for c:=1 to round(abs(b-a)/e) do begin
    write('iteracia ',c,chr(13));
    xx1:=fun(a+c*e);
    xx2:=fun(a+c*e+e);
    if xx2>xx1 then xx3:=xx1 else xx3:=xx2;
    i:=i+abs(xx2-xx1)*e+abs(xx3)*e;
 
  end;
   M_Trapecii:=abs(i);
end;
procedure TDec_Integral.Integration;
var a,b,e,result:real;
begin
clrscr;
 
 
  writeln('------------------------------------------------');
  writeln('Vvedite ishodnie znachenia:');
  write('Nachalnoe znachenie x (a)=');readln(a);
  write('Konechnoe znachenie x (b)=');readln(b);
  write('Tochnost vichislenia (e)=');readln(e);
 
  result:=M_Trapecii(a,b,e);
 
 
  writeln('------------------------------------------------');
  writeln('integral=',result);
  writeln('------------------------------------------------');
  writeln('-->metod trapecii.');
  writeln('vsego iteracii:',round(abs(b-a)/e));
  readln;
end;
 
constructor TDec_Integral.Init;
begin
end;
destructor TDec_Integral.Done;
begin
 clrscr;
end;
 
end.