function wynik = rprogresywna (x, h)
  wynik = (exp((x+h).^2+1) - exp((x).^2+1))/2*h;
endfunction