function [Y_int] = integ(Y, dx)
%Интегрирует поточечно, возвращает массив
Y_int = zeros(size(Y));
Y_int(1,:) = Y(1,:)*dx;
for i = 2:max(size(Y))
Y_int(i,:) = Y_int(i-1,:) + Y(i,:).*dx;
end

end