function [X1,Y1] = e_inter(Y, q)
% E_INTER - interpolacja sygnału funkcją e^(-pi t^2)
% Wejście:
%   Y - wektor próbek
%   q - skala (np. dla skali 2 otrzymamy 2x więcej próbek
% Wyjście:
%   X1 - czas dla 'nowych próbek',
%   Y1 - interpolacja.

p = width(Y);
X = 1:p;

X1 = voffrept(X, q - 1);
Y1 = zeros(size(X1));
for i=1:1:width(X)
    Y1 = Y1 + Y(i) * exp(1) .^ (-1 * pi * (X1 - i).^ 2);
end
X1 = (X1 - 1) * q + 1;
end