Octave / task_10_r3_system

function [x,y] = task_10_r3_system (f1, f2, xs, znach)
  x = xs(1) : 0.1 : xs(2);
  h = 0.1;
  u = zeros([1 size(x,2)]);
  u(1) = znach(1);
  v(1) = znach(2);
  for i = 1:size(x,2) - 1
    ku1 = h * f1(x(i),u(i),v(i));
    kv1 = h * f2(x(i),u(i),v(i));
    ku2 = h * f1(x(i) + h/2,u(i) + ku1/2, v(i) + kv1/2);
    kv2 = h * f2(x(i) + h/2,u(i) + ku1/2, v(i) + kv1/2);
    ku3 = h * f1(x(i) + h,u(i) + 2 * ku2 - ku1, v(i) + 2 * kv2 - kv1);
    kv3 = h * f2(x(i) + h,u(i) + 2 * ku2 - ku1, v(i) + 2 * kv2 - kv1);
    u(i + 1) = u(i) + 1/6 * (ku1 + 4 * ku2 + ku3);
    v(i + 1) = v(i) + 1/6 * (kv1 + 4 * kv2 + kv3);
  endfor
  y = [u;v]';
endfunction