disp('Wir behandlen jetzt gamma von 3 (es sollte 2 sein)');
disp('n         erg           diff');
klaus = [2 4 8 16]';
erg =[GaussQuad(inline('x^2'),2,2) GaussQuad(inline('x^2'),4,2) GaussQuad(inline('x^2'),8,2) GaussQuad(inline('x^2'),16,2)]';
diff=[abs(gamma(3)-erg(1)) abs(gamma(3)-erg(2)) abs(gamma(3)-erg(3)) abs(gamma(3)-erg(4))]';
schorsch = [klaus erg diff];
disp(num2str(schorsch));
disp('Und wir sind fertig');

disp('Wir behandlen jetzt gamma von 7 (es sollte 720 sein)');
disp('n      erg        diff');
klaus = [2 4 8 16]';
erg =[GaussQuad(inline('x^6'),2,2) GaussQuad(inline('x^6'),4,2) GaussQuad(inline('x^6'),8,2) GaussQuad(inline('x^6'),16,2)]';
diff=[abs(gamma(7)-erg(1)) abs(gamma(7)-erg(2)) abs(gamma(7)-erg(3)) abs(gamma(7)-erg(4))]';
schorsch = [klaus erg diff];
disp(num2str(schorsch));
disp('Und wir sind fertig');

disp('Wir behandlen jetzt gamma von 11 (es sollte 3628800 sein)');
disp('n      erg        diff');
klaus = [2 4 8 16]';
erg =[GaussQuad(inline('x^10'),2,2) GaussQuad(inline('x^10'),4,2) GaussQuad(inline('x^10'),8,2) GaussQuad(inline('x^10'),16,2)]';
diff=[abs(gamma(11)-erg(1)) abs(gamma(11)-erg(2)) abs(gamma(11)-erg(3)) abs(gamma(11)-erg(4))]';
schorsch = [klaus erg diff];
disp(num2str(schorsch));
disp('Und wir sind fertig');

disp('Wir behandlen jetzt gamma von 3/2 (es sollte 0,88622692554276 (sqrt(pi)/2) sein)');
disp('n      erg        diff');
klaus = [2 4 8 16]';
erg =[GaussQuad(inline('x^(1/2)'),2,2) GaussQuad(inline('x^(1/2)'),4,2) GaussQuad(inline('x^(1/2)'),8,2) GaussQuad(inline('x^(1/2)'),16,2)]';
diff=[abs(gamma(3/2)-erg(1)) abs(gamma(3/2)-erg(2)) abs(gamma(3/2)-erg(3)) abs(gamma(3/2)-erg(4))]';
schorsch = [klaus erg diff];
disp(num2str(schorsch));
disp('Und wir sind fertig');




































%copyright phx ltd