v0 = [ 1; 0; 0 ];
sigma = 10;
R = 28;
b = 8/3;
T = 1;
f = @(v) [ -sigma*(v(1)-v(2)); R*v(1)-v(2)-v(1)*v(3); v(1)*v(2)-b*v(3) ];

len = 10;

result = zeros(3,len+1);
h = zeros(1, len-1);
e = zeros(1, len-1);

for i = 0:len
	h(i+1) = 10^-2*2^-i;
	V = euler(f, v0, h(i+1), T);
	result(:,i+1) = V(:,end);
end

% calculation of error:
for i = 0:len
	e(i+1) = norm(result(:,i+1)-result(:,len+1));
end

loglog(h, e, '.-')