diff --git a/TP1b/DiogoEliseuHugo_TP2_10.m b/TP1b/DiogoEliseuHugo_TP2_10.m
new file mode 100644
index 0000000..665a657
--- /dev/null
+++ b/TP1b/DiogoEliseuHugo_TP2_10.m
@@ -0,0 +1,60 @@
+%% Inicialização do ambiente
+clear ; close all; clc
+
+%% Exercício 10
+n = 0:250; % janela
+
+w1 = 0.1*pi;
+w2 = 0.2*pi;
+w3 = 0.3*pi;
+
+x1 = 2*cos(w1*n);
+x2 = 2*cos(w2*n);
+x3 = 2*cos(w3*n);
+
+% Figura com os resultados de x1
+figure
+hold on
+% Plot de x1
+subplot(2,1,1)
+stem(n, x1)
+title('x1')
+% Plot do output para x1
+subplot(2,1,2)
+stem(n, fun_y(x1))
+title('output para x1')
+hold off
+
+% Figura com os resultados de x2
+figure
+hold on
+% Plot de x2
+subplot(2,1,1)
+stem(n, x2)
+title('x2')
+% Plot do output para x1
+subplot(2,1,2)
+stem(n, fun_y(x2))
+title('output para x2')
+hold off
+
+% Figura com os resultados de x3
+figure
+hold on
+% Plot de x3
+subplot(2,1,1)
+stem(n, x3)
+title('x3')
+% Plot do output para x3
+subplot(2,1,2)
+stem(n, fun_y(x3))
+title('output para x1')
+hold off
+
+% Definição do sistema y
+function y = fun_y(N)
+    x0  = N;  %% For x[n]
+    x1  = [N(1, 2:size(N, 2)) 0];  %% For x[n-1]
+    xn1 = [0 N(1, 2:size(N, 2))]; %% For x[n+1]
+    y = x0.^2 - x1 .* xn1;
+end