Adaptive Quadraturwird mittels Simpson umgesetzt

restart: Digits:=8:Funktion eingebenf:=x->sin(x);

NiM+SSJmRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc2luR0YlNiM5JEYlRiVGJQ==

Fehlererr:=0.0001;

NiM+SSRlcnJHNiIkIiIiISIl

Grenzena:=0; b:=Pi;

NiM+SSJhRzYiIiIhNiM+SSJiRzYiSSNQaUdJKnByb3RlY3RlZEdGJw==

Berechnungh:=(a,b)->(b-a) / 2: S:=(a,b)->h(a,b) / 3 *(f(a) + 4*f(a + h(a,b)) + f(b)): AdaptQuad:=proc(a, b) local fehler, m; global err; fehler:=evalf(abs(S(a,b) - S(a,(a+b)/2) - S((a+b)/2, b))); if fehler < 15 * err then return S(a,b); else m:=(a+b)/2; err:=err/2; return AdaptQuad(a, m) + AdaptQuad(m, b); end if; end proc:Ausgabeevalf(AdaptQuad(a, b));

NiMkIikjcC0rIyEiKA==

restart: Digits:=8:Funktion eingebenf:=x->sin(x);

NiM+SSJmRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc2luR0YlNiM5JEYlRiVGJQ==

Anzahl Unterteilungenn:=2;

NiM+SSJuRzYiIiIj

Bereicha:=0; b:=Pi;

NiM+SSJhRzYiIiIhNiM+SSJiRzYiSSNQaUdJKnByb3RlY3RlZEdGJw==

Berechnungh:=(b-a) / n: S:=(a,b)->h/3 *(f(a) + 2* sum( f(a+2*j*h) , j=1.. (n/2)-1)+ 4* sum( f(a+(2*j-1)*h) , j=1.. (n/2)) + f(b)):Ausgabeevalf(S(a, b));

NiMkIil4SCw/ISIo