Octave / task_10_r4

function [x,y] = task_10_r4 (f, xs, y0)
x = xs(1) : 0.1 : xs(2);
h = 0.1;
y = zeros([1 size(x,2)]);
y(1) = y0;
for i = 1 : size(y,2) - 1
k1 = h * f(x(i),y(i));
k2 = h * f(x(i) + h/2, y(i) + k1/2);
k3 = h * f(x(i) + h/2, y(i) + k2/2);
k4 = h * f(x(i) + h, y(i) + k3);
y(i+1) = y(i) + 1/6 * (k1 + 2*k2 + 2*k3 + k4);
endfor
endfunction