CONST
 Npoints=12;
VAR
 X,DX,S:Real;
 i:Integer;

Function Fun(xf:Real):Real;
 begin
  Fun:=SQR(Ln(xf))/xf;
 end;

BEGIN
 DX:=3/Npoints;
 S:=(Fun(1)+Fun(4))/2;
 For i:=2 to Npoints do
  begin
   X:=1.0+Dx*(i-1);
   S:=S+Fun(X);
  end;
 WriteLn('Approximate value = ',S*Dx:9:7);
 WriteLn('Accurate value = 0.8880657');
 Readln;
END.