Const
 a=1.0;
 b=4.0;
 N=40;

Var
 AV:Real;

Function F(x:Real):Real;
begin
 F:=x*x*x/(x*x*x*x+1);
end;

Function Antiderivative(x:Real):Real;
begin
 Antiderivative:=ln(x*x*x*x+1)/4;
end;

Function Trapezium:Real;
var
 i:Byte;
 Sum,h:Real;
begin
 h:=(b-a)/N;
 Sum:=(F(a)+F(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F(a+h*i);
 Trapezium:=Sum;
end;


Begin
 Writeln('Result: ',Trapezium:0:8);
 AV:=Antiderivative(b)-Antiderivative(a);
 WriteLn('Accurate value: ', AV:0:8);
 Readln
End.

Const
 x=1.75;
 e=0.001;

Var
 S,a:Real;
 n:Integer;

Begin
 a:=1;
 S:=a;
 n:=0;
 repeat
  n:=n+1;
  a:=a*Sqr(x)/((2*n-1)*(2*n));
  S:=S+a;
 until a<e;
 Writeln('Result: ',S:0:8);
 WriteLn('Accurate value: ', (exp(x)+exp(-x))/2:0:8);
 Readln
End.

Const
 a=2.0;

Var
 x:Real;
 i:byte;

function L_n(t:real):real;
var
 n:byte;
 c,S:real;
begin
 n:=0;
 c:=(t-1)/(t+1);
 S:=c;
 for n:=1 to 10 do
  begin
   c:=c*Sqr(t-1)/Sqr(t+1)*(2*n-1)/(2*n+1);
   S:=S+c;
  end;
 L_n:=S*2;
end;

function Z1(a1,x1:real):real;
begin
 if x1>=a1 then Z1:=a1*L_n(x1) else Z1:=x1*L_n(a1);
end;

function Z2(a2,x2:real):real;
begin
 if x2>=a2 then Z2:=a2*Ln(x2) else Z2:=x2*Ln(a2);
end;

Begin
 Writeln(' x         Result        Accurate');
 for i:=0 to 6 do
  begin
   x:=3.0+0.5*i;
   writeln(x:3:1,Z1(a,x):15:5,Z2(a,x):15:5);
  end;
 Readln
End.