%% Sygnały
x=[1 2 -1 3 -5]
y=[-2 2 10 -15 32 -41 26 -15]
%% Elementy
N=numel(x)
M=numel(y)
H = M-N+1
%% Równanie
 h=0*(1:H);
for n=1:H
 if n==1
    h(n)=y(1)/x(1);
 else
 rownanie_wew=0;
 for k=1:n-1
 rownanie_wew = rownanie_wew + h(k)*x(n+1-k);
 end
 h(n)=(1/x(1))*(y(n)-rownanie_wew)
 end
end