function wynik = rich()
  h = 0.1;
  x=1;
  A = zeros(3);

  for j=1:2;

    #metoda progresywna
    if j == 1
      for i=1:3;
        A(i,1)=(sin((x+h).^2)- sin(x.^2))./h;
        h=h/2;
      end;
      p=1;
     end;

     #metoda centralna
     if j == 2
       for i=1:3;
        A(i,1)=(sin((x+h).^2) - sin((x-h).^2))./2*h;
        h=h/4;
      end;
      p=2;
     end;

    A(2,2) = A(2,1) + (A(2,1) - A(1,1))/(2.^p -1);
    A(3,2) = A(3,1) + (A(3,1) - A(2,1))/(2.^p -1);
    A(3,3) = A(3,2) + (A(3,2) - A(2,2))/(2.^(2*p) -1);
    W(j)=A(3,3);
   end;

   display(W)
end;