function [y,e] = sigmadelta(x) N = length(x); y = zeros(1,N); a = y; b = y; c = y; d = y; e = y; for i = 1:1:N %for i = 2:1:N iprev = max(1,i-1); % obchodzi iterowanie dla ujemnych 'i' dla poczatkowych indeksow d(i) = y(iprev); %d(i) = y(i-1); a(i) = x(i) - d(i); c(i) = b(iprev); %c(i) = b(i-1); b(i) = c(i) + a(i); if b(i)>0 y(i) = 1; else y(i) = -1; end e(i) = y(i)-b(i); end end