Zlozona_Kwadratura_Trapeza - Calka

%{
function Wynik = oblicz_zkt(f,n,a,b)
  h = (b - a)/n;
  s=0;
  for i = 1:(n-1)
    s+=f(a+i*h);
  endfor
wynik_zkt = (h/2)*(f(a)+f(b)+2*s);
blad_teoretyczny = (h^2)/12*(b-a);
Wynik = [wynik_zkt,blad_teoretyczny];
endfunction
%}

format long
function y=f(x);
  y=2/(1+x);
endfunction
a = 0;
b = 0.8;
n = 12;
%Wynik rzeczywisty calki obliczony recznie
wynik_rz = 1.17557332980423801;

%Wynik  obliczony obliczony za pomoca wzoru
wynik_zkt = oblicz_zkt(@f,n,a,b)
%wynik(1) = oblicz_zkt(@f,n,a,b)
%wynik(1)

%Wynik obliczony za pomoca funkcji octave quad
q_wynik = quad(@f,a,b)

%Blad Teortyczny z funkcji
blad_teoretyczny = oblicz_zkt(@f,n,a,b)
%wynik(2) = oblicz_zkt(@f,n,a,b)
%wynik(2)

%Poprownanie wyniku_zkt(obliczonym za pomoca wzoru) z wynikiem_q  (obliczonym za pomoca funkcji octave quad)
%Szacowanie bledu teoretycznego
porow_w_zkt__q_w = abs(wynik_zkt - q_wynik)

%Porownanie wyniku_rz (obliczonego recznie) z wynikiem_zkt (obliczonym za pomoca wzoru)
%Szacowanie bledu rzeczywistego
porow_w_rz__w_zkt = abs(wynik_rz - wynik_zkt)

%Porownanie wyniku_rz (obliczonego recznie) z wynikiem_q  (obliczonym za pomoca funkcji octave quad)
porow_w_rz__q_w = abs(wynik_rz - q_wynik)

clear;