load C:\Users\123\fib.dat
plot(fib(:,1), fib(:,2), 'o')
x = fib(:,1);
y = fib(:,2);
%fit = polyfit(x, log(y), 1)
xlabel('Fibonacci Element Number [th]')
ylabel('Time [ms]')
title('Calculation of N-th fibonacci element')
%semilogy(x, y, 'o', x, exp(fit(2)).*exp(fit(1)*x))
%x = (0: 0.1: 5)';
%y = erf(x);
p=polyfit(x, y, 6);
f = polyval(p,x);
plot(x, y,'o',x,f,'-')
axis([0 50 0 7000])