lab 5 / my_linfit3

function my_linfit3(fn,N)
if(exist(fn)~=2)
disp('file not found')
return
end%if
M=read_matrix(fn);
x=M(:,1);
y=M(:,N+1);
MY_Fun=fittype(@(a,b,c,x) -2-8*sin(x)+2*x.^3)
MY_res=fit(x,y,MY_Fun)

T=(x(1):0.1:x(end));
close all
hold on
plot(T,MY_res(T),'-',x,y,'o')
grid on
title('plot of model and original Data')
xlabel('X');
ylabel('Y');
legend({'apprix.function','Table'},'Location','best');
end