program proef; function macht(x,p : double) : double; begin macht := exp(ln(x)*p); end; procedure vgl(x,y : double); var v,f : double; begin f := x/macht(4*(sqr(y)*y+1),1/3) + y/(x+1) + 1/(x+y) - sqrt(9/4+sqr(x-y)/(x*y+x+y)); v := x/macht(4*(sqr(y)*y+1),1/3); Write(100*(ln(v)-ln(f)):10:6,' % = '); Writeln(100*(1-f/v):10:6,' %'); end; procedure test(N : integer); var k,m : integer; begin m := 1; for k := 1 to N do begin m := m*10; Writeln('10^',k,' : '); vgl(m,1); vgl(m,2); Writeln; end; end; begin test(9); end.