Переробити код Матлаб

В мене є код для ітераційного процесу синхроної апроксимаціїї, потрібно змінити його для асинхронної 

ось це код для цієї формули

function rez = cursova()

    format longG;

    clear;

    e = 1e-6;

%---------------------------------------------------------------------------------- 

     x_k = [0;0.5];

     disp('---------Метод з синхронною апроксимацією--------')

     count = 0;

     F = [f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

         f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)];

     A=inv(F);

     x_next=x_k-(A*[f1(x_k(1),x_k(2), 1); f2(x_k(1),x_k(2), 1)]);

     A_next=A*([2*[1 0;0 1]-[f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

         f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)]*A]);

     count = count+1;

     disp(x_k)

     disp(abs(f1(x_next(1),x_next(2),1)))

     disp(abs(f2(x_next(1),x_next(2),1)))

     while (norm(x_next-x_k,"inf") >=e)

         A=A_next;

         x_k=x_next;

         x_next=x_k-(A*[f1(x_k(1),x_k(2), 1); f2(x_k(1),x_k(2), 1)]);

         A_next=A*([2*[1 0; 0 1]-[f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

                f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)]*A]);

         count = count+1;

         disp(x_k)

         disp(abs(f1(x_next(1),x_next(2),1)))

         disp(abs(f2(x_next(1),x_next(2),1)))

     end

    disp('К-ть ітерацій : ')

    disp(count)

    disp('Приблизний результат')

    disp(x_k)

    disp('f1(a)')

    f1_x=f1(x_next(1),x_next(2),1);

    disp(abs(f1_x))

    disp('f2(a)')

    f2_x=f2(x_next(1),x_next(2),1);

    disp(abs(f2_x))

%-----------------------------------------------------------------------------

% перше р-ня системи

    function rez = f1(x, y, fun)

if fun == 1

    %rez = cos(1/3*(x-y))-2*y;

     rez = 2*x-sin(0.5*x-0.5*y);

else

  %rez = cos(1/3*(x-y))-2*y;

     rez = 2*x-sin(0.5*x-0.5*y);

end

end

% друге р-ня системи

function rez = f2(x, y, fun)

if fun == 1

   %rez = sin(1/3*(x+y))-2*x;

    rez = 2*y-cos(0.5*x+0.5*y);

else

   % rez = sin(1/3*(x+y))-2*x;

   rez = 2*y-cos(0.5*x+0.5*y);

end

end

% похідна 1го р-ня по X

function rez =f1dx1(x, y, fun)

if fun == 1

    %rez = -1/3*sin((x-y)/3);

     rez = 2-(0.5*cos(0.5*x-0.5*y));

else

   % rez = -1/3*sin((x-y)/3);

     rez = 2-(0.5*cos(0.5*x-0.5*y));

  end

end

% похідна 1го р-ня по Y

function rez =f1dx2(x, y, fun)

   if fun == 1

     %  rez = 1/3*(sin((x-y)/3)-6);

        rez = 0.5*cos(0.5*x-0.5*y);

   else

     % rez = 1/3*(sin((x-y)/3)-6);

      rez = 0.5*cos(0.5*x-0.5*y);

   end

end

% похідна 2го р-ня по X

function rez =f2dx1(x, y, fun)

    if fun == 1

      % rez = 1/3*(cos((x+y)/3)-6);

         rez = 0.5*sin(0.5*x+0.5*y);

    else

      % rez = 1/3*(cos((x+y)/3)-6);

         rez = 0.5*sin(0.5*x+0.5*y);

    end

end

% похідна 2го р-ня по Y

function rez =f2dx2(x, y, fun) 

   if fun == 1

       % rez = 1/3*cos((x+y)/3);

       rez = 0.5*sin(0.5*x+0.5*y)+2;

   else

       %rez = 1/3*cos((x+y)/3);

       rez = 0.5*sin(0.5*x+0.5*y)+2;

   end

end

end

ПОТРВБНО ПЕРЕПИСАТИ ДЛЯ ЦІЄЇ