function wynik = rcentralna (x, h)
  wynik = tan((x+h).^2+1) - tan((x-h).^2+1);
endfunction