% badamy quad w prostym przypadku

quad('sin', 0, pi)
quad('sin', 0, pi)-2
quad('sin',0,pi,1e-10)-2
quad('sin',0,pi,1e-10)-2
quad('sin',0,11)
quad('sin',0,11) - (1-cos(11))

% a tu cos nie gra...

quad('sin',0,1111) - (1-cos(1111))
quad('sin',0,1111)

% kwadratura prostokatow
prostokaty('sin',0,1,100)- (1-cos(1))
% 10x wiecej wezlow zmniejsza blad... ile razy?
prostokaty('sin',0,1,1000)- (1-cos(1))

% kwadratura prostokatow kontra quad
quad('sin',0,1)- (1-cos(1))

% wyznaczamy wartość funkcji Gamma wprost kwadraturą
global x;
x = 4;

% jak wygląda funkcja podcałkowa?
t = linspace(0,1111,10000);
plot(t,podgamma(t))
semilogy(t,podgamma(t))

x
quad('podgamma',0,100)-6
x = 4.4;
quad('podgamma',0,100)
quad('podgamma',0,100)-gamma(x)

% czy to było trudne?
[I, nF] = quad('podgamma',0,100, 1e-3)

% a teraz Gauss-Laguerre...
% ...węzły...
tt = [
0.263560,
1.413403,
3.596426,
7.085810,
12.640801];
% ...wagi...
ww = [
0.521756,
0.398667,
0.759424e-1,
0.361176e-2,
0.233700e-4];
% kwadratura: tylko 5 wartości funkcji!!!
I = sum( (tt.^(x-1)).*ww)-gamma(x)

% ale: czy to jest poprawna kwadratura??...
% jak to sprawdzić?
