function cal = calka;
F = input ('\n\nPodaj wzor funkcji:\nf(x) = ','s'); 
fn = inline(F); %Przedstawienie pobranego ciągu znaków z F jako funkcji (x)
a = input ('\nPodaj poczatek przedzialu calkowania, a = '); 
b = input('\nPodaj koniec przedzialu calkowania, b = '); 
while a > b 
    fprintf ('\nPodaj poprawne przedzialy:\n');
    a = input ('\nPodaj poczatek przedzialu calkowania a = ');
    b = input('\nPodaj koniec przedzialu calkowania b = ');
end; 
n = input  ('\nPodaj ilosc przedzialow, n= ');
h = (b-a)/(2*n); 

pole=0; %wartosc calki
z=0; 
y=0;
for m = 1:n %m to jest licznik
    i=2*m
    x  = a + i*h 
    z  = z + fn(x - h)
    if i<(2*n) 
    y = y + fn(x)
    end;
end;

pole=(h/3)*(fn(a) + fn(b) + 4*z + 2*y)
fprintf('\n\n\nWarosc calki w podanym przedziale = %0.4f',pole'); 
%podanie wyniku z dokładnością do 4 miejsc po przecinku