Aufgabe 4 a, b, c, d

a)
Plotten des absoluten Fehlers:

ya = []; yb = []; y_b = [];
for i = 1:20;
  ya = [ya, 0.5*log(10) - romberg('int1','trapez',0,3,10,i)]
  yb = [yb, 0.704037753083331094342 - romberg('int2','trapez',0,(pi/2),10,i)]
  y_b = [y_b, 0.704037753083331094342 - romberg('int2','trapez',(pi/2),0,10,i)]
end
semilogy(1:i, ya, 1:i, yb, 1:i, yt)

a), b)
Werte der Integrale:
octave:163> romberg('int1','trapez',0,3,10,20)
ans = 1.15129254649701
octave:164> romberg('int2','trapez',0,(pi/2),10,20)
ans = 0.704037752073156
octave:165> (-1) * romberg('int2','trapez',(pi/2),0,10,20)
ans = 0.704037752073034
octave:166>

c)
Steigung ist unendlich im Punkt 0, siehe Graph

d)
siehe Werte der Integrale, von pi/2 nach 0 anstelle von 0 nach pi/2,
bringt aber nichts.